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The  articles  in  this  collection  contain  the  results  of 
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of  criteria!  formulas  for  calculating  the  heat-transfer 
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thermal  conductivity  are  solved  in  general  form  and  numerically 
for  cooled  hollow  blades  and  disks. 
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FOREWORD 

Recently,  ever  increasing  attention  has  been  paid  to  questions 
of  heated  and  loaded  parts  of  high-temperature  gas  turbines ,  since 
improvement  of  turbine  cooling  is  one  of  the  effective  ways  of  im¬ 
proving  their  performance.  The  materials  used  to  make  turbine  parts 
allow  raising  the  temperature  of  the  gas  ahead  of  a  cooled  turbine 
to  liJ00-1500°K.  However,  imperfections  of  the  cooling  systems  do 
not  allow  heating  the  gas  to  such  a  temperature. 

Turbine  cooling  can  be  improved  only  after  careful  study  of  the 
temperature  state  of  individual  stator  and  rotor  parts  under  operating 
conditions . 

In  this  collection  we  present  the  results  of  theoretical  and  ex¬ 
perimental  studies  of  tne  temperature  state  of  the  most  critical  parts 
of  turbine  stators  and  rotors  —  the  nozzle  and  rotor  blades  and  disks. 

The  statement  and  solution  of  research  problems  had  the  aim  of 
explaining  the  reasons  for  high  heating  of  these  parts  and  finding 
ways  of  Improving  methods  for  their  cooling. 

The  first  articles  analyze  the  boundary  conditions  in  the  rotors 
of  axial  gas  turbines  and  nozzle  blades.  The  rules  established  by  the 
authors  make  it  possible  to  refine  and  simplify  the  familiar  methods 
of  heat  calculations.  Then,  on  the  basis  of  a  generalization  of  the 
criterial  dependences,  used  to  determine  the  average  heat-transfer 
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coefficient  in  vane  cascades,  calculation  formulas  are  obtained  for 
calculating  this  coefficient  in  cascades  with  varying  reaction.  As 
a  result  of  an  examination  of  steady  and  transient  heat  conduction 
of  the  rotor  with  the  blades,  disks,  nozzle  blades,  and  rotor  blades 
and  other  parts,  finite  express.r.o*  s  are  derived  which  can  be  used  when 
analyzing  their  temperature  state. 

Subsequent  articles  study  the  Cooling  of  rotor  blades  by  blowing 
cooling  air  from  inside  onto  the  blade  surface  through  perforated 
holes  and  show  the  effectiveness  of  this  method,  and  also  give  a 
method  for  calculating  disk  cooling. 

Of  considerable  interest  are  the  results  of  an  experimental 
study  of  the  thermal  state  of  the  rotor  of  an  axial,  three-stage  gas 
turbine;  these  make  it  possible  to  evaluate  the  -ffectiveness  of  the 
cooling  system  and  note  ways  for  improving  it. 

The  last  article  presents  a  theoretical  study  of  the  process  of 
heat  transfer  in  the  flow  of  cooling  liquid,  considering  the  energy 
dissipation. 

The  results  of  the  studies  are  presented  in  the  form  of  mathe¬ 
matical  dependences  suitable  for  calculations;  the  most  complex 
solutions  are  presented  as  nomograms  or  are  reduced  to  numerical 
functions . 

The  authors  ask  that  all  comments  on  the  contents  of  this 
collection  be  sent  to  the  following  address:  Moscow,  B-66,  1-y 
Basmannyy  per.,  3,  Izdatel'stvo  "Mashinostroyeniye . " 
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BOUNDARY  CONDITIONS  IN  PROBLEMS  OF  KEAT  CONDUCTION  FOR  GAS-TURBINE 
ROTORS 


V.  S.  Petrovskiy  and  M.  I.  Tsaplin 


In  a  theoretical  study  of  the  temperature  states  of  gas-turbine 
rotors  the  boundary  conditions  are  formulated  on  the  basis  of  spe¬ 
cific  physical  premises  which  depend  on  the  structural,  gas-dynamic, 
and  operational  characteristics  of  the  turbine.  In  order  not  to 
complicate  the  solution  of  the  initial  differential  equations,  we 
usually  attempt  to  simplify  the  boundary  conditions  within  tolerable 
limits.  In  this  case  we  will  examine  only  one-dimensional  problems. 

Let  us  set  the  coordinate  system  as  shown  in  Fig.  1.  If  the 
heat  problem  is  solved  separately  for  a  rotor  blade,  at  the  base  of 
the  blade  (x  =  0)  the  boundax*y  condition  can  be 


MO)-/, 


C _ £ _ L, _ ^ 

l'-*'5  Joint.  *X 

1  digit _  I  t>  >«dt  •« 


Fig.  1.  Coordinate  system  used  for  stating 
heat  problems  for  the  cooled  rotors  of  gas 
turbines:  1  -  approximate  experimental 
curve;  2  —  theoretical  curve. 
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whdre  ).  is  tht  eat  conduction  of  the  blade  material. 
n 

In  both  cases  the  boundary  conditions  contain  values  which  can¬ 
not  be  determined  by  calculation,  and  therefore  a  theoretical  study 
of  the  temperature  field  is  possible  if  we  know  the  experimentally 
determined  values  of  temperature  t^  or  heat  flux  qfl  . 

Measurement  of  t, ,  to  say  nothing  of  q.  .  and  q*  involves 
laborious  experiment;  however,  this  is  Justified  since,  first,  by 
knowing  the  actual  temperature  on  the  boundary  of  the  blade  with  the 
disk  we  reduce  the  probability  of  error  when  calculating  the  tempera¬ 
ture  dependence  for  the  disk  or  blade,  while  information  on  heat  flows 

q  and  q'  makes  it  possible  to  predict  the  temperature  gradient 
•  A  A  •  A 

from  the  middle  to  the  root  of  the  blade  or  from  the  periphery  to  the 
center  of  the  disk.  Second,  measurement  of  t, ,  q  ,  and  q*  in 
conjunction  with  the  calculation  method  of  determining  the  tempera¬ 
ture  field  in  the  blade  and  disk  is  a  less  laborious  process  than 
purely  experimental  study  of  the  temperature  state  of  turbine  rotors. 

In  designing  a  turbine  there  is  no  need  to  measure  temperature 
t^,  and  therefore  we  assign  its  value  in  the  calculations.  As  has 
been  3hown  by  experimental  study  of  temperature  distribution  in  gas- 
turbine  rotor  blades  and  disks,  the  value  of  t-^  depends  to  a  con¬ 
siderable  extent  on  the  temperature  of  the  gas  at  the  blade  root  and 
on  the  cooling  of  its  root  section  or  footing  (if  such  exists).  If 
we  know  the  calculation  mean-mass  temperatures  of  the  gas  at  the 
inlet  to  the  rotor  cascade  and  at  its  outlet,  we  can  assume  that  t^ 
is  0.7-0. 8  of  the  arithmetic  mean  of  these  temperatures. 


We  can  avoid  the  need  for  measuring  t1,  and  q^  by 

solving  the  problem  for  a  turbine  rotor  for  the  system  disk/blade, 
i.e.,  by  solving  a  system  of  two  differential  heat-conduction  equa¬ 
tions.  In  this  case  we  must  formulate,  besides  the  boundary  con¬ 
ditions  for  x  *  l  and  r  *  0  (or  r  *  i’2)»  two  more  boundary  conditions 
on  the  blade/disk  joint. 
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In  final  analysis,  temperature  distribution  in  the  system  disk/ 
blades  is  caused  by  the  gas  flow  passing  through  the  rotor  cascade, 
and  by  disk  and  blade  cooling  conditions.  The  characteristics  of  the 
gas  flow  are  known  from  previous  thermal  gas-dynamics  calculations, 
while  the  characteristics  of  the  cooling  air  are  determined,  for  the 
given  design,  by  its  flow  rate.  The  flow  of  air  for  cooling  is 
usually  3-^%  of  its  total  flow  rate  through  the  engine.  Corresponding 
hydraulic  calculation  should  show  the  distribution  of  cooling  air  as 
it  move  through  the  entire  cooled  blading  of  the  rotor  before  passing 
into  the  gas  flow  or  to  the  atmosphere.  Knowing  the  distribution  of 
the  cooling  air,  we  can  calculate  its  speeds  of  motion  along  the 
cooled  surfaces,  and  then  the  values  of  the  heat-transfer  coefficient. 

Thus,  the  characteristics  which  determine  the  temperature  field 
in  the  rotor  can  be  obtained  from  calculation.  They  can  be  used  to 
determine  the  boundary  conditions  of  the  blade/disk  system  and  to 
find  the  temperature  field  as  functions  of  the  gas-flow  and  cooling 
air  parameters. 

First  let  us  examine  the  typical  boundary  conditions  for  turbine 
disks.  With  one-dimensional  statement  of  the  problem,  the  turbine 
disk  ir  a  symmetrical  body  whose  temperat>.-e  changes  along  the  radius. 
If  there  is  no  central  opening  (r  =  0)  the  heat  flow  in  the  center  of 
the  disk  (because  of  symmetry)  is  zero: 

dtA(0)/dr=0.  (3) 

For  a  disk  with  an  opening  in  the  center  (r  =  r2)  the  boundary 
condition  can  be  given  in  the  form  of  the  equation 

(*0 

where  r2  is  the  radius  of  the  opening  in  the  disk. 

Disks  have  variable  thickness  along  the  radius.  Their  profiles 
can  be  approximated  by  hyperbolic,  linear,  or  parabolic  functions 
(Fig.  2).  By  correcting  the  profile  or  heat-transfer  conditions  in 
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the  central  part  of  the  disk,  the  boundary 
conditions  (with  r  *  r2)  can  also  be 
written  in  the  form  [dt^(r2)/dr]  *  0. 

Here  the  solutions  of  the  initial  dif¬ 
ferential  equations  are  simplified. 


//"A  j  •'  r,  We  can  Sive  the  following  example  of 

^ 1  ;l'  Tj^^T —  correction.  If  the  disk  profile  can  be 

1  lz=i===sa|ij  *  approximated  by  a  hyperbolic  function, 

.  when  constructing  the  theoretical  profile 

a  disk  profile  by  linear,  there  can  possibly  be  too  strong  develop- 

hyperbolic,  and  parabolic  ment  of  its  lateral  surface  as  the  center 
functions : 

i  -  s  •  *0  ♦  Cl(r0  -  r)n,  n  »  lj  a  -  °f  the  disk  is  approached.'  Convection 

«  .  *0  *  t/rS  3  -  *  •  «/r.  heat  transfer  from  the  surface  of  the 

opening  will  be  compensated  by  this  in¬ 
creased  area  of  heat  transfer  with  corresponding  dimensions  of  the 
disk  opening.  For  disks  of  trapezoidal  profile  the  correction  can 
be  done  by  increasing  somewhat  the  heat-transfer  coefficient  on  the 
lateral  surface  of  the  disk. 


Fig.  2.  Approximation  of 
a  disk  profile  by  linear, 
hyperbolic,  and  parabolic 
functions : 

l  -  *  •  *0  ♦  -  «•)".  n  »  lj  a  - 

«  •  *0  ®  C/rs  3  -  *  *  «/*■ 


The  boundary  condition  for  a  disk  with  r 
the  form 

•.  _ • 


r1  can  be  written  in 


where  q'  is  the  heat  flow  to  the  disk  from  the  blades  and  the  vane 
A  •  A 

channel,  with  consideration  of  the  thermal  resistance  in  the  roots. 


The  value  of  q'  in  condition  (5)  is  determined  mainly  by  the 
heat  from  the  blades  to  the  disk.  But  this  heat,  in  turn,  is  deter¬ 
mined  by  the  boundary  condition  with  x  *  0,  which  can  be  found  from 
the  heat-balance  equation.  Consequently,  we  can  connect  the  tem¬ 
perature  fields  of  the  blade  and  disk  at  their  joint  using  the 
general  boundary  condition;  then  the  thermal  problem  for  the  general 
temperature  field  of  the  disk/blade  system  will  be  solved  using  only 
initial  data  from  gas-dynamics  and  aerodynamics  calculations  for  the 
gas  and  cooling  air,  respectively. 


Let  us  see  how  we  can  associate  functions  t  (r)  and  t  (x)  using 

A  n 

the  boundary  conditions. 

Let  us  assume  that  in  the  sector  from  t^(r^)  to  tfl(0)  (from  the 
disk  periphery  to  the  blade  root)  the  temperature  gradient  does  not 
change,  i.e.,  At/Ax  =  const.  This  makes  it  possible  to  represent 
the  temperature  differential  at  the  blade/disk  boundary  in  the  form 
of  the  difference 

Af*»ia(0)— (6) 

and  introduce  the  value  of  At  into  one  of  the  terms  of  the  heat- 
balance  equations  compiled  for  t^x),  with  x  *  0,  and  tfl(r),  with 
r  *  r, ,  since  the  conduction  heat  flow  from  the  blade  to  the  disk 
in  this  case  is  proportional  to  At. 

Thus,  if  we  solve  the  differential  equations  for  heat  conduction 
for  the  blade  and  the  disk  individually,  in  the  final  expressions  we 
get  the  unknown  At  which  we  introduced  into  the  boundary  conditions 
and  which  can  be  defined  by  supplementary  equation  (6).  Substituting 
into  it  the  values  of  t^(G)  and  ( r*^ )  a  found  from  the  general 

solutions  of  the  heat  conduction  equations,  we  get  a  formula  for 
calculating  At  in  terms  of  known  values. 

As  an  example,  let  us  examine  the  connection  in  terms  of  the 
temperature  differential  in  the  blade  and  the  disk.  Let  the  disk 
have  blades  with  roots,  and  also  an  upper  and  lower  shroud.  During 
assembly,  channels  are  formed  for  air  passage  during  turbine  operation 
between  the  surfaces  of  the  root,  the  lower  shroud,  and  the  disk 
crown.  The  diagram  of  such  a  blade,  mounted  on  a  disk,  is  shown  in 
Pig.  3- 

For  such  a  design  it  is  expedient  to  consider  the  change  in 
temperature  along  the  blade  root  from  t^(0)  to  t^(r1)  to  be  linear; 
the  coordinate  origin  along  the  blade  should  be  in  a  plane  which  coin¬ 
cides  approximately  with  the  middle  thickness  of  the  shroud,  while 
the  terminal  coordinate  of  the  radius  coincides  with  r^. 
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Pig.  3.  Schematic 
drawing  of  a  shrouded 
solid  rotor  blade 
with  root,  positioned 
on  a  turbine  disk: 

a  —  average  iiaat-tranafer 
coefficient*  of  the  corree- 
pondlng  surface* ,  Q  -  heat 
removed  froa  thea. 


is  the  heat  removed  by 


fer  surfaces  S'  and  S 
a5 


Considering  all  heat-transfer  surfaces 
when  x  *  0,  we  can  compile  the  following  heat- 
balance  equation  (see  Fig.  3): 

Qi+Qs^Qj+Qt*  (7) 

where  0.  =).  s,)i  is  the  amount  of  heat  fed 
per  unit  time  from  the  working  part  of  the 
blades  to  their  lower  sections  SQ,  which  we 
will  consider  to  be  located  in  plane  x  ■  0 
(n  is  the  number  of  blades)  j  (?*=»  —  u<t\ti{G)-~ir(0))S,,n 
is  the  heat  obtained  by  the  surfaces  of  the 
lower  shrouds  of  the  blades,  turned  toward  the 
gas  flow]  <?,=>.,  j-S*n  18  the  heat  removed  from 
boundary  sections  (SQ)  to  the  blade  roots  (h  is 
the  height  of  the  root,  SH  i3  the  area  of  its 
section) ;  <?<--- ^1 

the  cooling  air  from  the  convection  heat-trans- 

"  within  the  channels  in  the  lower  part  of  the 

“r, 


blades . 


Substituting  the  expressions  for  into  equation  (7)»  after 

transformations  we  get  (with  x  »  0) 


db,(  0) 
dx 


?  a +?  »• 


(8) 


Here 


M0)=/^0)-/a0): 
®3  ($.,  T^,J  *  «»$»,  . 


ri — 


■  Now  let  us  formulate  the  boundary  condition  for  the  disk  with 
r  -  r1.  The  heat-balance  equation 


Qi— 


(9) 


where  Q-« '/.3  —  SHn  is  the  quantity  of  heat  fed  per  unit  time  through 
the  root  cross  section  from  all  blades  to  the  disk  crown  (to  an  arbi¬ 
trary  coaxial  surface  of  radius  r  *  r^j  <!*(/,  (r, )—?„]£../*  is  the  heat 

removed  by  the  cooling  air  from  the  blade  surfaces  S"  of  convection 
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heat  transfer;  Q.=2nrt/yS/.1  is  the  heat  entering  the  disk  through 

dr 

the  outer  surface  of  the  crown  (H  is  the  width  of  the  disk  crown,  8 
is  a  coefficient  which  takes  into  account  the  thermal  resistance  ir. 
the  blade  joints). 


Substituting  into  equation  (9),  after  transformations  we 


j2iiOL_TA,r,i+?,. 

tlr 


Here 


?3—  — : 


'  •  >.a2 anH-i 

c  ^  1  Si/tllf 

™  *,2 * 

Boundary  conditions  (8)  and  (10),  together  with  boundary  con¬ 
ditions  (1)  and  (3)  or  (4),  make  it  possible  to  obtain  particular 

x  solutions  to  the  initial  differential  equations 

0f  for  the  blade  and  disk,  from  which  we  can  also 

find  the  formula  for  determining  At.  In  cal- 
^  x'l  culations  by  the  obtained  formulas  it  is  con- 

■fyc  •  ?  S  A  venient  first  to  find  the  numerical  value  of 

\w  (__5  'J  At  from  formula  (6),  and  then  calculate  the 

\  \  l\  temperature  distribution  in  the  blade  and  disk. 


vT  I  Let  us  examine  the  case  of  determination 

of  the  general  boundary  condition  at  the  joint 

Pig.  4.  Diagram  for  between  a  disk  and  blades  without  roots.  The 

compiling  heat-bal-  diagram  for  connecting  the  blade  to  the  disk 
ance  equations  for  ° 

lower  part  of  solid  in  this  version  is  shown  in  Pig.  The  origin 

shrouded  rotor  blade  f  coordinatu  x  (along  the  blade)  is  best  located 

without  root,  and  ° 

for  a  disk  with  r  =  in  the  root  section,  and  it  is  advisable  to  ex- 

—  v>  t 

amine  the  temperature  differential  throughout 
the  disk  crown  (coordinate  r).  Let  us  consider  that  plane  x  =  0 
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coincides  with  coordinate  r^.  For  blade  cross  section  SQ  coinciding 
with  this  plane  we  can  compile- the  following  heat-balance  equation 
for  the  amount  of  heat  passing  per  unit  time: 

Q»+Q#"Qio+Qu*  (11) 

Here  <?,=  /..  ihSSl s.jrt  is  the  heat  coming  to  the  cross  section  from 

dx 

the  working  parts  of  the  blades;  </,-=«,  f/^0)— is  the  heat  fed 

to  the  disk  crown  from  the  vane  channels;  Q10 * ,.x  is  the  heat 

In 

removed  from-  the  outer  surface  of  the  crown  (At  is  the  temperature 
differential  throughout  the  crown  from  r1  to  rp  ;  QncarSa*  /MJ  x 

x  the  heat  removed  from  the  side  surfaces  of  the  crown 

by  the  cooling  air. 

Substituting  Qg_11  into  (11),  after  transformations  we  get 


^-?A(0)+?« 


where 


^  tt|  'S.,  |  UQ  IXkrtfa  -T-hj,  2) 

'  ).j  SPn? 

_ 0«  4rl'f»(rl*^**«-2)  U  ,',n\  1 

- s7> - lWr(0)1“ 

_ |  ag_  2 1.1*  1 

Li,  S^n'i  So"?  I»  (»•,  r\)  J 

For  a  disk  with  r  *  r^  the  balance  equation  can  be  written  in  the  form 


Qi:— Qio* 


Here  Qu  =  i.t‘2xr,H — — —  is  the  heat  removed  per  unit  time  from  the  crown 
to  the  web  of  the  disk. 


After  substituting  Q12  and  Q]_q  lnt0  d3),  with  transformations, 
we  get 

Min 
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Where 


„ _ it 

ri  r'i  In  (rx<r\)  ' 

Obtained  boundary  conditions  (12)  and  (14)  include  the  tempera¬ 
ture  differential  At  *  tn(0)  -  tfl(r|).  Just  as  for  a  disk  having 
shrounded  blades,  we  can  first  find  the  particular  solutions  for  the 
blade  and  the  disk,  and  then  eliminate  the  unknown  At  from  them. 

The  boundary  condition  at  the  blade  tip  (x  *  l)  in  the  form 

(15) 

dx 

is  the  only  item  found  in  the  literature. 

Measurements  show  the  nonagreement  between  the  actual  temperature 
distribution  along  the  blade  and  that  calculated  with  boundary  con¬ 
dition  (15).  This  is  explained  by  the  fact  that  the  final  formula  for 
determining  t^(x)  does  not  take  into  account  changes  in  gas  tempera¬ 
ture  along  the  blade  (see,  e.g.,  [1,  2]).  In  addition,  boundary  con¬ 
dition  (15)  is  physically  incorrect,  since  there  is  actually  a  tem¬ 
perature  gradient  at  the  blade  tip. 

In  [2]  a  solution  is  giver  for  the  thermal  stationary  problem 
for  a  blade  with  a  profile  of  variable  cross  section,  in  which  the 
change  in  gas  temperature  along  the  blade  is  taken  into  account;  this 
same  source  recommends  finding  the  constants  in  the  general  solution 
of  the  initial  differential  equation  using,  in  particular,  boundary 
condition  (15).  The  function  for  the  change  in  gas  temperature  along 
the  blading  of  the  rotor  cascade,  leading  to  the  occurrence  of  a  tem¬ 
perature  gradient  at  the  blade  tip,  is  incommensurate  when  solving 
the  heat  problem  with  boundary  condition  (15).  To  ovei come  the  contra¬ 
diction  between  the  proposed  and  the  actual  boundary  conditions  at  the 
blade  tip  we  can  isolate,  at  the  blade  tip,  a  short  segment  of  length 
0.2 Z,  considering  its  temperature  constant  and  equal  to  the  temperature 
of  the  gas  on  the  periphery  of  the  blading  of  the  rotor  cascade.  Such 
an  assumption  is  not  far  wrong,  since  the  blade  tip  receives  a  cer¬ 
tain  amount  of  heat  from  its  middle  part  due  to  heat  conduction,  and 
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tihis  to  some  extent  compensates  for  heat  lost  to  radiation  in  the 
examined  segment.  For  rotor  blades  with  an  upp'er  shroud  we  can  dis¬ 
regard  the  heat  fed  to  the  shroud  through  the  blade  tip  section. 

Then  the  temperature  at  the  blade  tip  can  be  considered  equal  to  the 
temperature  of  the  shroud  and  defined  as  some  steady  temperature 
established  as  a  result  of  the  balance  of  the  heat  fed  to  the  shroud 
by  the  gas  flow  and  lost  by  it  due  to  radiation.  The  heat  lost  to 
radiation  can  be  defined  approximately  by  the  formula 

«•»■•••=»•  u6) 
where  e*  „  is  the  coefficient  of  blackness  of  the  shroud  surface:  a 
is  the  thermal-radiation  constant;  Tfi  is  the  steady  temperature  of 
the  shroud. 

When  compiling  equation  (16)  we  disregarded  the  radiation  heat 
transfer  between  the  gas  flow,  the  surrounding  parts,  and  the  shroud, 
considering  that  the  shroud  only  gives  off  heat.  This  assumption 
lowers  T,.  „  only  insignificantly. 

Based  on  conditions  of  convective  heating  of  the  shroud  and  the 
heat  lost  to  radiation,  we  can  write  the  heat-flow  equation  in  the 
form 

S|r;(/i-r..„]  ......  Ik.  (i7) 

Herp  a  =  0. 5(0^  +  «2)  is  the  average  heat-transler  coefficient  for 
the  upper  and  lower  shroud  surfaces;  T*(Z)  is  the  gas  temperature  at 
the  shroud. 

Numerical  solution  of  (17)  for  T*  can  be  easily  done  by  trial 

o  •  n 

and  error  with  subsequent  approximation. 


Thus,  for  the  tip  of  a  rotor  blade  of  a  gas  turbine  (x 
following  boundary  condition  is  more  precise: 

un-tt.u  or  w)=tL-tyi)< 


l )  the 
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Where  t«  •  is  the  steady  temperature  of  the  blade  tip,  equal  to  the 
>  o  •  n 

temperature  of  the  upper  shroud  6fl(Z)  a  t^U)  -  t*(Z.).  ' 


The  rotor  blacks  can  be  cooled  from  within  by  air,  liquid,  the 
vapors  from  some  liquid  or  other,  or  other  means;  .In  all  cases, 
when  formulating  the  boundary  conditions  for  x  .*»  0,  x  -  t,  r  t=  r^, 
r  *  rg»  and  r  *  r£,'  the  above  discussions  remain  valid.  Certain 
features  can  appear  only  when  compiling  the  equations  for  the  heat 
balance  at  the  blade/disk  Joint  and  when  compiling  the  heat-balance 
equation  for  the  determination  of  the  equilibrium  temperature  of 
the  rotor  blade  tip.  Thus,  boundary  conditions  (8),  (10),  (12),  (14), 
and  (18)  can  be  considered  as  typical  for  gas-turbine  rotor  blades. 

It  should  be  noted  that  for  rotta*  blades  with  an  upper  shroud, 
calculation  of  the  temperature  of  the^jlade  tip  using  boundary  con¬ 
dition  (15)  is  not  allowed,  since  an  erroneous  rise  in  temperature 
can  make  it  necessary  to  increase  the  area  of  the  profile  cross 
section  due  to  high  mechanical  loads  at  the  point  where  the  blade- 
touches  the  upper  shroud.  Increase  cross-sectional  area  increases 
the  dimensions  and  weight  of  the  blade  as  a  whole,  which  is  very 
undesirable  for  gas  turbines. 


One-dimensional  transient  heat  conduction  in  a  gas-turbine 
rotor  is  described  by  the  following  equations: 

for  the  disk  — 


Mr) 

dr  J 

C,0.;A  (r) 

for  solid  blades  — 

dx  S  (a)  dx  dx  J  0,0,5  (A)  1  •’ '  ’  1 


and  for  a  hollow  cooled  blade  — 


fX,(.x.T)  a,  d  fo, 

dx  ~  S(x)  dx  dA  j 


dA  j 


1 >■«  ^ 
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where  x  is  the  time;  z(r)  is  hair  the  disk  thickness j  a^,  a^,  p^, 

P/i*  Cfl*  and  c  i  are  thermal  conductivity,  density,  and.  specific 
heat  of  the  disk  and  blade,  respectively.  ° 

The  equations  which  describe  the  steady  heat  state  of  the  disk 
and  blade  can  be  obtained  from  equations  (19)-(21),  if  we  consider 
that  tfl  or  t^  do  not  depend  on  t,  i.e.,  3tfl(r,x)/3x  and  3t;)(x,t)/3T 
are  equal  to  zero. 

Solution  of  equations  (19)-(21)  in  general  form,  with  consider¬ 
ation  of  the  possible  dependences  of  z,  a'^,  cig,  a0,  S^,  P,  t»,  and  tQ 
on  x,  is,  as  a  rule,  impossible  due  to  insurmountable  mathematical 
difficulties..  In  this  case  the  simplicity  of  the  boundary  conditions 
allows  us  to  compile  an  initial  differential  equation  or  a  system  of 
equations  without  significant  simplifications,  and  thus  obtain  more 
precise  final  functions  tfl(r,x),  tfl(x,x),  or  tfl(r),  t^(x) . 

Valid  simplification  of  the  boundary  conditions  can  be  combined 
with  the  use  of  initial  differential  equations  in  the  most  general 
form,  if  we  consider  the  features  of  the  actual  temperature  distri¬ 
bution  along  the  blade  foil. 


Experiments  show  that  in  the  midsection  of  the  blade  the  tem¬ 
perature  distribution  curve  has  a  maximum.  This  is  natural,  since 
the  shape  of  the  dependence  t^Cx)  of  a  solid  or  hollow  cooled  rotor 
blade  is  determined,  to  a  considerable  extent,  by  the  temperature  of 
the  gas  in  the  blading.  Prom  this  it  follows  that  in  the  middle  of 
the  blade  foil  there  is  no  longitudinal  heat  flow.  Therefore,  if  we 
formulate  the  heat  problem  for  two  independent  sections  of  the  blade 
0  <  x1  <  and  0  <  x2  <  (assuming  that  w  l 2  and  +  l2  *  l) , 
on  the  boundary  of  sections  x-^  and  x2  (x1  =  x2  =  0)  we  can,  with 
sufficient  accuracy,  assume  that  3t(Z, ,x)/3x,  =  0  and  3t  (0,x)  ? 

J.  X  Jl  An 

t  3x2  =0. 


The  arb'  -ary  division  of  the  blade  into  two  sections  makes  it 
possible  to  i.  up  one  heat  problem  for  the  disk  and  the  root 


0,7  0,75  0,8 


Pig.  5,  Change  in  tem¬ 
perature  of  a  stagnant 
gas  flow  along  the  blade 
cascade,  in  relative 
coordinates : 

1  —  graph  constructed  from  measure¬ 
ment  results;  2,3-  graphs  of  the  type 

t*  •  1  -  CjXj.  whore  n  «  8  ar.d  6;  tf  * 

*  <x2)  •  t«(x2)/t»<0)x^;  5  -  graphs 

of  the  type  t*^)  «  1  -  C3O  -  Xj)", 


seccions  of  the  blades  (segment  x-^)  and 
another  independent  problem  for  the  blade 
tips  (segment  x2). 

Such  a  method  makes  it  relatively 
simple  to  take  into  account,  in  initial 
differential  equation  (20),  a  number  of 
functions  of  x,  including  the  most  im¬ 
portant  one  -  t#(x).  Figure  5  shows,  in' 
relative  coordinates,  the  most  typical 
change  in  gas  temperature  along  the  blade 
cascade.  Division  of  the  blade  into  seg¬ 
ments  x^  and  x2  makes  it  possible  to 
represent  function  t*(x)  In  the  form  of 
two  fourth-order  parabolas  on  segment  x^ 
and  two  eighth-order  parabolas  on  segment 


of  the  type  t •(x^)  •  1  -  1  -  *x) , 

whor«  n  -  3  and  H;  t»(x1>  -  tj(xx)  ♦  When  stating  transient  problems  for 

i  i  _■  z  gas  turbines  it  becomes  necessary,  as  a 

*  1-0U?.cP.nacc  <:iero  *i  •  W 

J2  -  x2/i2).  rule,  to  take  into  account  the  variable 

nature  of  the  process  not  only  along  the 
examined  direction  (with  one-dimensional  solutions),  but  also  on  the 
boundaries  of  the  corresponding  coordinates .  Let  us  examine  the 
possible  transient  boundary  conditions  and  formulate  them  in  more 
general  form. 


In  the  center  of  a  solid  disk  the  cona_tions  of  symmetry,  re¬ 
gardless  of  time,  allow  us  to  consider  that  3tfl(r,x)/3r  =  0.  This 
boundary  condition  is  not  always  acceptable.  If  the  central  part  of 
the  disk  is  in  contact  with  a  shaft  whose  temperature  is  always  kept 
lower  than  the  disk  temperature,  heat  will  be  removed  from  the  cen¬ 
tral  part  of  the  disk.  This  circumstance  must  be  considered,  even 
though  it  does  not  change  the  condition  3tfl/3r  =  0. 

In  this  case  we  can  state  the  problem  for  the  interval 
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if 


^2  t  r  1  ri»  where  r2  is  a  certa3;n  radius;  we  can  consider  approxi¬ 
mately  that  on  this  radius  the  direction  of  the  heat  flow  in  the 
disk,  from'  the  periphery  to  the  center,  coincides  with  the  direction 

of  reading.  The  value  of  the  heat  flow  q„  will  equal  Q/S„  (Q  is 

r2  r2 

the  total  quantity  of  heat  removed  to  the  turbine  shaft,  S  is  an 


arbitrary  coaxial  surface  on  the  disk  at  radius  r^) .  The  value  of 

q„  for  particular  cases  can  be  determined  from  equations  of  the 
r2 

conservation  of  energy  (heat  balance) 


For  a  disk  with  a  central  opening  the  inner  surface  of  the 
opening  can  be  cooled.  This  is  taken  into  account  by  the  boundary 
condition  (with  r  *  r2) 

Ur,  \ix (>.)—/>.)).  ( 22 ) 

In  the  absence  of  cooling  the  boundary  condition  will  be 

dl3(rj,x)ldr=0. 

At  the  point  where  the  blade  roots  join  the  disk  the  boundary 
conditions  consist  of  the  heat-balance  equation,  just  as  for  the 
case  of  steady  problems,  but  with  consideration  of  the  dependence 
of  At  on  time.  In  accordance  with  this,  stationary  boundary  con¬ 
ditions  (8)  and  (10)  assume  the  form  [3]: 


wnen  r  =  r. 


(0.  T)-ek,(r,.  xH-ri 


(23) 


and  when  x.^  =  0 


Here 


-- 


.  Sh»i  c* 

-  u,S,n 


>.,2  r\Hi 


(24) 


_•  •<.) « ,  _■  ,  _•  a^,  ■■"A, ,  s„  ,  <j 

?5=—  —  +  ?i.  ?•■»*=  (ri-riKlO:  =4  = - — — 1-~; 

A  /jS  o  ASy  2/.jSq 


*-  1% :  *-  °iSt  Cm,.  +SiS£&l_  ur,). 

/iSy  2 f.  ,S0  *  Wv  Mo 
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.  For  blades  without  roots  the  boundary  conditions  have  a  some¬ 
what  different  form.  It  is  not  advisable  to  use,  in  this  case,  the 
assumption  of  a  logarithmic  law  of  the  change  in  temperature  on  the 
disk  crown,  as  was  done  for  the  steady  process,  since  the  time  de¬ 
pendence  tfl(r,x)  is  of  a  different  nature.  Let  us  compile  a  heat- 
balance  equation  with  r  »  r1  and  x-^  =  0  (see  Fig.  5): 

<?;+<?;»  (25) 

where  Q s0n'}  is  the  heat  from  the  blades  to  the  disk  in 
1  second;  Q^  —  t^  -tx  (r}.  t)]  S,,ii  is  the  heat  to  the  disk  from  the 

vane  channels;  Q,  — Pi'JSLli  2 nrtH  is  the  quantity  of  heat  removed  from 

dr 

a  surface  of  radius  r^  to  within  the  disk  in  1  second. 


After  substituting  Q£,  Q^,  ar.d  into  equation  (25)  we  get 


0Mn.T)  „•  _x _ • 

— - ^ - *)=¥». 


or 


(26) 


where 


•  XrSonP  .  •  °|S«.n  ,  .  •  f  (t\\ 

•tl-T-Z - rr>  ?*=T" *T - —  •  — ri'r(0)r,. 

/  ,2 X,2nf|H 


Now  let  us  determine  the  boundary  condition  for  the  tip  of  the 
rotor  blade.  With  transient  heating,  the  temperature  of  the  blade 
tip  changes  somewhat  more  slowly  than  that  of  its  middle  section, 
due  to  temperature  nonuniformity  in  the  gao  flow.  First  let  us  ex¬ 
amine  a  shrouded  blade.  The  shroud  at  the  blade  tip  has  a 
sufficiently  large  neat-exchange  surface,  and  we  can  consider  that 
in  the  general  energy  balance  that  heat  which  is  fed  from  the 
middle,  more  heated,  part  of  the  blade  is  not  definitive.  Basically, 
the  temperature  of  the  shroud  depends  on  convection  heating  by  the 
gas,  and  to  a  lesser  extent  on  the  heat  lost  to  radiation  at  a 
high  temperature,  which  we  can  disregard.  Then  heating  of  the  shroud 
will  occur  in  accordance  with  a  regular  regime  of  the  first  kind, 
i .e. ,  when  x2  =  l 2 

(j  *)~  !^r(fj) — i$\  e  n  »  (27) 
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where  tfl(£2,i)  t^ie  'fcera'Perature  of  the  shroud  (and  the  blade  tip); 
tn  Is  the  initial  temperature  of  the  shroud;  »<=-£“-  is  tne  rate  of 

0  -‘U.hCa 

the  regular  regime  (S*  is  the  area  of  the  surface  and  is  the 

o.n  o.n 

mass  of  the  shroud) . 


According  to  boundary  condition  (27),  the  shroud  should,  in 
time,  be  heated  to  the  temperature  of  the  gas  (as  t  -*-<») .  This 
does  not  really  happen,  since  as  a  result  of  the  rise  in  temperature 
of  the  shroud,  equilibrium  gradually  sets  in  between  its  convection 
heating  and  radiation  cooling.  Here  we  cannot  use  boundary,  con¬ 
dition  (18)  for  the  steady  process  because  t^  will  depend,  under 
transient  conditions,  on  time.  For  precise  determination  of  the 
boundary  condition,  considering  the  radiation  heat  losses,  it 
would  be  necessary  to  compile  a  transient  differential  heat-balance 
equation  using  the  Stefan-Boltzmann  law;  because  of  the  difficulties 
involved  in  integrating  this  equation  we  can  use  the  following 
method.  Let  us  introduce  into  (27)  a  correction  factor  which  can 
depend  regularly  on  time;  as  t  -*■  ®  it  becomes  numerically  equal  to 

t *(l~)  -  t  ,  l.e.,  when  x~  =  we  should  have 
r  d  o.n  *  2  d 


/,  t)--- i* (/.)•—  fMLWol  ,.}  ( 1 


(28) 


For  blades  without  a  shroud  we  can  formulate  condition  (27) 
under  certain  assumptions,  which  nevertheless  lead  to  more  precise 
results  than  does  the  condition  3t^( 72,t)/3x2  =  0  for  the  blade  tip. 
In  this  case  the  blade  element  whose  temperature  changes  by  the  law 
of  a  regular  regime  will  be  a  certain  section  near  its  tip.  The 
length  of  this  section  should  be  about  10?  of  the  length  of  the 
working  part  of  the  blade.  Then  coordinate  will  pertain  to  the 
section  of  the  blade  at  the  boundary  between  the  entire  working 
part  and  the  isolated  segment.  Here  the  assumption  that  X  for 
the  blade  tip  section  somewhat  distorts  the  final  temperature 
function.  This  distortion  can  be  corrected  by  arbitrarily  changing 
the  tren..  of  the  temperature  function  t  (x)  when  graphically  analyz¬ 
ing  the  calculation  results.  The  corrected  graph  will  follow  along 
the  upper  boundary  of  the  field  of  scatter  of  the  results  of  calcu¬ 
lating  the  blade  temperature,  and  thus  correspond  to  more  rigid 
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conditions  for  operation  of  the  upper  part  of  the  blade  foil. 

The  examined  method  also  makes  it  possible  to  use  boundary  con¬ 
dition  (28)  for  a  blade  without  a  tip  shroud. 

This  analysis  of  the  boundary  conditions  for  a  transient  pro¬ 
cess  in  the  disk/blade  system  makes  possible  a  well-founded  selection 
of  that  combination  of  boundary  conditions  which  is  close  to  reality. 
Undoubtedly,  particular  cases  can  occur  which  require  a  certain  re¬ 
finement  of  the  examined  boundary  conditions.  These  refinements 
obviously  should  be  introduced  into  the  initial  heat-balance  equa¬ 
tions.  Then,  if  refinements  become  necessary,  the  final  form  of  the 
expressions  for  the  boundary  conditions  when  r  =  r2,  r  =  r^,  x^  =  0, 
and  x2  =  will  hardly  change  in  any  significant  manner. 

The  statement  and  solution  of  two  independent  heat  problems  for 
a  disk  with  blade  roots  (x^)  and  for  blade  tips  (x2)  result  in  the 
fact  that  function  t  (x,r)  for  the  entire  blade  can  have  a dis- 
continuity  at  point  x1  *  l ^  and  x2  *  0.  Practice  has  shown  that 
the  value  of  this  discontinuity  is  very  insignificant,  and  it  is 
much  less  than  the  error  of  the  calculation  method  as  a  whole. 
Therefore  we  can  eliminate  the  discontinuity  of  the  function  t  (x,t) 

J1 

on  the  blade-section  boundary  Xj  and  x2,  by  calculating  the  arith¬ 
metic  mean  temperature  / L-'-'J—  til ^°~  T^».  and  using  it  as  the  only 

value  for  the  coordinates  x^  =  l ^  and  x2  =  0.  The  entire 
function  t^Cxjj)  at  this  point  should  be  corrected  graphically  and 
represented  as  a  continuous  function. 

Thus,  the  known  forms  of  the  boundary  conditions,  encountered 
when  solving  the  problems  of  heat  conduction  for  gas-turbine  rotors, 
do  not  always  correspond  to  the  actual  forms.  The  data  given  in 
the  article  give  grounds  for  considering  that  in  certain  cases  we 
snould  solve  the  problem  of  heat  conduction  for  the  disk  and  blades 
simultaneously.  This  is  advisable  when  the  only  initial  data  for 
calculating  the  temperature  state  of  the  rotor  are  data  from  pre¬ 
liminary  calculation  of  turbine  blading  and  hydrodynamic  calculation 
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of  rotor  cooling. 

The  heat-balance  equations  compiled  for  the  blade  root  and  the 
periphery  of  the  disk  allowed  us  to  formulate,  and  represent  in 
general  form,  the  boundary  conditions  with  x  *  0  and  r  *  ^  both  for 
steady  and  transient  processes  of  heat  conduction.  For  shrouded 
blades  the  boundary  condition  at  the  blade  tip  can  be  compiled  by 
starting  from  the  thermal  state  of  the  upper  shroud  which,  when  the 
turbine  operates  under  steady-state  conditions,  has  a  certain  steady 
temperature  which  is  found  by  solving  the  corresponding  heat-balance 
equation. 

REFERENCES 

1.  )K it p ii uk u it  I'.  C.,  .1  o k a h  B.  ii..  M»kcj  rotJ  M.  K-.  Cipyii- 
K II h  B.  A.,  ra30DNe  T\p6ii!iM  aBiiamiomiux  amtraTe.ieii,  OGopjiinu,  1953. 

2.  IUbcu  H.  T.,  ii  k  G  a  h  E.  11.  Boiujwh'.'c  ox.ia/tuoiimi  poropoa  raicsux 
TypGllll.  iija.  KiieBCKoro  \  miBcpcinera,  1059. 

3  ricTpoBCKiiii  B.  C..  cAnaaunoHMan  •..\xmsHa*.  1069.  ,\Y  1. 


18 


BOUNDARY  CONDITIONS  IN  HEAT-CONDUCTION  PROBLEMS  FOR  COOLED  NOZZLE 
BLADES 

V.  S.  Petrovskiy 

Theoretical  study  of  the  temperature  state  of  nozzle  blades  of 
a  gas  turbine  involves  solution  of  the  boundary- value  problem,  which 
includes  the  boundary  conditions  at  the  blade  tips.  The  accuracy 
with  which  the  boundary  conditions  are  given  determines,  to  a  con¬ 
siderable  extent,  the  accuracy  of  the  final  result. 

Let  us  ex  mine  one  of  the  versions  for  formulating  the  boundary 
conditions  at  the  tips  of  a  hollow  cooled  nozzle  blade  when  solving 
a  one-dimensional  steady  heat  problem. 

At  the  blade  tips  there  are  shrouds,  and  the  tips  are  in  con¬ 
tact  with  the  inner  and  outer  shrouding  strips.  With  steady-state 
operation  of  the  gas  turbine  the  temperature  of  the  shrouds  depends 
little  on  the  design  of  the  cooling  system  within  the  blades.  If 
we  disregard  the  heat  which  is  fed  to  the  shrouds  from  the  working 
part  of  the  blades,  the  shroud  temperature  can  be  found  from  the  con¬ 
dition  of  thermal  equilibrium  between  the  heat  flows  from  the  gas  to 
the  inner  surface  of  the  shroud  and  from  its  inner  surface  to  the 
cooling  air.  Such  thermal  equilibrium  can  be  described  by  the  cor¬ 
responding  equation  for  the  balance  of  the  heat  flows. 
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Pig.  1.  T^t) feuding  of  nozzle 
blades  into  the  shrouds.  The 
arrows  show  the  path  of  the 
cooling  air  within  the  blades : 
a  —  emergence  into  the  tur¬ 
bine  blading;  b  —  emergence 
to  the  base  of  the  rotor 
blades . 


Let  us  examine  a  nozzle  blade 
to  which  cooling  air  is  fed  from  an 
annular  cavity  in  the  upper  shroud 
(Pig.  1).  After  cooling  the  blade, 
the  air  enters  the  blading  of  the 
turbine  through  openings  near  the 
rear  edge  on  the  concave  part  of  the 
profile  (Pig.  la)  or  combines  with 
the  air  flow  which  cools  the  turbine 
rotor,  through  an  opening  in  the 
lower  shroud  (Fig.  lb).  In  both 
cases  the  temperature  of  the  shroud 
is  determined  by  convection  and  ra¬ 
diation  heating  from  the  gas  flow 
and  by  cooling  by  the  air  flow.  To 
compile  the  heat-balance  equation 
let  us  schematize  the  design  of  the 
shroud,  representing  it  in  the  form 


Fig.  2.  Schema¬ 
tized  nozzle  blade. 


of  a  plane  wall  (Pig.  2). 

Let  us  also  give  the  coordinate  system 
in  which  we  will  determine  the  boundary  values 
of  blade  temperature  and  gas- flew  stagnation 
temperature.  The  coordinate  origin  (x  *  0) 
is  made  to  coincide  with  the  inner  surface  of 
the  lower  shroud  adjacent  to  the  inner  shroud 
strip. 


In  the  heat-balance  equations  we  consider  radiation  heat  trans¬ 
fer,  proportional  to  the  fourth  power  of  the  absolute  temperature. 
This  makes  it  necessary  to  express  the  temperature  differential  for 
convection  heat  flows  in  the  form  of  the  difference  between  the  ab¬ 
solute  temperatures  of  the  m°dium  and  the  wall. 


Ultimately  the  equation  for  heat-flow  balance  for  the  upper 
shroud  can  be  written  as  follows: 


q\—ij2=qi+qA, 


where  qt—a,  [7'r(/)-7'I]  is  the  convection  heat  flow  from  the  gas  to  the 

shroud  (a  is  the  coefficient  of  heat  transfer  of  the  gas);  qj—-—’ 
r  ii 

•  l7”j — 7\)  is  the  conduction  heat  flow  through  the  shroud  (X  is  the 
coefficient  of  heat  conduction  of  the  blade,  6^  is  the  shroud 
thickness);  q3-u,{7'2—  T„)  is  the  heat  flow  from  the  outer  surface  of 


the  shroud  to  the  cooling  air  (a  is  the  coefficient  of  heat  trans- 

ri 

fer  from  the  shroud  to  the  air);  q4^=t- — —  is  the  heat  lost  to  ra- 

*  10CM 

diation  from  the  cooled  surface  of  the  shroud  (e  is  the  coefficient 
of  blackness  of  the  shroud  surface,  a  is  the  thermal  radiation  con¬ 
stant)  . 


Equation  (1)  does  not  take  into  account  radiation  heat  trans¬ 
fer  between  the  inner  surface  of  the  shroud,  the  heated  surfaces  of 
adjacent  parts,  and  the  gas  flow.  We  can  assume  that  the  value  of 
this  flow,  determined  by  the  Stefan-Boltzmann  law,  will  be  very 
small  compared  with  the  intensity  of  convection  heat  transfer,  since 
the  difference  in  temperatures  of  the  inner  surface  of  the  shroud 
and  the  surfaces  of  the  surrounding  parts  is  not  great,  while  the 
gas  flow  has  a  low  radiation  coefficient.  In  radiant  heat  flow 
no  account  is  taken  of  the  component  of  radiation  of  the  heated  sur¬ 
face  of  the  shroud,  since  the  temperature  of  this  surface  is  con¬ 
siderably  lower  than  temoerature  T2  of  the  outer  surfs  ■  of  the 
shroud. 


Substituting  the  values  of  q1_i<  into  expression  (1)  and  solving 
the  resulting  equations  for  T^,  we  get 


X  I  (KM 


(2) 


where  Bir— and  —  Z.  are  the  Biot  criteria  for  the  inner  and 
x  x 

outer  surfaces  of  the  upper  shroud. 
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To  determine  T2  for  the  corresponding  values  of  all  quantities 
entering  into  equation  (2),  we  can  solve  the  latter  by  selecting 
that  value  of  T2  for  which  the  equation  becomes  identical.  This 
method  allows  us  to  find  T2  from  the  3rd  or  4th  approximation. 

Using  the  found  value  of  T2,  from  equation  (1)  in  the  form 

[T'r i/)  —  7-,]  =-  A-  I/*,  —  (3) 

we  can  calculate  the  value  of  T^. 

Temperature  T.^  will  thus  be  a  certain  steady  temperature  of 
the  blade  cross  section  when  x  =  l,  which  makes  it  possible  to  wr"- 
the  boundary  conditions 

/.•>(/)=/..  (4) 


If  the  cooling  air  emerges  into  the  cavity  of  the  lower  shroud 
the  steady  temperature  of  the  inner  surface  of  the  lower  shroud  ii ^ 
determined  analogously.  Using  as  our  basis  the  same  concepts  as 
when  compiling  equation  (2),  we  get 


Gi>r>)  .  -  t* 

-■V  ,  Ta  =  — r„)-f 51  —  — 

Rir — 1  4  ).  J00« 


Solving  equation  (5)  for  T^  by  selecting  and  substituting  the 
obtained  numerical  value  of  T^  in  equation  (1)  in  the  form 


r.r{7'U0)~r,]^±  cr. -T,), 

we  get  the  formula  for  calculating  T^. 


(6) 


If  the  cooling  air  emerges  into  the  turbine  blading,  the  heat- 
balance  equation  has  the  following  form: 

(7) 

where  in—  u,  1/7  —  r3j  is  the  convective  heat  flow  from  the  gas  to  the 

inner  surface  of  the  lower  shroud;  <7?=  — « 7~, —  r«>  is  the  heat  flow 
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across  the  shroud;  g\  —  t s  — i-  is  the  heat  lost  to  radiation  by  the 

«00* 

outer  surface  of  the  shroud. 


< 

f, 

| 

J 
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From  equation  (7),  by  substitution  of  the  values  of  qj ,  qi,  and 

Jl  C. 

q/,,  we  get 


_IL. 

(8) 

Bif -rl  *  A 

too* 

B«;(ra0)-7'J=r,- 

T<, 

(9) 

from  which  we  can  calculate  and  then  T^. 


Thus,  in  both  cases  for  the  lower  end  of  a  nozzle  blade  [x  =  0) 
the  boundary  condition  can  be  written  in  the  form 


MO  )=fe  do)  1 

Boundary  conditions  (4)  and  (10),  obtained  as  a  result  of  these 
discussions,  are  definitely  not  the  only  possible  ones.  They  can 
change  depending  on  the  design  of  the  fitting  of  the  nozzle  blades 
into  the  shrouds.  However,  for  the  design  shown  in  Fig.  1  these 
boundary  conditions  give  quite  good  coincidence  between  the  theoreti¬ 
cal  distribution  of  t^(x)  and  the  experimental  distribution:  the 
results  of  calculating  t  (x)  from  the  equations  obtained  when  solving 
the  corresponding  thermal  problem  differed  by  no  more  than  4-5/5  from 
data  on  direct  measurements  at  the  nozzle  blade  tips. 


The  boundary  conditions  obtained  can  be  used  when  setting  up 
steady-state  problems  of  heat  conduction  relative  to  cooled  nozzle 
blades . 
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THE  QUESTION  OF  THE  AVERAGE  HEAT-TRANSFER  COEFFICIENT  IN  VANE 
CASCADES 

N.  A.  Petrovskaya  and  V.  S.  Petrovskiy 

/ 

/ 

/ 

The  heat-transfer  coefficient  a  in  a  vane  cascade  is  one  of  the 
basic  values  describing  the  distribution  of  heat7 in  gas-turbine 
rotors  and  stators.  The  initial  differential  equations  for  transient 
and  steady  heat  conduction,  compiled  for  solving  the  appropriate 
problems,  include  the  coefficient  of  heat  transfer  to  the  blade 
surfaces . 

/ 

/ 

The  intensity  of  the  heat  transfer  oh  blade  surfaces  determines, 
to  a  considerable  extent,  the  temperatui;4  field  in  turbine  rotor 
disks,  since  heating  of  the  disk  is  caused  basically  by  the  flow  of 
heat  from  the  blades.  The  influence  of  the  heat-transfer  coefficient 
in  a  rotor  vane  channel  on  the  temperature  field  of  the  blades  and 
disk  of  the  rotor  is  particularly  noticeable  when  simultaneously 
solving  the  initial  equations  of  heat  conduction  for  the  blades  and 
disk.  In  this  case,  in  final  analjrsis  we  get  temperature  depen¬ 
dences  (linear,  two-dimensional,  or  three-dimensional)  for  the 
blades  and  disk  which  are  connected,  at  the  point  where  they  come 
together,  by  appropriate  boundary  conditions.  These  conditions  in¬ 
clude  the  functions  tn(x)  or  dtJ1(x)/ax  when  x  =  0  (see,  e.g.,  [5]) 
which  also  include  the  value  of  the  heac-transfer  intensity  coef¬ 
ficient  a  on  the  blade  surface  (x  is  a  coo- iinate  along  the  blade). 
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For  two-  and  three-dimensional  heat-conduction  problems  the  distri¬ 
bution  of  heat-transfer  intensity  along  the  blade  profile  is  of 
interest.  Such  data  have  been  obtained  in  certain  investigations, 
but  it  is  difficult  to  use  them  in  calculations  since  it  is  im¬ 
possible  to  describe  the  dependence  a  =  f(x)  by  a  function  which 
would  coincide  with  the  actual  distribution  of  a.  If  this  can  be 
done,  however,  it  is  extremely  complex  to  solve  the  problem  by  intro¬ 
ducing  such  a  function  into  the  Initial  equations’.  Therefore  it 
becomes  necessary  to  replace  it  by  simpler  dependences  which  only 
approximately  correspond  to  the  true  distribution  of  a  along  the 
blade  profile. 

The  distribution  of  heat-transfer  intensity  along  the  blade 
foil  should  have  a  substantial  influence  on  the  distribution  of  tem¬ 
perature  along  the  blade. 

As  shown  in  the  work  by  V.  M.  Kapinos  and  A.  G.  Knabe  [1],  the 
intensity  of  heat  transfer  in  a  vane  cascade  depends  on  the  nature 
of  the  gas  flow  along  the  blade  profile,  which  in  turn  is  determined 
by  the  reaction  of  the  cascade.  Consequently,  the  dependence  of  a 
on  x  can  be  determined  if  we  know  how  the  reaction  of  the  cascade 
changes  along  the  blade. 

The  authors  of  L 1 U  qualitatively  substantiate  the  dependence 
of  heat  transfer  on  the  shaping  of  the  cascade  —  its  reaction.  Act¬ 
ually,  as  shown  in  the  works  by  Ainley  [2,  3]  and  Andrews  and  Bradley 
[4],  the  intensity  of  convection  heat  transfer  is  determined  by  the 
nature  of  the  gas  flow  in  the  blade  channel,  more  precisely,  by  the 
position  of  the  point  of  transition  of  laminar  into  turbulent  flow 
on  the  convex  and  concave  sides  of  the  blade. 

On  active  profiles,  on  profiles  with  small  leading-edge  radius, 
the  length  of  the  segment  with  turbulent  flow  increases,  since  the 
point  of  transition  of  laminar  into  turbulent  flow  is  displaced 
toward  the  leading  edge  of  the  profile.  The  probability  is  especially 
high  of  turbulence  of  the  gas  flow  on  the  convex  side  of  the  profile. 
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where  stability  of  the  laminar  flow  is  easily  lost  due  to  centrifugal 
forces.  The  experiments  of  K.  Bammert  [6],  L.  P.  Lozitskiy  [7];., and 
M.  N.  Bodunov  [8]  on  determining  the  dependence  of  heat  transfer  on 
the  angle  of  attack  (leak  angle)  are  very  indicative  in  this  sense. 
They  show  that  there  is  a  leak  angle  for  which  heat  transfer  will  be 
mini-num.  Thi3  angle,  measured  between  the  flow  direction  and  the 
front  of  the  cascade,  is  somewhat  greater  than  the  calculated  angle 
of  inlet  of  the  gas  flow  into  tne  cascade,  since  when  the  flow  im¬ 
pinges  against  the  convex  side  of  the  blade,  movement  along  the  con¬ 
vex  surface  of  the  profile,  up  to  certain  limits,  is  easier  to 
establish  than  for  the  calculated  inlet  angle.  The  authors  of  [1] 
generalize  the  studies  on  convection  heat  transfer  in  various  cas¬ 
cades  performed  by  various  authors.  They  explain  the  nonagreement 
of  their  obtained  experimental  results  by  the  incompleteness  of  the 
similarity  criteria,  and  to  expand  the  range  of  application  of  the 
generalizing  expression  they  introduce  the  value  y  =  sin  02/sin  0i, 
where  0i  and  02  are  the  calculated  inlet  and  outlet  angles  of  the 
gas  flow  to  and  from  the  cascade.  Angle  0i  is  measured  between  the 
frontal  line  of  the  cascade  and  the  tangent  to  the  convex  part  of 
the  profile  at  the  inlet. 

We  cannot  agree  with  the  authors  of  [1]  that  the  discrepancy 
of  the  data  obtained  in  various  studies  does  not  depend  on  experiment 
errors.  Undoubtedly,  the  errors  arising  when  processing  the  experi¬ 
mental  data  and  which  depend  on  selection  of  the  determining  param¬ 
eters  do  not  play  a  major  role  in  the  spread  of  coefficient  C  and 
exponent  n  in  the  final  criterial  formulas.  Analysis  has  shown  that 
the  ready  formula  in  [1]  is  also  valid  only  within  specific  limits, 
although  it  was  obtained  on  the  basis  of  analysis  of  criterial  form¬ 
ulas  reduced  to  the  identical  determining  magnitudes.  In  our  opinion, 
the  main  reason  for  the  discrepancy  in  the  experimental  data  is 
peculiarities  in  the  experiments.  These  include  differences  in  the 
shapes  of  the  blade  profiles,  relative  pitch,  and  other  structural 
data  of  the  cascade,  and  also  differences  in  various  details  of  the 
experiment  methods:  selection  of  the  direction  of  the  heat  flow, 
ways  of  treasuring  the  calculated  temperatures,  the  amount  of  heat. 
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methods  for  scavenging  the  cascade,  features  in  the  design  of  the 
test  apparatus,  the  forms  and  characteristics  of  the  heating  agents, 
etc.  For  these  reasons  it  is  quite  unlikely  that  tests  of  specific 
cascades  by  various  authors  using,  e.g.,  the  method  presented  in  [12], 
would  yield  identical  final  calculation,  formulas.  We  should  note 
that  the  use  of  any  empirical  formula  can  give  a  precise  result  only 
under  very  specific  particular  conditions,  although  it  might  have 
the  form  of  a  criterial  dependence. 

The  idea  of  generalizing  criterial  formulas  for  an  entire  range 
of  change  of  parameter  y  is  of  great  practical  significance,  but  it 
is  also  advisable  to  find  criterial  relationships  for  narrower  ranges 
of  change  of  the  cascade  characteristics.  In  this  case  we  can  reduce 
the  error  arising  due  to  nonagreement  of  the  generalizing  expression 
with  the  individual  particular  formulas.  It  makes  sense  also  to  find 
generalizing  criterial  expressions  with- constant  values  of  coefficient 
C  and  exponent  n  in  a  formula  of  the  type  Nu  =  C  Re11  for  cascades  of 
the  impulse  or  reaction  type.  Such  formulas  are  suitable  for  pre¬ 
liminary  calculations  of  heat  transfer. 

In  this  article  we  have  performed  repeated  analysis  of  existing 
criterial  formulas  that  describe  convection  heat  transfer  in  gas- 
turbine  vane  cascades.  Here  we  considered  certain  of  the  inaccuracies 
admitted  in  [1].  As  an  additional  similarity  criterion  for  heat 
transfer  in  cascades,  one  which  takes  into  account  the  geometric  and 
hydrodynamic  disparity  of  comparable  processes,  we  used,  unlike  the 
authors  of  [1],  a  convergence  coefficient  k  -  sin  Si/sin  62,  where 
B>  and  are  the  calculated  inlet  and  outlet  angles  of  the  gas  flow 
betvet-n  the  front  line  of  the  cascade  and  the  tangent  to  the  middle 
line  of  t profile  at  the  leading  and  trailing  edges  of  the  blade. 
This  coefficient  is  more  often  used  in  gas-dynamic  and  thermal  cal¬ 
culations  and,  in  our  opinion,  is  more  expedient  as  an  auxiliary 
criterion.  The  coefficients  k  and  y  and  approximately  reciprocals. 

All  initial  criterial  formulas  are  given  in  Table  1  and  presented  in 
order  of  increasing  convergence  coefficients  of  the  cascades  1 
which  they  were  obtained.  Likewise,  the  reaction  of  the  cascades 
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A,  p,  p  —  heat  conduction,  density,  and  dynamic  viscosity  of 
the  gas  ; 

t  ■  0.5(ti  +  t2);  p  »  0.5(pi  +  P2); 

w  =  0.5(wi  +  w2);  t  *  O.SCt^  +  t*); 

tl’  t2  ~  statlc  temperatures  of  the  gas  flow  at  the  cascade 
inlet  and  outlet; 

t*  -  temperature  of  stagnant  gas  flow; 
tj  —  the  same,  measured  at  cascade  inlet; 

P1J  p2  ~  static  pressures  at  cascade  inlet  and  outlet; 

wl5  w?  ~  relative  speeds  of  gas  flow  at  cascade  inlet  and 
outlet; 

t^  —  temperature  of  blade  surface; 

t’  -  cascade  pitch;  b*  -  cascade  width. 


will  increase  from  the  beginning  to  the  end  of  the  table.  For  form¬ 
ulas  6,  12,  16,  17,  and  18,  which  were  obtained  from  study  of 
heat  transfer  as  a  function  of  leak  angle,  the  values  of  C  and  n  were 
taken  for  zero  angle  of  attack.  When  processing  the  measurement  re¬ 
sults,  as  the  determining  dimension  in  various  cases  we  used  the 
profile  chord  b  or  the  equivalent  diameter  d  .  The  determining  tem- 
perature  also  varied.  In  order  to  be  able  to  compare  the  results  of 
the  investigations,  all  criterial  formulas  were  reduced  to  a  single 
determining  dimension  and  a  single  determining  temperature.  As  the 
determining  dimension  we  used  the  equivalent  diameter  d  =  P/tt, 
which  is  the  one  most  often  encountered  in  studies  on  heat  transfer 
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In  vane  cascades;  as  the  determining  temperature  t  we  used  the  arith¬ 
metic  mean  of  the  gas  temperatures  at  the  cascade  inlet  and  outlet. 
This  determining  temperature  is  obviously  more  advisable  that  deter¬ 
mining  temperatures  which  include  the  blade  temperature;  usually  the 
latter  is  an  unknown  value  and  therefore  it  must  be  specified  for 
calculating  heat-transfer  coefficients  in  a  theoretical  study  of  the 
temperature  states  of  a  turbine  stator  or  rotor.  The  use,  in  the 
calculations ,  of  an  approximate  value  for  the  determining  tempera¬ 
ture  leads  to  an  error  in  the  value  of  the  heat-transfer  coefficient 
and  thus  reduces  the  accuracy  of  the  desired  temperature.  The  gas 
density  led  to  the  selected  determining  temperature  and  pressure  p 
which  is  the  arithmetic  mean  of  the  static  pressures  of  the  gas  at 
the  cascade  inlet  and  outlet.  As  the  determining  velocity  for  the 
Reynolds  number,  in  acco'.-dance  with  recommendations  in  [1],  we  used 
the  arithmetic  mean  of  the  velocities  of  relative  motion  of  the  gas 
in  the  inlet  and  outlet  sections  of  the  vane  cascade  (w) . 


The  initial  formulas  are  reduced  to  these  determining  parameters 
with  consideration  of  the  concepts  presented  in  [1],  As  the  initial 
transfer  equation  we  used 


Basically,  a  change  in  C  involves  transfer  to  a  new  determining 
dimension  and  new  determining  speed.  Therefore,  dropping  the  ratio 
of  determining  temperatures,  equation  (1)  can  be  written  in  the  form 


The  ratio  p  /p  in  all  cases  was  equal  to,  or  close  to,  one.  Hence 
hcx  np 


The  flow  equation  for  the  inlet  and  outlet  sections  of  the  cascade 
allows  us  to  represent  expression  (3)  in  the  following  form: 
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Here  k  is  the  convergence  coefficient,  and  Vg  s.re  the  specific 
volumes  of  the  gas,  and  and  Zg  are  ^he  heights  at  the  cas¬ 

cade  inlet  and  outlet. 


To  obtain  criterial  expressions  which  generalize  all  formulas , 
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ces  n(k).  The  points 
show  the  values  of  n 
from  individual  experi¬ 
ments  . 


t,  p,  and  w  from  Table  1,  we  first  found 

the  dependence 


n(k)  (Fig.  1). 

Then,  from  these 
values  of  n  for 
each  formula  in 
Table  1  which  gave 
the  dependence 
n(k),  and  from  the 
corresponding  aver¬ 
age  values  of  cri¬ 
teria  Nu  and  Re, 
we  found  the  thus 
corrected  values  of 
constants  C.  The 


obtained  constants  C  were  used  to  find  the  dependence  C(k)  (Fig.  2). 
As  a  result,  for  convection  heat  transfer  in  any  type  of  cascade  we 
obtained  the  following  criterial  formula: 


Nu-CWRc^1. 

where  C  =  0,328k  -  0.282,  and  n  =  0 . 736k-0 .  As  the  calculations 
showed,  the  difference  in  the  values  of  the  heat-transfer  coefficients 
calculated  from  formula  (5)  and  from  the  appropriate  formulas  in 

it  5 

Table  1  depends  on  Re.  For  the  interval  Re  a  7*10  to  3*10  it  does 
not  exceed  3-8?,  and  only  in  four  cases  does  it  reach  15? •  For  Re  = 

=  10**  to  7*10^  and  3*10^  to  10^  the  discrepancy  increases  to  10-15?, 
and  to  20-30?  in  individual  cases.  The  obtained  formula  can  be  used 
to  calculate  the  average  coefficient  of  heat  transfer  along  the  blade 
profile  and  along  the  vane  cascade.  For  this,  if  k  depends  on  the 
cascade  height,  it  is  necessary  to  find  its  average  value  a3orig  the 
entire  blade  and  then,  from  this  value  of  k,  determine  n  and  C  in 


calculation  formula  (5)  using  graphs  (Figs.  1  and  2).  As  the  average 
value  we  can  use  the  half-sum  of  the  convergence  coefficient  at  the 
blade  root  and  tip,  since  the  dependence  of  k  on  the  coordinate  along 
the  blade  is  close  to  linear. 

When  compiling  the  boundary  condition  for  the  blade  surface  with 
hermetic  twisting,  i.e.,  with  variable  reaction  along  the  blade,  we 
must  take  into  consideration  the  change  in  heat  transfer  along  the 
blade  foil.  For  this  we  should  know  the  dependence  k(x),  which  we 
then  substitute  into  equation  (5),  from  which  we  can  obtain  the  de¬ 
pendence  «(x).  To  simplify  the  equation  for  the  boundary  conditions 
we  can  recommend  the  following  method.  For  specific  values  of  k 
corresponding  to  the  coordinates  x  =  0  (blade  root),  x  *  l  (blade 
tip),  and  certain  intermediate  values  of  x,  we  calculate  the  heal 
transfer  coefficients  a  from  formula  (5).  From  the  found  values  of 
a  we  can  construct  the  dependence  ct(x)  and  describe  it  by  the  most 
suitable  and  simple  equation. 

Equation  (5),  written  for  the  entire  range  of  change  of  k,  leads 
to  considerable  spread  of  the  results  ax  the  end  of  this  interval, 
i.e.,  in  those  regions  where  cascades  of  a  specific  type  are  grouped. 

To  refine  equation  (5),  if  there  is  no  need  for  using  it  over  a  broad 

range  of  change  of  k,  we  can  propose  linear 
dependences  of  the  change  in  C  and  n  for 
impulse  and  reaction  blades,  respectively. 

As  such  dependences  (Figs.  1  and  3)  we  will 
have,  for  impulse  cascades  (0.9  <  k  <  1.5) 

C=0,IG5/.’— 0,104;  0,835-0,125*  (6) 

and  for  reaction  cascades  (1.5  <  k  <  2.8) 

C= 0,1 77« -i  0,012;  0,032-0,04 7ok.  (7) 

Both  the  initial  as  well  as  the  general¬ 
izing  criterlal  formulas  allow  us  to  cal¬ 
culate  the  heat-transfer  coefficient  only  with  an  error  whose  magnitude 


Fig.  3.  Dependence  C(n) 
for  0.9  <  k  <  1.5  (im¬ 
pulse  cascade)  and  for 
1.5  <  k  <  2.8  (reaction 
cascade ) . 
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depends  not  only  on  the  accuracy  of  the  calculation  formula.  It  is 
determined,  in  addition,  by  the  degree  of  agreement  of  the  gas- 
dynamic  parameters  used  as  the  initial  ones  for  calculating  a  from 
the  appropriate  formulas,  their  true  value,  the  incomplete  similarity 
of  the  actual  process  of  heat  transfer  in  the  cascade  to  that  de¬ 
scribed  by  the  calculation  formula,  and  a  number  of  other  factors 
mentioned  above.  In  this  connection  it  is  advisable  to  examine  how 
the  inaccuracy  of  the  value  of  a  for  a  vane  cascade  influences  the 
error  of  the  calculated  value  of  the  blade  or  disk  temperature. 

As  calculations  have  shown,  the  decisive  factors  for  the  accu¬ 
racy  with  which  the  temperature  field  in  blades  or  a  disk  can  be 
determined  are  the  correctness  of  the  formula  for  the  boundary  con¬ 
ditions  and  consideration,  in  the  initial  differential  equations, 
of  the  change  in  gas  temperature  along  the  vane  cascade.  In  addition, 
the  coefficient  a  often  enters  to  the  0.5  power  in  the  calculation 
equations  used  for  theoretical  determination  of  the  temperature  of 
gas-turbine  rotors  or  disk.  This  reduces  the  value  of  the  relative 
error,  introduced  by  an  inaccurate  value  of  the  heat-transfer  coef¬ 
ficient  into  the  final  result,  by  one-half. 

From  this  it  follows  that,  in  addition  to  generalizing  formula 
(5),  it  makes  sense  to  find  approximate  criterial  relationships  which 
would  allow  us,  with  acceptable  accuracy,  to  calculate  relatively 
simply  the  average  heat-transfer  coefficient  on  the  blade  surface 
for  cascades  with  impulse  or  reaction  profiles.  As  can  be  seen  from 
Fig.  1,  the  finding  of  such  relationships  is  simplified  by  the  fact 
that  in  the  entire  field  of  scatter  of  points  n(k)  we  observe  an 
expressed  tendency  toward  grouping  around  two  average  values:  with 
k  <  1.5  and  with  k  >  1.5.  This  can  be  explained  by  the  fact  that 
the  nature  of  the  gas  flow  In  the  vane  channel  of  a  cascade  with  a 
value  of  k  of  about  1.5  is  such  that  the  position  of  the  point  of 
transition  of  laminar  into  turbulent  flow  is  unstable.  Therefore, 
the  relationshop  between  the  extent  of  the  sections  of  laminar  and 
turbulent  flow  continually  changes  and,  consequently,  the  value  of 
the  exponent  n  for  Re  will  also  be  unstable. 
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As  the  criterial  formulas  with 
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constant  values  of  C  and  n  we  can 
propose  the  following  (Figs .  1  and 
4):  for  impulse  cascades  (k  <  1.5) 

Nu-  0,085  ReW»;  (8) 

for  reaction  cascades  (k  >  1.5) 

N'u= 0,326  Rec'-1 2 * 4 5 6 7 8 9 10 11 12 13*.  '  (9) 


In  conclusion  we  should  note  that  in  this  article  we  have 
generalized  the  results  of  studies  of  convection  heat  transfer  in 
nonmoving  cascades.  In  a  moving  cascade  the  heat-transfer  coefficient 
coefficient  increases  1.3-1. 4- fold  compared  with  a  nonmoving  one. 
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HEAT-TRANSFER  COEFFICIENTS  ON  THE  SURFACES  OF  THE  BLADED  ROTOR  OF 
AN  AXIAL  GAS  TURBINE 

V.  Z.  Kitayev  and  V.  S.  Petrovskiy 


The  diversity  of  conditions  for  convection  heat  transfer  in 
gas  turbines  and,  in  particular,  on  rotor  surfaces,  and  the  complex¬ 
ity  of  a  clear  determination  of  the  hydrodynamics  of  this  process 
make  the  finding  of  the  heat-transfer  coefficient  values  an  inde¬ 
pendent,  specific,  important  problem.  Existing  recommendations  for 
calculating  heat-transfer  coefficients  are  based  both  on  ie  use  of 
generally-known  criterial  formulas,  obtained  for  certain  comparatively 
simple  cases  of  convection  heat  transfer,  as  well  as  on  the  results 
of  individual  experimental  studies  under  conditions  close  to  those  of 
heat  transfer  in  real  operating  turbine  constructions.  Up  to  the 
present,  convection  heat  transfer  in  vane  cascades  has  been  studied 
most  completely. 

As  a  rule,  the  calculation  formulas  are  obtained  from  study  of 
an  immobile  object  consisting  of  several  blades.  For  an  operating 
gas  turbine,  however,  these  formulas  are  applicable  only  after  care¬ 
ful  analysis  of  the  similarity  of  the  specific  case  of  heat  transfer 
to  that  for  which  the  criterial  equation  employed  was  obtained. 

The  need  for  conducting  the  operations  discussed  below  arose  in 
connection  with  solution  of  the  problem  of  heat  conduction  for  a 
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system  turbine  disk/rotor  blade.  Unlike  problems  solved  individually 
for  a  disk,  blade,' or  other  part  of  a  gas  turbine,  it  was  necessary 
to  introduce,  into  the  boundary  conditions,  the  heat-transfer  coef¬ 
ficient  a  on  various  surfaces  of  a  bladed  rotor.  Naturally,  the 
accuracy  of  the  calculation  results  was  determined,  to  a  considerable 
extent,  by  the  sum  of  the  errors  which  can  be  introduced  into  the 
problem  by  coefficients  a.  This  necessitated  a  more  accurate  valida¬ 
tion  of  selection  of  the  criterlal  formulas,  explanation  of  the  gen¬ 
eral  picture  of  convection  heat  transfer  over  the  entire  gas-turbine 
rotor,  and  determination  of  the  dependence  of  the  heat-transfer 
coefficients  on  the  turbine  operating  conditions  and,  what  is  par¬ 
ticularly  important,  on  the  ambient  conditions  characterized  by 
specific  tH  and  pH. 

Errors  in  the  values  of  the  heat-transfer  coefficients  can 
occur  due  to  discrepancies  in  the  real  conditions  of  heat  transfer 
from  those  conditions  under  which  the  calculation  formula  was  ob¬ 
tained,  and  due  to  errors  in  the  values  of  the  initial  determining 
parameters  (the  temperature  of  the  flow  medium,  its  pressure,  flow 
rate,  and  determining  dimensions). 

With  combined  study  of  convection  heat  transfer  in  the  entire 
rotor  it  was  possible,  by  comparing  various  coefficients  a,  to  ex¬ 
plain  how  correctly  the  values  of  the  initial  determining  parameters 
were  chcsen. 

Figure  3  of  the  first  article  in  this  collection  (see  p.  6) 
shows  a  diagram  of  a  shrouded  rotor  blade  with  base  attached  to  a 
disk,  and  gives  values  for  the  average  heat-transfer  coefficients 
entering  into  the  initial  equations  of  heat  conduction  and  into  the 
boundary  conditions  when  solving  the  appropriate  problem  for  a 
specific  three-stage  axial  turbine. 

Let  us  examine,  in  sequence,  each  heat-transfer  surface  of  in¬ 
terest  to  us,  and  validate  the  selection  of  some  criterlal  expression 
for  calculating  the  average  value  of  a. 
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The  upper  blade  shrouds,  closing,  create  two  heat-transfer  sur¬ 
faces:  an  external  cylindrical  surface  and  an  internal  vane  surface 
on  the  side  toward  the  gas  flow.  When  stating  the  problem  of  heat 
conduction  the  boundary  condition  at  the  blado  tip  is  as  follows 
(when  x  =  l) : 


f.i(A  t)  — /i.n(f) 


for  a  transient  process,  and 


for  a  steady  process, 


t.i  0)  —  ^n.n- 


(1) 


(2) 


The  value  of  t  is  determined  from  the  heat-balance  equation 
a .  n 


for  v^e  upper  shroud 


(3) 


where  a  =  0.5(0^  +  a2),  and  T*  and  Tq  n  are  in  °K. 


The  coefficient  a  also  appears  in  the  expression  for  t  _(x). 
Convection  heat  transfer  on  the  outer  and  inner  surfaces  of  the  upper 
shroud  occurs  under  varying  conditions.  We  can  assume  that  the  de¬ 
termining  temperature  and  pressure  will  be  identical  on  the  outer  and 
inner  surfaces  of  the  shroud.  Heat  transfer  on  the  outer  cylindrical 
surface  moving  relative  to  the  stator  surface  can  be  examined  as  con¬ 
vection  in  a  slot  channel  or  as  heat  transfer  in  a  plane  of  length  ud. 
The  high  degree  of  turbulization  of  the  gas  in  the  gap  between  the 
outer  shroud  surface  and  the  stator  surface  does  not  allow  us  to 
assume  laminar  flow  in  the  gap.  Even  if  we  consider  that  flow  of 
the  gas  in  the  axial  direction  creates  relatively  low  washing  speeds, 
the  high  peripheral  speeds  of  the  rotor  result  in  high  Reynolds  num¬ 
bers  in  the  gap  and,  consequently,  high  gas  turbulence.  Such  con¬ 
ditions  o'^  convection  do  not  allow  us  to  use  the  criterial  formulas 
recommended  for  slot  gaps  and  labyrinth  seals:  this  would  lead  to 
a  lowering  of  the  heat-transfer  coefficient.  In  this  case  it  is 
advisable  to  use  a  formula  that  describes  convection  heat  transfer 
in  a  plane  with  turbulent  motion  of  the  washing  medium: 


* 
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Nu»  0,032  Reow. 


(4) 

Theoretical  gas-dynamics  calculation  shows  that  at  the  upper 
shrouds  the  temperature  of  the  gas  is  somewhat  higher  than  at  the  • 
level  of  the  middle  section  of  the  blade.  This  was  not  actually 
observed.  The  temperature  of  the  gas  in  the  turbine  blading  is 
sharply  reduced  in  the  lower  and  upper  parts  of  the  vane  cascades. 
Consequently,  as  the  determining  temperature  we  must  use  not  the 
gas  temperature  obtained  from  gas-dynamics  calculation  but  a  some¬ 
what  lower  value. 

In  the  absence  of  expe_j mental  data  on  gas-temperature  distri¬ 
bution  along  the  vane  cascade  we  can,  baaed  on  numerous  measurements, 
consider  it  equal  to  0.85-0.90  of  the  gas  temperature  in  the  middle 
section  cf  the  blade.  Then  the  determining  temperature  in  (4)  is 

(5) 

Here  t-jj  and  t 2II  are  the  temperatures  in  the  middle  section 
of  the  blade  at  the  inlet  to  and  outlet  from  the  rotor  cascade  for 
the  appropriate  stage,  taken  from  gas-dynamics  calculation. 

The  determining  pressure  and  the  relative  speed  will  be 

p-0$(Pii+P:n)  and  H'-ndit.  (6) 

(Subscripts  1  and  2  pertain  to  the  cascade  inlet  and  outlet,  re¬ 
spectively;  the  values  of  pj  and  p-j-j  are  taken  from  gas-dynamics 
calculation;  n  is  the  number  of  turbine  rotor  revolutions  per  unit 
time;  d  is  the  diameter  of  the  upper  shroud  of  the  rotor  cascade.) 

The  determining  dimension  in  formula  (4)  should  be  plate  length 
l.  In  this  case  it  is  advisable  to  use  the  length  of  the  generatrix 
of  the  cylindrical  surface  ird  as  the  determining  dimension. 

The  coefficient  of  heat  transfer  on  the  inner  surface  of  the 
upper  shrouds  (in  the  vane  channel)  can  be  determined,  based  on  the 
following  concepts.  The  bottom  of  the  vane  channel  os.n  be  represented 
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as  a  flat  surface  whose  length  is  equal  to  the  length  of  the  cascade 
blading,  as  part  of  a  short  channel  of  variable  section,  and, 
finally,  as  one  of  the  surfaces  of  the  vane  cascade  for  which  the 
value  of  ap  is  determined  from  known  criterial  expressions ,  There 
are  no  equations  for  convection  heat  transfer  for  short  channels  of 
variable  cross  section.  The  use  of  the  cril  ‘•'I  relationships 
recommended  for  a  vane  cascade  is  also  unacceptable  here.  The  flow 
of  gas  along  the  blade  profile,  to  which  the  formulas  for  convection 
heat  transfer  in  vane  cascades,  obtained  in  various  studies,  obtain, 
is  of  a  different  nature  than  flow  along  the  inner  surface  of  the 
upper  shroud.  Correspondingly,  there  is  no  hydrodynamic  similarity 
of  the  processes,  i.e.,  the  criterial  formulas  for  convection  heat 
transfer  should  be  different.  We  see,  in  the  given  case,  just  as 
for  the  outer  shroud  surface,  that  the  most  suitable  formula  for 
calculating  a  will  be  formula  '4)  in  which  the  determining  tempera¬ 
ture  and  pressure  remain  the  same,  while  the  determining  velocity 
is  equal  to  the  arithmetic  mean  of  the  gas  velocities  at  the  cascade 
inlet  and  outlet  in  the  upper  section  of  the  blade  (subscripts  I, 

II,  and  II  in  Pig.  1,  p.  56): 


u?==0,5(»ui  +  ^’2n). 


(7) 


The  determining  dimension  will  be  the  length  of  the  vane  chan¬ 
nel  along  its  middle  line. 


■a 


The  next  heat-transfer  coefficient  is  a^,  which  describes  heat 
transfer  in  the  vane  cascade  itself.  To  determine  this  coefficient 
we  have  many  criterial  relationships  of  the  form 

Nu-CRe",  (8) 

differing  in  the  values  of  C  and  n.  These  relationships  were  ob¬ 
tained  in  experiments  with  various  cascades,  with  various  gas-flow 
leak  angles,  for  a  broad  range  of  change  of  temperatures  and  other 
parameters  describing  heat  transfer.  At  present  attempts  are  being 
made  to  find  a  more  or  less  universal  criterial  relationship  for 
calculating  coefficient  a.  The  first  results  obtained  require  ad¬ 
ditional  verification  and  comparisons.  Therefore,  it  is  most  advis¬ 
able  to  use  a  formula  which  takes  into  account  a  sufficient  number 
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of  similarity  criteria  for  processes  of  heat  transfer  in  vane  cas¬ 
cades.  Such  a  formula 


has  been  obtained  at 
-  T  „/T 


n  r 
ratio. 


is  a  temperature  factor;  *  y/b  is  the  cascade  pitch-chord 


Nu  =  C  Rc0t/jKT0SA'p° 12  (9) 

the  Kazan’  Aviation  Institute  [1].  Here  K, 


T 


The  value  of  C  .n  (9)  depends  on  the  radius  of  the  leading  edge 
of  the  blade.  The  dstermining  temperature  and  pressure  are  found 
as  the  arithmetic  rae&n  for  the  cascade  inlet  and  outlet: 

/ 


i=0,5(<in+/;n)  and  pn^QMpm—pm)- 
dimension  we  use  the  equivalent  diameter  d. 


3K 

ning  velocity  is  found  from  formula  (7)»  except 
ies  are  taken  at  the  level  of  the  middle  section 
fmula  (9)  was  obtained  for  a  static  cascade.  Ac- 
recommendatlons ,  the  coefficient  a  for  the  rotor 
turbine  should  be  increased  by  a  factor  of 


As  the  determining 
a  F^/ir.  The  determJ] 
that  the  gas  velocilj 
(subscript  II).  Fo 
cording  to  existing 
cascade  of  a  working 
1.25-1.30. 

The  heat-transfer  coefficient  on  the  inner  (turned  toward  the 
gas  flow)  surface  of  the  lower  shroud  is  determined  based  on  the 
same  concepts  as  for  the  inner  surface  of  the  upper  shroud.  Here 
the  determining  temperature  will  be  the  same  as  for  a2  (i.e., 

=  tj).  The  determining  pressure  and  velocity  are  calculated  from 
formulas  (6)  and  (7),  but  the  Initial  values  of  p  and  w  are  taken 
for  the  lower  blade  section  (subscript  III). 


Stages  I  and  II  of  a  turbine  have  unique  heat  chokes  in  the  form 
of  cooling  channels  in  the  blade  roots  between  the  lower  shroud  and 
the  disk  crown.  The  Influence  of  the  passage  of  air  through  these 
channels  on  the  total  temperature  field  of  the  blade  and  disk,  as 
appropriate  calculations  have  shown,  is  extremely  great.  In  this 
connection  the  coefficient  a^,  which  characterizes  heat  transfer  in 
the  cooling  channels,  is  of  particular  importance. 


iJO 
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An  error  in  calculating  coefficient  ot^  occurs  mainly  not  due  to 
the  selected  criterial  relationship* 3  not  corresponding  to  the  exam¬ 
ined  heat-transfer  conditions,  but  due  to  the  low  accuracy  in  deter¬ 
mining  the  rate  of  flow  of  cooling  air  through  the  channels .  The 
air  speed  here  can  be  determined  by  using  the  continuity  equation, 
but  the  initial  values  of  the  flows  through  each  channel  are  not  re¬ 
liable.  It  should  be  noted  that  for  any  turbine,  in  the  absence  of 
well-organized  air  flov?,  its  flow  rate  for  cooling  the  individr'J 
rotor  elements  can  only  be  determined  with  high  error. 

At  the  same  time,  gas-dynamics  calculation  for  the  blading  gives, 
with  sufficient  accuracy,  values  of  the  gas  pressure  at  the  vane 
cascade  inlet  and  outlet  at  the  base.  If  we  consider  that  the  pres¬ 
sure  differential  from  the  cavity  ahead  of  the  cooling  chanel  to 
the  turbine  blading  is  approximately  equal  to  the  differential  behind 
the  cooling  channel,  the  pressure  differential  of  the  cooling  air  in 
the  channel  of  a  turbine  of  the  given  design  is  approximately  equal 
to  Ap  =  P-jjjj  -  Pjin*  Using  this  value  for  the  differential  we 
can  obtain  an  expression  for  determining  the  velocity  through  an 
individual  channel  in  the  form  of  the  Bernoulli  equation  for  arbi¬ 
trary  sections  ahead  of  the  inlet  to  the  channel  and  at  its  outlet: 


Here  w,  =  u  sin  9  is  the  speed  of  the  cooling  air,  occurring 
u 

due  to  the  compressor  effect  during  peripheral  movement  of  the  blade 
ro\  positioned  at  angle  9  to  the  flow  axis,  u  is  the  peripheral 
speed  of  the  cooling  channels,  p  and  p  are  air  pressure  and  density. 
Subscripts  1  and  2  pertain  to  sections  at  the  channel  inlet  and  outlet. 


If  we  assume  that  px  -  p2  *  Pmj  “  ^III’  then  Ap  is  known* 
Speed  w2  is  taken  as  the  determining  one  when  calculating  a^.  We 
get  the  calculation  formula  for  w2  from  equation  (10): 


/~  \  t 

= y  —  V-yp-Qi^uj  ■ 


(ID 
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^he  resistance  coefficient  c  must  be  determined  considering  the 
recommendations  presented  in  [2].  It' depends  on  the  speed  of  rotation 
of  the  disk  and  is  found  from  the  formula 


where  ' - - - —  is'  the  flow  factor;  vr‘  is  the  flow 

It  0.5a  . 

factor  for  a  static  disk;  is  the  hydraulic  resistance  of  a  short 

channel  with  different  inlet  and  outlet  cross  sections;  K  «  u/w2  is 
a  parameter  which  defines  the  dependence  of  £  on  the  peripheral 
speed  for  known  w2. 

f' 

c 

,  Prom  the  formulas  derived  we  see  that  before  finding  the  value 
of  w2  we  must  know  the  coefficient  C  which,  in  turn,  depends  on  w2< 
Here  we  can  proceed  as  follows.  First,  from  the  flow  equations  find 
'the  approximate  value  of  w2<  Then  from  formula  (11)  find  the  more 
accurate  value  for  the  speed  in  the  cooling  channel  and  repeat  the 
calculation  using  the  corrected  value  of  C.  In  determining  w2  using 
the  continuity  equation  we  used  the  results  of  a  preliminary  tenta¬ 
tive-hydraulic  calculation.  An  initial 'approximate  value  of  w2  is 
required  hpth  for  parameter  K  as  well  as  for  determining  coefficient 
T ’ •  This  coefficient  depends  on  the  ratio  of  the  areas  of  the  channel 
inlet  and  outlet  (S^/S2),  and  on  tne  Reynolds  number,  which  is  cal¬ 
culated  from  the  initial'  value  of  w2  [3].  The  refined  speed  of  the 
cooling  air  in  the  channels  is  then  used  to  calculate  the  Reynolds 
number;  as  the  determining  dimension  we  used  the  equivalent  channel 
diameter  d  =  4S/P,  where  S  is  the  channel  section  and  P  is  its 

3h 

perimeter. 


As  calculations  have  shown,  the  Reynold's  number  for  cooling 

'  h. 

channels  of  stage  I  of  a  rotor  does  not  exceed  10  ,  while  for  stage 
II  it  does  not  exceed  10  .  Consequently,  the  airflow  here  can  be 
considered  laminar,  or. at  least  transitional,  out  not  turbulent. 
Then  the  formula  for  determining  a^-  will  have  the  form 

A’o=0,!7  Re033  c;'  (13) 
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here  Kq  =  NuPr  is  a  complex  that  depends  on  the  Reynolds  number; 

Sj  is  a  coefficient  which  takes  into  account  the  dependence  of  heat 
transfer  on  the  ratio  1/ d  for  short  channels. 

H 

The  values  of  Kq  for  the  appropriate  Reynolds  number  can  be 
found  from  graphs  given  in  [4].  The  coefficient  a,  is  defined  in  this 
same  work.  The  determining  temperature  and  pressure  when  using  for¬ 
mula  (13)  will  be  the  ambient  temperature  and  pressure. 

The  next  neat-transfer  coefficient  of  interest  to  us  is  ctg, 
which  describes  convection  heat  transfer  on  the  side  surface  of  the 
disk.  To  calculate  ctg  it  is  expedient  to  use  the  formula 

Nu =0,0348  Rc0'77,  (l1*) 

obtained  at  the  Kazan'  Aviation  Institute.  As  parameters  we  use 
here  the  air  temperature  and  pressure  and  the  diameter  of  the  disk 
bed.  The  Reynolds  number  is  calculated  from  the  relative-  air  speed. 
The  air  flows  for  disk  cooling  are  relatively  slow,  and  therefore  the 
relative  speed  found  from  the  continuity  equation  can  be  determining 
only  for  the  central  part  of  the  disk.  For  the  peripheral  region, 
however,  It  is  advisable  to  consider  the  relative  air  speed  equal  to 
the  peripheral  speed  of  the  rotating  disk. 

The  examined  heat-transfer  coefficients  were  calculated  for  two 
turbine  operating  modes  -  idle  and  nominal  (rated)  —  at  ambient  air 
temperature  tH  =  -50,  -30,  >15,  and  +50°C.  The  data  obtained  are 
represented  in  the  form  of  graphic  dependences  in  Figs.  1-5.  Opera¬ 
tion  of  a  turbine  in  the  nominal  regime  is  characterized  by  the  depen¬ 
dences  of  coeificiont  ct^_^  on  t^,  shown  in  Fig.  3.  The  same  order  of 
magnitude  of  these  coefficients  is  observed  for  most  gas  turbines. 
Convection  heat  transfer  is  most  intense  in  the  vane  cascade.  This 
has  the  result  tnat  with  low  heat-transfer  intensity  to  the  disk, 
the  temperature  of  solid  blades  is  almost  the  same  as  that  of  the 
gas  throughout  the  length  of  the  rotor  section,  except  for  a  snail 
area  near  the  root. 


The  decrease  in  the  heat-transfer  coefficients  with  an  increase  1 

in  tjj  is  explained  by  the  associated  decrease  in  air  density  and,  .„J 

correspondingly,  the  gas  density.  The  pressure  of  the  gas,  and  its  i 

temperature,  decrease  from  stage  to  stage.  The  combined  influence  s 

of  these  factors  also  results  in  decreased  gas  density.  Therefore  j 

the  heat-transfer  coefficients  become  less  in  each  successive  stage.  : 

The  change  in  speed  of  the  gas  flow  has  a  lesser  influence  on  the  H 

reliction  in  coefficients  a  from  stage  to  stage.  •  5 

Fi sures  2a  and  2b  show  the  generalized  results  of  calculating 
the  coef  "icients  for  all  rotor  stages  and  both  turbine  opera¬ 
ting  modei  in  the  form  of  the  dependences  ! 


respectively,  where  TQ  «  273°K  (or  0°C).  When  constructing  these 

graphs  the  values  of  a™  were  found  from  Fig. 

i0 

lation.  All  calculated  points  lie  quite  well 

“r/,«‘,r.(!.6->-0.0023r//l 

(Fig.  2a)  and 

«r/#=l.l6ar^2,»7-^-)  ’ 

(Fig.  2b).  Dependence  (16)  is  more  accurate, 
scatter  of  the  calculation  da^a  is  narrower. 

Formulas  (15)  and  (16)  can  be  used  to  determine  the  heat-trans¬ 
fer  coefficients  for  various  TH,  but  we  must  ilrst  know  the  coef¬ 
ficient  a,p  of  each  stage,  which  complicates  the  calculation.  It 

is  considerably  simpler  to  use  the  dependence  of  the  heat-transfer 
coefficients  of  each  successive  turbine  stage  on  the  heat-transfer 
coefficients  of  stage  I  (Fig.  3). 

In  constructing  the  graph,  along  the  horizontal  we  plotted  the 

values  of  cu,  for  stage  I,  and  along  the  vertical  -  the  corresponding 
0 
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1,  as  well  as  by  calcu- 
on  the  line 

(15) 

(16) 

since  the  field  of 


values  of  cu  for  stage  II.  Then  along  the  horizontal  we  plotted 
i0 

cu  for  stage. II,  and  along  the  vertical  -  the  values  for  stage  III. 
x0 

The  connection  between  the  heat-transfer  coefficients  with  Tq  « 

=  273°K  for  various  stages,  based  on  this  graph,  can  be  represented 
by  the  formula 


where  n  is  the  sequence  number  of  the  turbine  stage. 


Substituting  this  expression  into  (16)  in  place  of  a.  ,  we 

i0 

obtain  the  following  dependence,  which  is  more  convenient  for  prac¬ 
tical  use: 

U6ar,i  1.34‘~".  (18) 

The  change  in  coefficient  a,.  in  the  cooling  channels  of  stages 
I  and  III  of  a  turbine  rotor  is  shown  in  Pig.  4.  From  the  absolute 
values  of  a,.  we  see  that  this  coefficient  is  approximately  one  order 
of  magnitude  smaller  than  c^,  a2,  a^,  and  a^.  Because  of  the  ex¬ 
tremely  weak  air  flow  through  the  cooling  channels  of  stage  II,  the 
value  of  a-  here'  is  practically  independent  of  T^.  For  stage  I  this 
coefficient  is  very  sensitive  to  a  change  in  temperature  of  the 
ambient  air,  and  decreases  abruptly  when  the  latter  is  increased. 

As  was  noted  above,  coefficient  has  a  very  slight  influence  on 
the  nature  of  temperature  dependence  t^x)  and  t^(r);  consequently, 
its  sharp  decrease  with  an  increase  in  T^  will,  to  a  considerable 
extent,  determine  the  nature  of  the  temperature  field  In  the  blades 
and  disk  of  turbine  stage  I. 

ar  /  T  \ 

Figure  5  gives  the  generalized  dependence  — '-L  j  — f£_)  for  a,-.  It 

“r.  \  T0  I  5 

is  described  by  the  equation 

ar//5  =  «r..l3.9-0,01irw).  (19) 

The  value  of  a™  in  formula  (19)  is  determined  from  the  graphs 
05 


given  in  Fig.  4. 


Thus,  in  this  article,  using  the  criterial  formulas  which  are 
most  expedient  for  calculating  convection,  we  have  obtained  calcu¬ 
lation  dependences  of  a  on  tH  for  the  examined  axial  turbine.  Prom 
them  it  follows  that  the  intensity  of  convection  in  the  turbine  is 
reduced  with  a  rise  in  t^  and  the  coefficients  of  heat  transfer  on 
the  surfaces  of  the  vane  cascades  are  reduced  from  stage  I  to  III. 
The  value  of  t^  has  an  especially  great  influence  on  reducing  the 
heat-transfer  coefficient  in  the  cooling  channels. 

Prom  the  results  of  calculation  of  heat-transfer  coefficients 
for  the  given  turbine  we  have  found  formulas  by  which  we  can,  rela¬ 
tively  easily,  calculate  the  value  of  the  heat-transfer  coefficient 
for  any  turbine  operation  mode  and  for  any  temperature  tH. 

The  method  used  to  process  the  results  of  this  theoretical 
study  of  convection  heat  transfer  in  a  gas  turbine  can  be  useful 
in  heat  calculations  of  other  types  of  turbines. 
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SOLUTION  OF  THE  PROBLEM  OF  STEADY  HEAT  CONDUCTION  FOR  HOLLOW  COOLED 
NOZZLE  BLADES  OF  GAS  TURBINES 

V.  S.  Petrovskiy 


Let  us  examine  the  dy  heat  problem  for  a  cooled  nozzle 
blade  having  construction  as  shown  in  Pig.  1.  Cooling  air  enters 
the  annular  cavity  of  the  upper  shroud  and  then  the  inner  blade 
channels.  To  increase  the  cooling  effect,  a  deflector  can  be  in¬ 
stalled  in  the  inner  cavity  of  the  blade.  After  cooling  the  blade 
the  air  goes  to  the  rotor  cavity  and  then  to  the  turbine  flow 
channel. 

The  influence  of  the  deflector  on  blade  cooling  can  be  taken 
into  account  by  a  certain  coefficient  k.  Then  a  =  ka' ,  where  o' 
is  the  heat-transfer  coefficient,  calculated  by  known  [1,  2]  cri- 
terial  relationships.  We  can  set  k  =  1  for  cooling  without  a  de¬ 
flector,  k  =  1. 3-1.8  for  cooling  with  a  deflector,  and  k  *  2. 1-2. 6 
for  jet  cooling  using  a  deflector. 

To  facilitate  solution  of  the  problem  and  obtain  finite  expres¬ 
sions  in  a  form  convenient  for  practical  use,  let  us  divide  the 
blade  lengthwise  into  two  approximately  equal  parts  x1  and  x^,  (Pig. 
2).  We  will  seek  the  solution  for  each  part  individually. 

The  initial  equations  for  section  x-^,  considering  the  change  in 
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fig.  2.  HrrwjeMnt  of  the  coordinate 
/itsm  when  stating  the  thermal  prob¬ 
lem,  and  the  approximate  behavior  of 
th-  dependences  tJSfx^)  and  t*(Xj)  from 

aquations  (3)  and  (23):  1  -  experi¬ 
mental  curve;  2  -  exponents  Ct*(0)  * 

*1 

«  (o.so.o.asjt'dj);  t*(ix)  «  i.05  * 

*  t*e  (t*c  -  mean-mass  temperature  of 
the  gas)]. 


temperature  of  the  gas  along  the  blade,  can  be  written  as  follows: 


••  Jfl  y- “  -  *fr  K): 


(1) 


Here 


t>n  (  t|) 
rf.r, 


b  g»p" "  gl~pr  • 
1  XS 


(2) 


(P  and  P  are  the  perimeters  of  the  blade  profile,  washed  by  the 
b  r 

air  and  gasj  S  is  the  blade  section;  G  is  air  flow).  For  section  x-^ 


<,-C,(r+.-.)]> 


(3) 


Solution  of  system  of  equations  (1)  and  (2)  for  t^Cx^)  gives  a  linear 
nonhomogeneous  third-order  differential  equation: 


dx\  c/.t,  <tx\ 


where  A\-k2kAk-,:  B\- (£|— kA)k}.ks. 


Setting 


we  get 


£3 kit-,  (v, )  +A\—y, 


(5) 
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(6) 


Sr-*.  •£— 


•  ’  dx{  *  rfX! 

We  will  seek  the  solution  of  equation  (6)  as  the  sum  of  the  solutions 
of  the  homogeneous  and  nonhomogeneous  parts  of  the  equation,  i.e., 

g-9+u o-  (7) 

For  y  we  can  write  the  characteristic  equation 

r»+ft4r*- Jt,r-JkA«0.  (8) 

For  the  interval  of  change  of  coefficients  k^,  k^,  and  which 

is  possible  under  the  examined  conditions,  the  discriminant  of 
equation  (8)  is  less  than  zero.  Consequently ,  it  has  three  real 
roots.  On  this  basis  the  solution  of  the  homogeneous  part  of  the 
equation  can  be  represented  in  the  form 


Here 


t):  i)-f* 

r v<:os  ^60°-t  ; 

v—  —  |  |/;l:  c  —  aroros  —  ; 

\3 

X  .J  ■  9  /  1  27  '  G  2  ' 


The  solution  of  the  nonhomogeneous  part  can  be  obtained  from  the 
equation 

//o  ~  D\C~’ x*.  (10) 

Successively  differentiating  it,  and  substituting  the  obtained 
expressions  into  equation  (6),  we  find 

O'— 4- (11) 

Substituting  (9)  and  (10)  into  (7),  and  considering  (5),  we  get 
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/J(.v)=(C,e,'‘>  -fC;e'<r-  -f  C^»r>) 


1  ,  ^>i* 


,-t|X, 


•^1 


kiki 


k3k  4 


*i*4 


(12) 


Constants  C2,  and  are  defined  by  the  following  boundary 
conditions : 


when  x1  *  0  i 


when  x. 


W 


^.i(^).fi  — Ai.n* 

<//,(/,  )/</*,  =0. 


In  addition,  from  equation  (1)  it  follows  that  when  =  0 


rf.V| 


*A(0U,  —  £-/r  (O.lr, • 


(13) 

(14) 
i 

(15) 


In  boundary  condition  (13)  t  is  the  steady  equilibrium  tem- 
perature  of  the  surface  of  the  upper  shroud  turned  toward  the  gas 
flow.  Temperature  t  can  be  determined  from  the  heat-balance 
equation  for  the  condition  of  steady  heat  conduction  through  the 
shroud  from  the  cavity  with  the  coding  air  to  the  flow  channel  of 
the  turbine  vane  cascade.  This  equation  for  a  shroud  represented 
in  the  form  of  a  flat  wall  will  be  as  follows: 


(16) 


where  q.^  is  the  convective  heat  flow  from  the  gas  to  the  shroud 
(t0  n  enters  into  the  expression  for  q^),  q2  is  the  conduction  heat 
flow  across  the  shroud,  the  convective  heat  flow  from  the 

cooled  surface  of  the  shroud  into  the  air  flow,  and  is  heat  lost 
to  radiation  from  the  air-cooled  shroud  surface. 


Simultaneous  solution  of  equations  (13)— (15)  gives 

Dr 


Ci  — - 


'c, 

D\j  j 


Dr  Dr 

C.,= — — :  Cs=-£- 


where 


Dc,  —  —  r:er‘'<  I r*ryK i  -~(rl—rl)  AT2-j  (r 

A\  =  {r^'t  —  /< ,  —  (rl -  r{)  AT2  H  —rjer*l>)  Kj 

DCt—(r.:er'1'  —  r,r:A(,  -f  (r]~  rl)  K3-i-(r.ter<,Txcr>t')  Kj 

A».i—  (r;  —  rl)  rxer> '■  -f  (rl—  rl)  r.i +  (r\ -  r\)  rae'*'. . 
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(17) 

(18) 

(19) 

(20) 
(21) 


1 


In  expressions  (18)-(21) 


v', - . -A- i4»;  | 

*^(oU  a3*«-c?a.  f 


(22) 


The  initial  differential  equations  for  section  x2  of  the  blade 
be  t 
pression 


will  be  the  same  as  for  x^.  The  function  is  described  by  the  ex- 


^r('Vj)  —  4lQ)r,—  [^r ^0)j.  MW]  * 

Repeating  the  series  of  discussions  for  section  x2  we  get 

1  ,  />!*'•*•  A2 


<*  tv,)  -  + C,?'*'  -f  C6 «v.|  —  4 


Mi 


(23) 


(24) 


Here 


.  • 

D  mi 

s |*t  ^1^2  ^1*2  ~ 

~':t0)*-Jw 

Constants  C^,  C^3 

and  Cg  can  be  found  from  the 

following  boundary 

conditions : 

when  x2  =  0 

«A(0),,/i/xa*=.0; 

(25) 

when  x2  = 

(26) 

when  x2  *  0 

f.i  r  »  *=  (^1 )  r ,  - 

^7) 

Here  t  is  the  steady  equilibrium  temperature  of  the  surface 
of  the  lower  shroud  turned  toward  the  gas  flow.  It  is  found  in  the 
same  way  as  t  n.  Temperature  t^Z.^)  is  calculated  using  (24). 


Prom  conditions  (25)- (27)  we  find 


where 


(28) 
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Lc,~(cr>:> — c'*'-'  ,V.,  ~ir, — r,)  K  (r  *•'*'>  —  #M  A'6:  (  29 ) 

Kx  -{/■, -r.l A'; -r,c'W.) A'c:  (  30 ) 

Dc, = (e'.i.  -  *'.*-)  -H/-J  -  r,;  A"b;  (  31 ) 

=  -f  (r2-r,;^.'.,  (  32) 

A', =  (33) 

.3,. 

Thus  we  obtain  all  the  required  calculation  expressions  for 
calculating  the  temperature  along  the  blade. 

The  proposed  calculation  method  makes  it  possible  to  take  into 
account  relatively  easily  the  changes  in  gas  temperature  along  a 
vane  cascade.  This  method  can  be  used  to  calculate  the  steady  tem¬ 
perature  state  of  cooled  nozzle  blades. 
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SIMULTANEOUS  SOLUTION  OF  STATIONARY  PROBLEMS  OF  HEAT  CONDUCTION  FOR 
GAS-TURBINE  BLADES  AND  DISKS 


V.  S.  Petrovskiy 

In  this  pa£er  we  examine  stationary  heat  problems  for  the 
shrouded  rotor  blades  and  the  disk  of  a  gas  turbine  with  a  central 
opening,  where  the  profile  of  the  disk  can  be  approximated  by  a 
hyperbolic  function. 

Solution  of  stationary  neat  problems  for  gas-turbine  blades 
and  disks  ha^  been  examined  in  [1-33*  In  this  article  we  aim  to 
show  how  we  can  associate  the  independent  solutions  of  heat  problems 
for  the  blades  and  disks  into  a  single  temperature  function.  In 
addition,  we  show  how  to  formulate  the  boundary  conditions  for 
shrouded  blades  with  a  root  and  for  a  disk  containing  such  blades. 
The  examined  method  for  associating  the  distributions  of  temperature 
in  the  blade  t^(x)  and  in  the  disk  tfl(r)  in  terms  of  the  linear 
temperature  differential  in  the  root  showed  that  it  yields  a 
comparatively  simple  finite  system  of  expressions  for  calculating 
the  temperature  of  the  disk  and  blade  in  the  form  of  a  continuous 
function.  The  need  is  also  eliminated  for  giving  the  temperature 
at  any  point  of  the  disk  or  blade  a  priori,  in  order  that  further 
calculation  be  done  based  on  this  temperature.  The  only  Initial 
data  which  we  use  in  calculations  by  the  proposed  formulas  are  the 
results  of  pre-thermal  gas-dynamics  calculation. 


First  let  us  examine  the  thermal  probler  for  the  blade.  In 
general  form  it  can  be  described  by  the  equation 

~  [>■, (/,')  5  (-V)  W  I’', “  1  (*']  P (-V).  (  ].  ) 

In  equation  (1)  the  values  of  heat  conduction  of  the  blade 
material  (A^),  the  area  of  its  section  (S),  and  the  perimeter  i.P), 
and  also  the  average  heat-transfer  coefficient  (a  )  and  the  stagnant 

B 

gas-flow  temperature  (t#),  depend  on  coordinate  x  along  the  blade. 

Considerable  difficulties  are  involved  in  solving  equation  (1)  with 

substitution  of  functions  A  ,(t  ),  S(x) ,  a  (x),  t*(x),  and  P(x).  It 

•i  f\  b  r 

was  obtained,  using  certain  assumptions,  in  [1]  and  [2]. 

Let  us  examine  the  case  when  A^,  aQ,  and  t*  do  not  depend  on  x. 
The  area  of  the  blade  cross  section  will  be  considered  to  change  by 
-ka/L 

the  law  S  =  SQe  ,  where  k  is  a  constant  which  uefines  the  steep¬ 
ness  of  curve  S(x),  and  L  is  the  blade  length. 

The  function  P(x)  can  also  be  represented  by  an  exponent,  but 
a  more  sloping  one,  since  the  perimeter  of  the  blade  changes  more 
slowly  than  the  area  of  its  section.  Calculations  have  shown  that 
this  makes  it  possible,  when  P(L)/P(0)  >_  0.8,  to  introduce  into  (1), 
in  place  of  function  P(x),  the  value  of  the  average  perimeter  P, 
which  does  not  depend  on  x  and  is  equal  to  the  arithmetic  mean  perim¬ 
eter  of  the  base  and  tip  of  the  blade. 


Here  equation  (1)  assumes  the  form 


_ k  ./I,  <kp  t i.'i (J 

d.x-i  L  . lx  t.iS.j 


0;  (o, 


C2) 


where  t*jj  is  the  temperature  of  the  gas  at  the  level  of  „.ie  average 
blade  cross  section  (Fig.  1).  Equation  (2)  is  solved  by  substituting 
e(k/z)x  =  z: 


G,=-  S,lL  [,!,/, (2i:i,/x/:'-)  r 


f 


where  ]/  I0()  and  K0()  are  modified  Bessel  functions, 

I  /. ,.Sjk ■  2  2 
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Constants  A^  and  are  found  from  the 
boundary  conditions  at  the  blade  base  and  tip. 


jti&u&uiusv  r 


JL 


wm^ 


/ 


Fig.  1.  Dlagran  of  a 
jhroudtd  blade  on  a- root 
(a  daalgr.ates  the  mean 
heat -t  ran  a  Ter  coef¬ 
ficients). 


Let  us  superpose  the  reference  origin  x 
and  the  plane  of  the  blade  base  .section.  Let  U3 
assume  that  the  temperature  along  the  root  changes 
linearly,  while  it  remains  conrtanfc  throughout 
the  lower  shroud,  equal  to  the  temperature  of  the 
blade  base.  Then  the  boundary  condition  with 
x  =  0  can  be  obtained  from  the  heat-balance 
equation  (Fig.  2) 


/,(r})  is  the  temperature  of  the  di 
+  a,  is  the  heat 


(4) 

where  50is  the  heat  entering 

section  SQ  per  unit  time  from  the 
working  part  of  the  blade;  Q,  » 

=  -  tic  (0,  f  is  the  heat  ob¬ 

tained  by  the  surface  of  the  lower 
shroud,  turned  toward  the  gas  flow; 
Q,za\A  —  S„  is  the  heat  removed  from 
section  SQ  across  the  blade  root  to 
the  turbine  disk  (here  « 

sk  when  r  =  r^j  /.|S.rf+ 

removed  from  surfaces  S’  and  S" 

“d  °d 


by  the  cooling  air. 


As  can  be  seen  from  the  last  expression,  the  intensity  of  con¬ 
vection  heat  transfer  depends  not  on  but  on  the  gas  temperature 
at  tiic  butt  section  (tpjjj)*  The  difference  in  the  values  of  t*^ 
anu  takes  into  account  the  drop  in  gas  temperature  at  the 

>)ade  boundary  (when  x  =  C) .  The  gas-temperature  distribution 
v  sion  thus  given,  with  consideration  of  its  lowering  at  the  blade 
,  is  of  rectangular  shape  (.Fig.  3)-  Actually,  the  change  in 
temperature  will  be  smoother.  Above  it  was  noted  that  t*  along 
the  b^ade  foil  is  assumed  to  be  constant,  to  facilitate  solution  of 


equation  (.1).  Thus,  assuming  the  gas- 
temperature  distribution  function  to 
be  rectangular  can  reduce,  to  the 
best  of  our  ability,  the  unfavorable 
role  of  the  approximation  made  when 
stating  the  problem. 

pig.  3.  oiati Lbution  of  t*  along  tin  We  must  also  explain  why  the  second 

blsda,  used  to  state  the  problem. 

term  appears  in  the  equation  for  Q1  . 

“d 

Although  the  area  of  cooling  S"  does  not  correspond  to  coordinate 

ad 

x  =  0,  nonetheless  convection  cooling  of  the  side  surfaces  of  the 
root  occurs  and  it  can  be  taken  into  account  by  this  te  .  In  final 
analysis  it  also  influences  the  value  of  t^(0).  Thus  the  boundary 
condition  for  x  =  0  will  be  refined. 


Unlike  the  boundary  condition  dt  /dx  =  C,  where  x  =  L,  used  by 
many  authors,  we  can  assume  that  the  temperature  of  the  upper  ohroud 
and,  consequently,  the  blade  tip,  is  constant.  Actually  the  upper 
shroud,  having  a  developed  surface,  easily  reacts  to  a  change  in  the 
dynamic  equilibrium  thermal  state,  and  its  temperature  is  determined 
by  the  neat-balanoe  equation,  in  which  the  convective  influx  of  heat 
from  the  gas  to  the  shroud  is  equal  to  the  heat  lost  to  radiation, 
i.e. , 


(6) 


It  is  also  essential  that  the  boundary  condition  dt^/dx  *  0,  when 
x  =  L,  is  completely  unconfirmed  by  experimental  data.  These  data 
indicate  that,  the  temperature  has  a  maximum  value  in  the  middle  of 
the  blade  and  dropa  off  toward  the  base  and  the  tip. 

It  is  not  possible  to  precisely  express  the  components  Q  and 

“a 

Q£,  since  the  first  of  these  is  determined  by  the  heat-cransfer 
coefficient  and  the  gas  temperature,  which  will  differ  at  various 
points  on  the  shroud  surface.  The  value  of  Q£  depends  on  the  tem¬ 
perature  of  the  parts  surrounding  the  shroud  and  the  mutual  ir¬ 
radiation  surface,  which  will  change  as  the  rotor  rotates. 

To  calculate  the  equilibrium  temperature  of  the  blade  in  this 
case  we  can  use  the  equation  of  the  values  of  the  convective  and 
radiant  fluxes  which  are  average  for  the  entire  shroud  surface: 


where 


q-,a -  «,  W  H  -  Tj.  <?.-*;  —  jjgp- .  “o  ~  •  7i 


By  temper 'ture  T^eT  we  mean  uhe  tentative  average  temperature  of 
the  surfaces  of  all  parts  surrounding  the  upper  shroud.  It  will  be 
approximately  150-200°  lower  than  T^q.  For  practical  calculations 
we  can  disregard  and  then,  substituting  q-  and  qg  into  equa¬ 
tion  (7),  we  get  a 

(8) 


The  boundary  condition  assumed  when  x  =  L  is  written  in  the 


form 


Using  boundary  conditions  (5)  and  (9)j  let  us  find  the  partic¬ 
ular  solution  of  equation  (2): 

V-  1C, / .  CX. + 

^5  /  -I  A  \ 


-  [C  J  'fatin'**  :L)~  C,  ,V ?  )}  . 
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Here 


Ct  =0.-.0  —  ?1^2  (— /,,1.,1* 

c2= (/,  (2/ji  ,) — s,  (2»iii)]; 

C,=  i>*/C;;2;n C4=//,(2m/2) 

and 

C5=<?*  ?,/,  (2m,)]-{- 

+  //5  (2m,^-5)  [/Ci  (2m,)  +  e,ff.(2m,)]. 


Thus  we  have  solved  the  first  part  of  the  general  thermal  problem, 
for  the  turbine  rotor. 


Now  let  us  find  the  particular  solution  for  the  differential 
equation  of  stationary  heat  conduction  for  the  turbine  disk.  The 
equation  for  an  elemental  annular  section  of  the  disk  will  be  (Pig. 

4) 


stationary  heat  conduction 
for  a  gas-turbine  disk. 


Qa -  Qi:--Q\~rQi'  (ID 

where  Q,—  Q:~~ |>v.Cn/-2y(rj-!~j(/r  is  the  quantity  of 

heat  propagating  along  the  radius  of  the  disk  in 
1  second  by  means  of  heat  conduction;  Q{  --Qs  = 

=  («,•  (d  |^-4^r)]  -I  -  -  /,i(r)])  2.m/c  is  the  quantity 

of  heat  given  off  by  the  lateral  surfaces  of  the 
annular  section  cf  the  disk  (2irr  is  the  length  of 
an  elemental  ring  of  width  dr;  y(r)  is  half  the 
thickness  of  the  disk  as  a  function  of  r;  a  -(r) 
are  heat-transfer  coefficients  on  the 

(r) 


and  aeij(r) 


lateral  surfaces  of  the  disk,  t^Cr)  and  tB^ 


are  the  temperatures  of  the  air  which  cools  the 
lateral  surfaces  of  the  disk). 


In  further  discussions  we  will  set 
«„(r) Mr)  and 

As  shown  by  comparison  of  the  results  of  calculation  with  ex¬ 
perimental  data,  we  can  consider  approximately  that  a  and  t  do  not 

S  B 

depend  on  r. 
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Substitution  of  -  Q2  and  into  equation  (11)  with 

consideration  of  the  approximations  made  gives 

—  / 4a>  v th)  dr  =--  2a, </s --  O  2ar dr .  ( 12  ) 

rir\  Jr  ) 

By  ag  we  can  also  mean  the  arithmetic  mean  of  the  heat-transfer 

coefficient  between  a  and  a  . 

s3  b4 


After  differentiation  of  equation  (12)  we  get 

..  d-'>i  i  u  ■  du  •">*  r.  _n 
^  rfi*  r  (//■  rfr  '/r 


where  0  =  t  -  t  , 

ft  fl  B 


As  an  example,  let  us  examine  a  disk  having  a  hyperbolic  pro¬ 
file  with  a  central  opening  (y  «  a/r,  where  a  is  a  constant). 
Usually,  the  crown  of  the  disk  is  thicker  than  its  bed,  but  for 
solving  the  heat  problem  we  can  replace  the  stepped  transition  from 
the  crown  to  the  bed  by  an  approximate  monotonic  hyperbolic 
function,  correcting  it  with  a  profile  coefficient  such  that  the 
change  in  the  average  value  of  the  thermal  resistance  of  a  real 
disk  is  the  same  as  that  of  the  present  hyperbolic-profile  disk. 


Equation  (16)  is  one  of  the  particular  cases  of  the  generalized 
Bessel  equations,  and  its  general  solution  will  be 


0J  -  V~:  M Aa (-* J). -5- 5;/-i:a (  '3/?)J  • 


(17) 


Re; burning  to  the  initial  independent  variable  we  have 
'h  —  f A./ 1 .1  m/1 :) 5j/ -i 3 (*^* ') |  • 


(18) 


Constants  A2  and  are  found  from  the  boundary  conditions  when  r  = 
=  r2  and  r  =  r^  (see  Pig.  4) .  For  coordinate  r  •  ^  wt  can  compile 
the  following  heat-balance  equation: 


Q'*-Q\=Q>-  (19) 

Here  Q,— b,  —  Car,/-/;:  is  the  amount  of  heat  fed  within  the  disk 

dr 

per  1  second  from  the  arbitrary  cylindrical  generating  sur¬ 
face  of  a  disk  with  radius  ^  and  width  H  (it  is  assumed  that  the 
isothermal  surfaces  in  the  body  of  the  disk  are  positioned  coaxially); 
3  is  a  coefficient  which  takes  into  account  the  thermal  resistance 
of  the  blade  roots;  O'. Sjt  is  the  heat  which  enters  the  disk 
through  the  above-mentioned  arbitrary  cylindrical  generating  sur¬ 
face  (SH  is  the  cross  section  of  the  blade  root);  n  is  the  number 
of  blades;  <?’  =- a/j ,S'  a  is  the  heat  removed  by  the  cooling  air  from 
the  sections  of  the  surface  of  the  disk  crown  included  between  the 
grooves  of  the  blade  roots  (see  Pig.  4). 


Substituting  the  values  of  Q  ,  Q’ ,  and  Q,f  into  (19)  we  get, 

fl  a  °td 

when  r  =  r^j 


where 


dr 


=?35jTrt' 


(20) 


a-!Sc.n 
_  J  • 

linin' 


/.,2-iri//> 


(20*  ) 


The  boundary  conditions  when-r  =  r2  are  compiled  based  on  the 
following  concepts.  Since  we  are  using  a  hyperbolic  disk  profile, 
the  theoretical  profile  of  its  stepped  part  will  be  too  strongly 
developed  in  the  axial  direction.  This  greatly  increases  the  area 
of  the  central  part  of  the  disk,  air-cooled  with  heat-transfer 
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coefficient  aQ.  To  compensate  for  the  artificially  elevated  heat 
removal  in  the  stepped  section  we  can  stipulate  that  there  is  no 
heat  transfer  along  the  inner  surface  of  the  opening,  and  then  the 
boundary  condition  when  r  =  r2  will  b-- 

</0A/t/r=0.  (21) 


Boundary  conditions  (20)  and  (21)  determine  the  values  of  constants 
A2  and  B2  in  equation  (18),  whose  particular  form  in  this  case  will 
be  as  follows: 


where 


O' 

li 

/z 

rj 

[iMi 

2  s:\ 

—  HI  .r  ■  — 

3  *  j 

-/** 

/  2  i  s' 

l  3  *  . 

c«  -Ml 

irt.rl' 

)  f«yi- % 

•It 

/  .2  3 

-fV-lii  yM/i 

In  equations  (10)  and  (22)  the  constants  q>2  and  <p^  include  the  tem¬ 
perature  differential  along  the  root  At  =  t  (0)  -  t  (r, ) . 


First  let  us  find  0^(0)  and  6fl(r1) : 


M'S  I .1  {~  «»/? /-l ^ y  til jr  j : )  - 

_  ^  i.  m/f'j  /,  3  (y  m  .rf3  *)]  • 

Let  us  introduce  the  designations 

MOi-Cni)5.,  t«-*(sv“'s»J)  /,•  f  i . 

C7-- - ~  ;  ~  >‘«:i — 'J  ’ 


i.  ,S$ 


1 

S,r. 


a  <; "  - —o-5.  n 

n  J  V.  />  _  ft 

Og  * - -  •  — 


ftSj  2>.,S„ 
■j 


o 


:•  my 5 2)  / -i '  m:r* ; j  - 

t  ?  i  ?\  #  ,*  2  ?  2\ 


(23) 

(24) 
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Then,  considering  the  values  of  qp2  ar-d  *Pij  f*1*0111  formulas  (5)  and 
(20),  equations  (23)  and  (24)  can  be  rewritten  in  the  form 

! C:/, .>!,)-:•  C.AC.(2/Mi)1  C:, 1 4- 

+  (C-/; i2m, ) -  CKKn (2/m, '.]  Cj 1  (C7 (25) 
tx — /,  (26) 

Subtracting  (26)  from  (25)  we  get 

/.,-rlCj/j  (-•.'!  ;)-,-C..V;  C~!—  C-,K:  (2m,\}CfC7 

”  ”  1  -  ICVs  (-?'«,)  -  C<Ki (2«,)J  Cr'/^-CJ-'CiC,.,  '  (27' 

The  right  side  of  (27)  does  not  depend  on  x  or  r,  and  therefore  the 
value  of  At  =  t  -  t  can  be  determined  on  the  basis  of  initial  data 
of  gas-dynamics  calculations  before  determining  functions  t^(x)  and 
tfl(r).  This  does  not  contradict  physical  sense,  since  At  in  final 
analysis  depends  on  blade  temperature  when  x  =  0  and  disk  tempera¬ 
ture  when  r  =  r^.  In  other  words,  the  values  of  t/](0)  and  t^r-^) 
are  interconnected  unambiguously  by  conditions  of  heat  transfer  and 
heat  conduction,  which  also  determine  the  nature  of  functions  t^(x) 
and  tfl(r)  themselves.  Thus,  substituting  (27)  into  (10)  and  (22) 
we  get  interrelated  dependences  t^(x)  and  t^(r)  in  final  form.  In 
practice  it  is  advisable  to  substitute  the  numerical  value  of  At 
directly  into  (10)  and  (22). 

Theoretical  calculation  of  blade  and  disk  temperature  using 
these  calculation  functions  begins  with  determination  of  the  heat- 
transfer  coefficients.  For  this  we  must  have  data  from  gas-dynamics 
calculation  and  data  on  air  used  to  cool  the  turbine.  The  next  step 
in  calculating  the  temperature  state  will  be  solution  of  the  inter¬ 
mediate  equations  in  numerical  form:  first  to  find  the  value  of  At 
and  then  to  calculate  the  constants  C^<p2/C^ 

...  n:-r: ')  ..  1  '"'A1) 

- - — „ - and  -j? - - in  equations  (10)  and  (22),  reduced  to 

c,.  >  r,  C(J  V  n 

comparatively  simple  form,  changing  coordinate  x  within  limits  from 
0  to  L,  and  r  within  limits  from  r2  to  r^ ,  we  can  find  the  radial 
distribution  of  temperature  in  a  gas-turbine  rotor  disk  and  blade. 
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Thus,  solution  of  the  simultaneous  problem  of  heat  conduction 
for  a  turbine  disk  and  blades  allows  us  to  calculate  the  temperature 
field  cf  a  turbine  rotor,  using  only  initial  data  from  gas-dynamics 
calculations  for  the  blading  and  hydrodynamic  calculations  for 
cooling. 


The  use  of  the  temperature  different-  *vt  the  root  makes  it 
possible  to  compile  boundary  conditions  for  the  blade  base  and  the 
disk  periphery  such  that  after  colution  of  individual  initial 
equations  for  the  disk  and  blade  we  get  a  single  temperature  function 
for  the  rotor  disk  and  blades. 
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NUMERICAL  SOLUTION  OF  THE  PROBLEM  OF  NONSTATIONARY  HEAT  CONDUCTION 
FOR  A  GAS-TURBINE  ROTOR  WITH  ROOT-SHROUDED  BLADES 

Ye.  Ye.  Denisov  and  V.  S.  Petrovskiy 


The  one-dimensional  nonstationary  process  of  heat  conduction 
for  a  gas-turbine  rotor  with  blades  will  be  described  by  the  follow¬ 
ing  system  of  equations: 


for  the  disk 

6t-(r,\)  d.  6  r  .  ,cMr,T)  I  «.(•')  .  ,  .  N,. 

— : - —  — —  —  re  ;r>  — e— i- - -i-i—  ir,  t — i,  ini; 

01  rz(n  C' [  dr  J  f >.V(r)  * 

for  a  solid  blade 

[/. ,A..  , 

(.V)  OV  I  *  Ox  j  C  { v)  * 


C*T 


(1) 


(2) 


We  will  examin-  a  disk  with  a  central  opening.  The  designation 
of  the  heat- transfer  coefficients  on  the  appropriate  surfaces,  the 
direction  of  the  heat  flow,  the  quantities  of  heat  passing  across 
specific  heat-transfer  surfaces  in  3  second  and  used  in  the  heat- 
balance  equations  when  formulating  the  boundary  conditions  are 
given  in  Fig.  3  on  p.  6. 


To  obtain  equations  convenient  for  numerical  analysis,  let  us 
use  the  feature  of  the  actual  distribution  of  temperature  in  the 
blade  and  formulate,  for  the  system  disk/blades,  two  heat  problems 
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flu  Ylh'  f*i  t  f  '  « t  i  i  f  i»  <■  y  <  »-o\a 


- 


opening 

disk 


root 


Mg.  1.  Position  of  coordinates  and  approximate 
form  of  the  temperature  dependence  for  the  sys¬ 
tem  dlsk/tlades:  1  -  disk  without  opening;  2  - 
disk  with  opening. 


(Pig.  1):  the  first  for  the  disk 


(r2  -  r  - 


r^  and  the  root  sections 


of  the  blades  (0  <  _<  l^)t  and 

the  second  for  the  blade  tip  sec¬ 
tions  (0  <  x2  1  l*)'  Here  the 
total  length  of  the  blade  l 
+  Irs  and  l-y  s  With  such  indi- 


h  + 


2  cmv-t  ~  ^2  • 

vidual  statement  of  the  heat  prob¬ 
lem  we  can  formulate  relatively 
simple  boundary  conditions  for  the 
entire  system  and  obtain  the  function 


t#(x)  which  is  of  decisive  signifi¬ 


cance  for  the  distribution  of  temperature  in  the  blade  in  the  follow¬ 
ing  form: 


for  section  0  _<  Z 

C  (.V,)~  £  \./t) — Cg  ( 3) 

for  section  0  1  >'2  £  ^ 

.  t;(x,)^C(0)-Uxl  (4) 

We  will  solve  the  problem  relative  to  a  gas  turbine,  for  which 
functions  x  in  initial  equations  (1)  and  (2)  can  be  written  thusly: 


o+^'W  — r0)-\  <r,), 

uc  (r)  ---  (/-,)  (r ,  <  r  .  '  r,), 


(0 


■  0  •*<:  X.  < 


(5) 

(6) 

(7) 

(8) 


The  temperature  of  the  cooling  air  t  and  the  average  heat- 

6 

transfer  coefficient  in  the  vane  cascade  will  be  considered  to 
be  constants . 


Let  us  x'ormulate  boundary  conditions  for  the  first  heat  prob¬ 
lem.  We  will  consiuer  the  surface  of  the  central  opening  to  be 
adiabatic.  Then 


i 

% 

l 


\ 

i 

i 

s 

s 


The  time  variable  varies  with-?- 
in  the  limits  0  <  x  <  <®. 


As  the  initial  condition 
let  us  assume  (when  x  «  0) 

/*=  (r,  0)  =/.,<*,.  0)  *=/o-  (17) 

We  will  seek  the  solution 
in  the  rectangles  r2  <  r  <  r-^ 
0  <  t  <  T,  and  0  <  xx  £ 

0  <  t  <  T  (Fig.  2)  which  are 

associated  with  the  appropriate  boundary  conditions. 

Let  us  examine  a  node  grid  formed  by  the  points  of  intersection 
of  two  families  of  parallel  straight  lines.  The  first  family 

|V«,ls«//i4  r:%  (i~- 0,  1,  2 . «),  ; 

.vfW'-,  « / — - 0,  l.  2 . *, =/,'//]; 

the  second  family  (v  --/At,  (i-0,  1,  2....)). 

Step  h  along  axes  x  and  and  time  step  At  along  axis  x  are 
taken  tne  same  for  the  disk  and  the  blade. 

The  nodes  lying  on  the  lines  r  =  r2,  r  =  r-^,  x^  =  0,  and  'x-^  = 

=  wil]  be  boundary  nodes  for  the  examined  problem;  all  others  are 
internal  nodes . 
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u. _ _  hUdt _ 

- 

H —  root 

Fig.  2.  Calculation  diagram  using  tht  grid  method. 


If  some  node  in  the  rectangle  for  the  disk/blade  pertains  to 
function  t^(r,x)  or  t^Oi^x),  respectively,  for  it  we  can  write  the 
difference  which  approximates  the  real  values  of  the  derivatives 
at  this  point.  For 

dtA(r,  x)idr  and  0t^(xu  x)/dxt 


this  is 


— 


(18) 


while  for 
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this  is 


o%(r,  r)/dr:  and  d^fx-..  T).dxt 


6il+'Ji-.2r,“'J>+b‘I-l‘H 


The  derivative  9t/3r  in  node  [i,j]  can  be  represented  in  three 


forms : 


fli/./i _ aW-7— n 

-"-S— ;  -  (21) 

6i/J-n_5i/..-n  (22) 

2ir 

In  accordance  with  ( 20 ) — ( 22 )  we  can  use  three  schemes  for 
approximating  initial  equations  (1)  and  (2).  To  solve  the  posed 
problem  on  a  computer  with  limited  storage,  the  first  difference 
scheme  is  more  convenient,  since  the  initial  condition  gives  the 
value  of  the  approximating  functions  0fl  and  6^  for  all  nodes  of  the 
first  series  (x  =  0).  Using  these  data  we  can  calculate  the  values 
corresponding  to  the  second-series  nodes  (x^  q  +  Ax). 

The  values  in  che  nodes  located  on  the  boundaries  of  the  rect¬ 
angles  can  be  obtained  using  boundary  conditions  (9),  .’’1),  (13), 
and  (1*0.  Then  all  calculation  operations  are  repeated  ±  r  the  next 

series  of  nodes  x.  ,  +  Ax,  etc. 

1,1 

The  use  of  a  difference  scheme  involving  equation  (21)  makes  it 
necessary  to  solve  a  system  of  [(n  +  1)  +  (k^  +  1)]  equations  with 
the  same  number  of  unknowns,  which  is  too  much  of  a  load  on  the 
computer  memory. 

The  difference  scheme  using  equation  (22)  is  unstable.  The 
calculation  errors  in  this  case  accumulate  extremely  rapidly,  and 
after  only  several  time  steps  in  many  nodes  they  become  commensurate 
with  the  value  of  the  approximating  function. 

Using  the  first  scheme  for  approximating  initial  differential 


equations  (1)  and  (2),  we  get  the  following  system  of  difference 
equations: 


where 


•  \  H* 

+»,  +,i  v"  +/i,|)at  +*s,ji, 


fiiL=Lii=L  +  iL;  fj'l _ _ • 

•'/  2A  r/  C,04f, 

/l'l--fl,|/L,i,  (/«**  1,  2 . // -  1  h;=0,  1.  2,...): 


(23) 


‘I0’  — 


where 


-rtf-'1.  (21) 

i  « 

5 1  2//  r.,o.»S/  1 

The  boundary  conditions  must  also  be  represented  in  difference 


form.  For  the  line  r  *  r. 


-• — — or  o|,OJ,«0!,’JI; 


for  lines  r  -  r1  arid  x.^  «  0 


fc. /•••») 

i.i'u-m  ..  l|o./tU  .  .. . 

- . - “Vl'l  1  TTji 

h 

V _ ..  ..  f.io.y  «i..„  f.i'.y  .n  _  . 

-  “*T3'A  *0> 


for  the  line  x^  «* 


(!*••/!.  .  1,1  *!*•*•/ 1  ..  ,  .. 

_ - nO,  or  Vi ’  •  . 


(25) 


(26) 

(27) 


(28) 


The  initial  conditions  for  the  examined  problem  can  be  written 


as  follows: 


/«.(>,  1.  2 . //: 

V,,,n  -  />«0,  1.  2 . 


(29) 


The  problem  is  formulated  such  that  the  rectangular  grid  for  a 
ilsk  on  lira  r  "  r,  has  no  independent  boundary  condition.  The 
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values  of  approximating  function  0fl  in  nodes  belonging  to  line  r  = 

=  are  determined  from  the  associated  boundary  conditions  for 
lines  r  =  r^  and  x.^  =  0.  In  turn,  each  of  the  values  of  the  numeri¬ 
cal  solution  in  the  nodes  on  these  lines  is  defined,  on  the  one  hand, 

by  and  0^°’^  according  to  (26),  and  on  the  other  hand,  by 

and  according  to  (27)  - 


Consequently ,  to  associate  equations  (23)  and  ( 2 h )  by  means  of 
such  difference  relationships,  in  which  there  would  be  successive 
transition  from  known  values  in  the  nodes  on  the  li*  ,  ?  and 

x™  to  unknown  values  of  the  approximating  functit.-n.r-  in  nodes  on 
the  lines  r  =  r1  and  x^  =  0,  we  must  solve  equations  (-'Sj  and  (27) 

for  0Ln  j  J  ]  and  qCOjJ],  we  ge|- 
fl  n  *  * 

<};■  "=--^1 ;  0«°-/,--=4*  ( 30) 

A 

where 


?«tj)  *  //(; 

i:=/6(ryr 


2)-/HrA-r;)-l: 


The  set  of  initial  difference  equations  (23)  and  (24)  and 
boundary  conditions  (25)— (28) ,  where  (26)  and  (27)  are  presented 
in  the  form  of  (30),  allow  us  to  determine,  on  a  computer,  the 
approximate  numerical  solution  of  the  thermal  problem  for  the  disk 
and  the  lower  part  of  the  blades. 


Let  us  formulate  the  problem  for  the  blade  tips  (0  £  <  Z0) . 

The  node  grid  for  the  approximate  solutions  is  also  shown  in  ^ig.  3 

It  is  adjacent  to  the  node  ? r  )  for  the  lower  part  of  the  blades 

(0  <  x,  <  I,),  and  is  formed  uy  the  intersect-'  n  of  the  families  of 
1  1  r-n 

parallel  lines  xlj  =  ih  (i  =  0,  1,  2,  . . . ,  kp)  and  =  JAx . 

The  process  of  heat  conduction  in  this  case  is  described  by 
equation  (2),  which  in  difference  form  will  look  lixe  (24),  except 
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that  i  changes  from  1  to 
k2  -  1.  Functions  a^»  P, 
and  have  the  same  form 
as  for  .the  blade  root  sec¬ 
tion  (x-^),  but  the  numerical 
values  of  P(0),  Sfl(0), 
and  will  differ  Func¬ 
tion  t*  is  represented  by 
parabolic  equation  (JJ). 

As  the  boundary  con¬ 
ditions  we  have :  when 
x2  «  0 

(31) 

0X2 

when  c 

1^2*  T)=/r  (/ 2)  “• 

-  (32) 

-  u  — 

The  last  equation  is 
based  on  the  ideas  that 


Fig.  3.  Results  of  numerical  solution  for  stage  I  the  temperature  of  the 
of  a  gas  turbine:  a  -  experimental  curve  for  the 

steady  process.  blade  tip  is  equal  to  that 

of  the  upper  shroud,  whose 

heating  is  determined  mainly  by  convection  and  radiation  heat  transfer 
with  the  ambient  medium  and  occurs  by  a  law  similar  to  those  for  reg¬ 
ular  heat  regimes . 


In  difference  form,  conditions  (31)  and  (32)  are  written  as 
follows : 

for  line  =  0 

(33) 


for  line  x9  «  Z2 
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(34) 


The  boundary  conditions  for  =  l ^  and  Xj  =  0  do  not  provide 
for  equality  of  the  temperature  functions  on  the  left  and  righuj 
therefore  it  is  possible  that  the  approximate  numerical  solutions 
for  sections  x-^  and  x2>  obtained  for  the  same  moment  of  time  jAt, 
will  differ  somewhat  from  one  another.  This  difference  does  not 
exceed  the  error  of  the  calculation  method  as  a  whole,  and  therefore 
we  can  use  as  the  final  result  one  of  the  two  temperature  values  ob¬ 
tained  for  the  boundary  of  blade  sections  x^  and  x2.  Obviously, 
to  increase  reserve  blade  strength  it  is  best  to  take  the  higher  of 
the  two  adjacent  temperature  distributions. 

Figure  shows  the  results  of  numerical  integration  by  the  pro¬ 
posed  method,  and  gives  data  on  direct  measurement  of  temperature 
for  stage  I  cf  a  turbine  rotor.  We  see  the  good  agreement  between 
the  theoretical  and  experimental  results.  The  theoretical  dependence 
for  a  prolonged  time  (approaching  the  stationary  mode)  is  somewhat 
above  the  experimental  dependence  because  the  calculation  was  done 
for  the  mean-mass  temperature  in  the  vane  cascade  of  stage  I,  t*  = 

=  1176°K,  while  the  experimental  data  were  obtained  at  1133°K. 


With  numerical  integration  by  the  grid  method  of  the  heat  -con¬ 
duction  problem  examined  in  this  article,  we  used  an  explicit  scheme. 
In  this  case  it  is  possible,  since  under  stable  conditions 


at  - - and  at  - — 

dma\.S(.  )  2maxP(i) 


functions  S(x)  and  P(~)  =  rz(r)  do  net  change  rapidly. 


The  integration  steps  were  selected  as  follows.  We  first  gave 
the  step  h  ~  0.15r^.  The  selected  step,  identical  for  the  disk  and 
for  the  blade  (the  problem  was  posed  for  a  disk  of  small  diameter), 
allowed  us  to  reduce  the  time  step  At  to  0.5  seconds.  Trial  cal¬ 
culations  with  three  values  of  At,  from  maximum  to  minimum,  ow'-d 
the  good  convergence  of  the  ^elution. 
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INTENSIFYING  THE  COOLING  OF  THE  LEADING  EDGE  OF  ROTOR  BLADES 
V.  S.  Petrovskiy  and  M.  I.  Tsaplin 

With  internal  cooling  of  gas-turbine  rotor  blades  it  is  not 
always  possible  to  guarantee  a  lowering  of  the  temperature  of  the 
hottest  parts  of  the  blade.  We  know  that  one  of  the  hardest-to-cool 
parts  is  the  leading  edge  which  receives  about  25$  of  the  heat  fed 
to  the  blade. 

At  present,  great  interest  has  arisen  in  so-called  penetration 
cooling.  However,  the  development  of  such  a  cooling  method  has  been 
hampered  by  the  lack  of  porous  materials  having  the  requisite 
strength  and  hydraulic  characteristics.  Therefore,  i  number  of 
organizations  have  studied  cooling  methods  similar  to  penetration 
cooling  -  perforation  of  a  finished  blade  by  electroerosion,  laser, 
or  electron  beam.  The  figure  shows  a  blade  in  which  cooling  of  the 
leading  and  trailing  edges  is  intensified  by  this  method. 

Naturally,  the  operating  conditions  of  a  wall  perforated  with 
a  specific  spacing  and  a  wall  made  of  a  porous  material  are  differ¬ 
ent,  particularly  relacive  to  the  protective  action.  When  air  is 
blown  in  in  the  form  of  individual  jets,  the  decrease  in  heat  flow 
at  the  point  of  intake  is  of  a  clearly  expressed  local  nature  which 
leads,  on  the  whole,  to  a  lowering  of  the  effectiveness  of  the  p^o- 


tective  effect.  However,  with 
exhausting  of  the  air  there  de¬ 
velops  a  noticeable  aftereffect, 
protecting  certain  parts  of  the 
blade  in  the  direction  of  gas  flow. 
Thus  we  can  consider  that  cooling 
by  blade  perforation  is  midway 
between  purely  penetration  cooling 
and  film  cooling.  In  this  regard 
it  is  ad\isable  to  examine  the 


Blade  with  perforated  leading  is  advisabie  to  examine  the 

eage. 

possibility  of  cooling  the  leading 
edge  of  a  hlade  which  combines  the  high  effect  of  internal  cooling, 
characteristic  of  a  permeable  wall,  with  the  protective  aftereffect 
at  the  point  of  inlet.  The  diagram  of  such  leading-edge  cooling  is 
shown  in  the  figure,  which  also  shows  the  location  of  the  coordinate 
system  with  an  indication  of  the  reference  origin  and  the  heat-flow 
scheme. 


Let  us  examine  the  distribution  of  temperature  on  the  leading 
edge  of  a  blade.  Let  us  isolate  a  sector  on  the  inlet  edge  which  is 
bounded  by  radius  of  curvature  rQ  and  lines  connecting  the  coordinate 
origin  with  the  points  of  discharge  of  air  onto  the  blade  surface. 

A  heat  flow  with  known  value  a  acts  on  a  cylindrical  surface  of 
radius  rQ.  The  influence  of  the  discarded  walls  on  the  temperature 
of  the  leading  edge  is  replaced  by  a  heat  flow  of  intensity  q,  uni¬ 
formly  distributed  along  the  boundary  lines,  considering  that  this 
value  is  known. 


Under  such  conditions  the  heat-balance  equation  for  an  elemental 
volume  will  have  the  form 

—  1-  a-Ju.1  4*.  ( 1 ) 

Jr  l  Jr J 

where  t  is  the  current  value  of  the  temperature  of  the  leading  edge 
along  the  radius;  6  is  the  temperature  of  the  cooling  medium;  r  is 
the  current  radius;  F(r)  *  q>rh  is  the  area  of  a  cylindrical  section 
of  an  element  of  the  leading  edge;  h  is  the  thickness  of  the  isolated 
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element,  equal  to  the  distance  between  the  channels  along  the  blade; 
d  is  the  channel  diameter;  aB  is  the  heat-transfer  coefficient  inside 
the  channel  to  the  cooling  air;  5  =  Z/rg  is  the  proportionality 
factor,  which  takes  into  account  the  difference  between  channel 
length  and  radius  rQ. 

Let  us  introduce  the  dimensionless  parameter  /» — and  the 

Ac// 

relative  coordinate  y  =  r/rg. 

Let  us  rewrite  the  initial  equation  in  the  form 

i/  _ m(t—  0)=2</-^  .  (2) 

'  ,l,j1  dtj  \  i  i  ?x 

We  take  into  account  the  heating  of  the  cooling  medium. 

For  the  given  cooling  scheme  we  can  assume  that  heat  transfer 
on  the  inner  radius  of  curvature  of  the  inlet  edge  is  practically 
equal  to  zero.  Then  the  current  value  of  the  temperature  0  of  the 
cooling  medium  can  be  associated  with  the  initial  value  eQ  by  the 
relationship 

(3) 

where  G  is  the  air  flow  arriving  per  element  of  length  of  the  blade. 

Substituting  the  value  of  0  into  equation  (2)  we  get 

//-TTT-vH-  wi  f  - ~  V’ (4) 

ay  V'* 

where  s—~vr—  is  a  dimensionless  parameter.  Let  us  introduce  the 

(iC  p 

new  variable  z  =  sy  and  ret 

S.w-Vr^'.  (5) 

Then  equation  (4)  can  be  rewritten  in  the  form 

0—  (6) 

d  c-  ..  - 


where  is  a  dimensionless  parameter. 
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For  integral  values  of  n,  the  general  solution  can  be  found 
using  exponential  series. 

The  general  solution  for  the  second-order  equation  will  be 

ti  =  C]0n,l(s) 

while  the  first  particular  solution  [1] 


*~0 

Since  for  boundary  conditions  of  the  form 

</0/j2=Bir(/’r  —  /)  when  2=s; 


— =0 
dz 


when  2  =  0 


the  integration  constant  C 2  should  be  zero,  it  is  not  necessary  to 
seek  a  second  particular  solution. 

Determining  integration  constant  C 1  from  conditions  (8)  and 
converting  to  the  initial  variables,  we  can  obtain  a  ca.  ulation 
expression  for  determining  the  temperature  across  the  in  edge 

'"M-  (9> 


Here 


’j  • 


is  the  equivalent  temperature  of  the  cooling  medium; 


n 

«„{(/)=  - - - (si,',*; 

Li  (*!,2 («-*)!  v  •  ' 


$n(l)  is  the  value  of  function  $n(y)  when  y  =  1;  t1  is  the  tempera¬ 
ture  on  the  surface  of  the  edge,  i.e.,  when  v  =  1, 


n:  /•  .  p  1 

•>htw-r  — — 

r- 

8*  (I)  d„ 

where  Bi_  =  u  rr,/A  is  the  Biot  criterion  for  the  leading  edge. 
Thus  we  have  obtained  an  expression  for  calculating  the  temperature 
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of  the  inlet  edge  which  contains ,  in  addition  to  the  known  values , 
the  value  of  q  which  is  still  unknown.  Let  us  find  the  value  of  q 
from  the  condition  that  on  the  blade  wall  there  is  a  protective  film 
of  air  which  passes  through  the  opening  on  the  leading  edge. 

Under  the  assumption  of  a  one-dimensional  scheme  for  heat  flows, 
the  initial  heat-conduction  equation  vrill  be 


Ti  [Lf  17 )  +  a‘H a*H 


where  6  and  f  =  H6  are  the  wall  thickness  and  area;  H  is  the  dis¬ 
tance  between  the  channels  emerging  on  one  side  of  the  blade;  l  is 
the  distance  from  the  point  of  air  exhaust. 

Introducing  tht  concept  of  the  average  temperature  of  the 
washing  medium 

(13) 

where  tg^  is  the  adiabatic  wall  temperature,  and  designating  the 
relative  coordinate  as  x  =  1/8,  we  get  initial  equation  (12)  in  the 


(i4) 

where  k-  —  1?-  is  a  dimensionless  parameter;  u/p—— '-—u  is  the  average 
x  ~ 

heat-transfer  coefficient. 


The  adiabatic  temperature  t  is  equal  to  the  temperature  of 
the  gas  at  the  heat-insulated  wall.  It  depends  on  the  flow  and  tem¬ 
perature  of  the  coola— 1  the  temperature  of  the  gas  flow,  and  the 
distance  from  the  point  of  injection.  Usually,  the  adiabatic  tem¬ 
perature  is  determined  experimentally  and  is  the  basic  factor  which 
takes  into  account  the  aftereffect  of  injection.  The  experimentally 
determined  adiabatic  gas  temperature  at  the  wall  approximated 
quite  well  by  the  formula 

( 15 ) 


Heffe  p^,— •  ,  where  (yw)  is  the  mass  velocity  of  the  gas  flow  at 

2H  ■ 

the  point  of  injection;  GQXJ1/2  is  the  flow  of  coolant  to  one  side 
of  the  blade  profile;  6n  is  the  air  temperature  at  the  outlet  from 
the  cooling  channels;  k  is  a  coefficient  which  takes  into  account 
the  direction  of  injection:  for  tangential  injection  k  «  0.337, 

B 

for  normal  inkection  k  a  0.186. 

8 

Since  T  is  a  complex  function  of  x,  to  simplify  the  calcula¬ 
tions  dependence  (13)  can  be  approximated  by  the  expression 

T  T0  -  <•/*  -  -  7VK1  ( 16 ) 

where  TQ  is  the  initial  temperature  of  the  washing  medium  when  x  **  0; 
T^  is  the  temperature  of  the  hing  medium  as  y  -*•<*>;  6  is  an  approx¬ 
imation  coefficient. 


Setting  AT  =  T  -  t  and  bearing  in  mind  that 

^L=-;:^r.-r0)e-u, 

dx 2 

we  get 


JliHL 

dx 2 


(17) 


(18) 


The  general  solution  of  the  given  equation  is  obtained  as  the  sum 
of  the  general  solution  of  the  corresponding  homogeneous  equation 

{\T)}=Cxc^  -rC,e~tx 
and  the  particular  solution 

Substituting  the  numerical  solution  into  the  initial  equation 
and  solving  it  for  D,  we  get 

& 


(S7  liutu, - 


•  (T  7“0)e_?x. 


Then  the  general  solution  will  have  the  form 

(Ar)= cv** + c,c-«  — (r„  -  ra)  *-»*. 

v  —  Jb« 


(19) 


,v»- 
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The  integration  constant  is  determined  from  the  following  boundary  j  | 

conditions : 


t—  /1)0when  x—0,  ■ 
when x  =  co 


(20) 


where  t^  Q  is  the  average  temperature  of  the  inlet  edge, 

t 

My-  (21) 

o 

Thus,  the  distribution  of  temperature  along  the  blade  profile 
near  the  inlet  edge  can  be  defined  by  the  formula 


/7\-r0)jc' 


—kx  .1.. 


•2 


■  & 


(22) 


Let  us  find  the  heat  flow  from  the  shaped  part  of  the  blade  when 
x  =  0: 


(23) 


where  — T±  is  the  average  relative  temperature  of  the  leading 

T  v  — 

edge . 


Consequently,  the  unknown  intensity  of  the  heat  flow  q  =  Q/2rnh. 

Solving  equations  (9),  (10),  (11),  (21),  and  (23)  simultaneously 
we  can  define  9n  ,  t, ,  and  Q,  corresponding  to  the  given  cooling 

U  j  3  K  1 

1 

conditions.  Setting  /,,0  =  ^  tjyss  *l±i« ,  where  tQ  is  the  temperature  of 

.»  ‘2 
o 

the  leading  edge  of  bhe  blade  at  radius  y  =  0 ,  we  can  write  the 
expression  for  the  average  temperature  of  the  inlet  edge 

1,0  ~  23n(l) 

Setting  . 2'Hk  ■  _ »I».  we  get 
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.‘IS-j.-Ui*! If  - 


and 


Qt  --'iJn  -l— &&  --  2kHk  (T--  T(\  ('  ■ — — — ■*  -  7: .o')  . 
•  »„U)  ‘iy  W-t* 


(26) 


(27) 


(28) 


It  shou  d  be  noted  that  the  second  term  in  expression  (28)  can 
be  both  positive  and  negative,  j.e.,  there  can  be  additional  heat 
removal  to  the  wall  adjacent  t  /  the  leading  edge  or  there  can  be 
heat  feed. 


The  direction  of  the  heat  flow  is  determined  by  the  sign  of 
the  complex 

When  3/ ( B  +  k)  >  t^  q  there  will  be  heat  feed,  while  when  8  t 
i  (8  lc)  <  t^  q  there  will  be  heat  removal.  Thus,  with  an  increase 
in  the  value  of  8  the  ratio  8/(8  +  k)  will  increase,  and  this  means 
that  the  temperature  of  the  washing  medium  as  a  function  of  the  dis¬ 
tance  from  the  point  of  injection  will  change  very  rapidly.  There¬ 
fore,  to  decrease  the  heat  fed  from  the  wall,  in  all  cases  it  is 
necessary  to  strive  to  increase  the  coefficient  k  *  /2Si,  whei*e 

Bi  =  (a  /X)5. 
cp 

Below,  by  way  of  illustration,  we  give  an  example  of  calcula¬ 
tion  of  cooling  the  leading  edge  of  a  rotor  blade  for  the  following 
initial  data: 


8. 


7V*=1 243s  K;  00=  738s  K;  «rti  =  4060  V/mi0K; 
ro=2,25- 10-3n;  X=21  W/m°Kj  ar=210O  W/ni0K; 
d=0,3.|0-3m;  an=2G7  V/ja*°K;  6=  1,5 -  !0~3  m; 
H--li~l0-3  m;  p  =  0,55:  |~0,94;  <f=1.2  rad; 

C*=6,5*10-5  kg/a  (air  flo»*  through  one  channel 
on  the  leading  edge). 


On  the  basis  of  these  initial  daca  we  calculate  the  following 
complexes  which  appear  in  the  calculation  expression: 

Bir=2£uH =0,437:  «  =  — =2,76; 

5=  =  0,0555;  ccn=  .c't-°--  =  1220  w/m2oK; 


^=1/2^-0,418;  <I>  =  J-£*-  =  5.91. 

V  >•  n^a.jrj 


Since  we  obtained  n  =  2.76,  to  calculate  the  values  of  the  functions 
$n(l)  and  d$n(l)/dy  w^  must  round  n  off  to  the  nearest  whole  number. 
When  the  calculation  must  be  precise,  we  must  calculate  the  initial 
values  (i.e.,  t^,  0q  gH,  t1  Q)  for  a  number  of  values  of  n  including 
the  interval  of  interest  to  us,  and  by  interpolation  calculate  the 
true  values  of  t^,  eQ  gK,  and  t1  0,  corresponding  to  the  fractional 
value  of  n.  Let  us  take  n  =  3;  then 


Here  we  assume  that  in  first  approximation  »  eQ  +  30,  since  heating 
of  the  air  in  the  channel?  if  still  unknown. 


Setting  n  =  3,  we  can  calculate,  from  formula  (25),  0  =  654°K, 
from  formula  (26)  T1  •  1054°K,  and  from  formula  (24)  T],  Q  =  908°K. 

Then  let  us  find  *~r°=0,0  and  -^2 - ,\0=  0,569—0,-' ''<i\  i.e., 

r  ,  —7*o  ■  •  ~  k 

we  obtained  the  case  of  heat  removal  to  the  wall  adjacent,  to  the 


leading  edge. 

Then,  using  formulas  (28)  and  (27)  we  get 

Qs=2,02»10-S  w  and  0i=-771*K. 

This  result  corresponds  to  the  case  when  the  protective  after¬ 
effect  of  the  air  discharged  from  the  openings  on  the  leading  edge 
was  estimated  from  the  dependences  for  normal  injection  (at  right 
angles  to  the  incoming  flow).  If  we  consider  that  at  an  injection 
angle  of  less  than  90°  (in  practice  it  can  be  about  60°)  the  pro¬ 
tective  aftereffect  is  intensified  and  will  be  somewhere  between 
the  aftereffect  with  tangential  and  normal  injection  (at  injection 
angles  from  0°  to  30°  the  protective  aftereffect  will  be  identical), 
with  injection  at  an  angle  of  60°  the  protective  effect  can  be  es¬ 
timated  by  the  coefficient  k  ,  which  is  approximately  equal  to  the 
arithmetic  mean  of  the  values  for  normal  and  tangential  injections. 

Considering  the  Influence  of  the  injection  angle  and  taking 
coefficient  8  *  0.376,  we  get  QQ  3K  *  583°K,  T^  »  1031°K,  and  T1  Q  * 
■  985°K,  i.e.,  the  maximum  temperature  of  the  leading  edge  is  re¬ 
duced  from  1054  to  1031°K. 

The  relative  depth  of  cooling  U  =  0.331.  where  t  is  the 

t'r-U 

temperature  of  the  air  behind  the  compressor. 

Thus,  the  depth  of  cooling  is  very  great. 

The  method  of  cooling  the  leading  edge  of  a  rotor  blade,  ex¬ 
amined  in  tnis  article,  is  very  expedient  since  with  the  appropriate 
materials  It  allows  us  to  reduce  the  blade  temperature  and  raise  the 
temperature  of  the  gas  in  a  turbine. 

The  given  calculation  example  shows  that  as  a  result  of  the 
simultaneous  action  of  internal  (in  the  channels)  and  external  air 
cooling  (by  creating  a  low-temperature  protective  film  on  the  blade 


surface),  the  temperature  of  the  leading  edge  is  noticeably  reduced 
compared  with  that  of  an  uncooled  blade. 
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CALCULATING  THE  COOL i NG  OK  GAS-TURBINE  DISKS 
V.  S.  Petrovskiy  and  Me  I.  Tsaplin 


The  reliability  and  longevity  of  gas-turbine  installations  de¬ 
pend,  to  a  considerable  extent,  on  the  efficiency  of  one  of  the  most 
important  parts  —  the  turbine  disks.  The  efficiency  of  the  disks  is 
greatly  influenced  by  the  level  of  permitted  temperatures,  which  is 
guaranteed  by  using  cooling. 


When  calculating  air  cooling  we  must  solve  the  heat-balance 
equation  for  the  turbine  rotor  as  a  whole,  using  dependences  which 
associate  the  heat  flow  with  the  temperatures  of  individual  character¬ 
istic  points  of  the  rotor.  In  particular,  to  perform  such  calcula¬ 
tions  we  must  establish  the  dependence  of  heat  transfer  to  the  disk 
on  the  temperature  of  the  crown  for  the  given  geometric  dimensions, 
the  thermophysical  properties  selected,  and  the  boundary  conditions. 


The  question  of  calculating  the  temperature  in  cooled  disks  has 
already  been  examined  in  the  technical  literature.  For  example,  [1] 
contains  a  survey  of  methods  for  solving  the  heat-conduction  equation 
compiled  on  the  assumption  of  a  one-dimensional  heat -flow  scheme: 


jt 

*/.V 


\'h  —  77  ^  A"  ~ '' 


(l) 


where  F(x)  is  the  area  of  a  cylindrical  disk  section;  \  and  o  are 
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coefficients  of  the  heat  conduction  of  the  disk  material  and  heat 
transfer  from  its  lateral  surface;  t  and  0  are  the  temperatures  of 
the  disk  and  the  coolant;  r  and  x  »  r/rQ6  are  current  and  dimension¬ 
less  radii;  rQ(3  is  the  radius  of  the  disk  crown. 

On  the  basis  of  the  survey  the  conclusion  is  drawn  that  simul¬ 
taneous  consideration  of  variability,  along  the  radius,  of  the  area 
of  the  cylindrical  section,  the  temperature  of  the  coolant,  and  the 
heat-transfer  coefficients  leads  to  a  substantial  complication  of  the 
calculation  expressions,  which  hinders  their  practical  application. 
This  same  work  poses  the  problem  of  developing  sufficiently  accurate 
and  simple  engineering  calculation  methods.  In  particular,  a  number 
oi  simplifying  assumptions  are  proposed  which  reduce  tne  problem  to 
one  of  solving  the  heat-conduction  equation  for  a  disk  of  constant 
thickness  when  the  temperature  of  the  coolant  and  the  heat-transfer 
and  heat-conduction  coefficients  are  constant  along  the  radius. 


Under  these  assumptions,  the  expressions  for  calculating  the 
temperature  distribution  along  the  radius  and  heat  transfer  from 
the  disk  with  boundary  conditions 

l  —  '.jh  when  -v—  I,  )  /  0  • 


dljiix — •  0  when 


•v—  I,  ' 


assume  the  form 


t = o  -f  — — —  /0  (« jck  ( 3 ) 

Q=a2n/7,lf/.H (/o5  —  5^ .  (U) 

where  IQ(n),  I0(nx),  and  I-^(n)  are  Bessel  functions  of  the  first  kind 
of  an  imaginary  argument;  it=|/  is  a  dimensionless  parameter; 


t  6  and  HoS  are  the  temperature  and  width  of  the  disk  crown. 


The  error  in  the  method  of  calculating  temperatures  with  such 
assumptions  is  too  great,  and  therefore  the  proposed  simplifications 
raise  a  number  of  objections. 
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First,  the  disks  of  modern  gas-turbine  installations 3  particu¬ 
larly  transport  Installations,  are  of  complex  shape  due  to  attempts 
to  achieve  uni form- strength  designs.  Here  the  width  of  the  disk 
crown  near  the  root  joints  is  considerably  greater  than  that  of  the 
shank  *'  with  which  the  crown  is  attached  with  smooth  fillets  (Fig. 
1).  Obviously,  approximation  of  such  a  disk  by  a  profile  of  con¬ 
stant  thickness  leads  to  a  substantial  increase  in  heat  transfer. 


Fig.  1.  Disk 
profile. 


rig.  2.  Chang*  tn  relative  haat  transfar 
*  dla*  or  constant  thleknaaa  at  a 
function  oi  paraEatex-  n,  whan  tha  taapar- 
atur*  of  th*  cooling  aadlua  la  replaced 
by  th*  arithmetic  mean  for  various  value* 
of  Kq  -  0c5/t«Ho#X. 


Second,  when  organ¬ 
izing  the  cooling,  the 
relative  flow  of  air  to 
one  stage  is  usually 
0.3-0.5*.  With  such 
low  flows  the  heating  of 
the  air  becomes  notice¬ 
able.  Introduction  into 
the  calculation  of  the 
average  temperature  of 


the  coolant  does  not  allow  proper  correction  of  the 
heat-transfer  value.  By  way  of  illustration.  Fig. 


2  shows  curves  of  the  change  in  relative  heat  transfer  from  a  disk 


of  constant  thickness  when  the  coolant  temperature  is  replaced  by 


the  arithmetic  mean  temperature.  As  we  can  see,  with  a  decrease  in 


flow  the  error  can  reach  40)t,  and  increases  heat  transfer. 


Third,  when  n  >,  10  the  peripheral  parts  of  the  disk  have  the 
greatest  influence  on  heat  transfer.  In  this  case,  use  in  the  cal¬ 
culations  of  heat-transfer  coefficient  values  which  are  average  over 
the  entire  disk  surface  will  also  yield  an  error,  since  the  local 
heat- transfer  coefficients  in  the  region  of  the  crown  will  be  greater 
than  the  average  coefficients. 


Thus,  the  problem  posed  in  [1]  of  developing  a  simple  and 
sufficiently  accurate  method  for  calculating  disk  temperature  can¬ 
not  be  considered  as  having  been  satisfactorily  solved. 


--  ■ 


Let  us  examine  the  problem  of  distribution  of  temperature  in  the 
disk  shown  in  Pig.  1.  Let  us  divide  a  disk  of  complex  shape  into  two 
parts:  peripheral  part  1  adjacent  to  the  crown,  characterized  by  an 
abrupt  contraction  in  profile,  and  central  part  2,  having  a  slight 
taper. 

For  the  contour  which  approximates  the  contracting  part  of  the 
disk,  the  change  in  area  of  the  cylindrical  section 

F(x)  -Fnr,x>\  (5) 

where  Fq6  =  2irro6HoC5  is  the  area  of  the  cylindrical  section  of  the 
crown,  and  p  is  a  certain  constant. 

Equation  (1)  can  be  solved  in  quadratures  only  for  the  case 
p  =  3  and  constant  temperature  of  the  coolant.  This  solution,  for 
boundary  conditions  (2),  has  the  form 

/-on- (6) 
where  m  =  -1  +  /I  +  n2  is  a  dimensionless  parameter. 

For  real  turbine  disks  the  contracting  peripheral  part  cannot 


always  be  reproduced  by  the  contour  p  =  3  and,  in  addition,  the 
temperature  of  the  cooling  air  is  not  constant  along  the  radius. 

Let  us  assume  that  the  temperature  difference  t  -  6  will  change 
along  the  radius  by  an  exponential  law 

t — 0=A/.vs,  (7) 

where  s  is  a  certain  constant. 

Then  the  solution  of  equation  (1)  with  boundary  conditions  (2) 
will  be 

*  *  —  Onx  +  ( /o5 — Ocx)-**.  (8) 

where  k  is  a  dimensionless  parameter,  6  is  the  initial  temperature 
of  the  cooling  air,  and  A6  is  air  heating  due  to  heat  transfer  with 
the  disk. 
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The  parameter 


(9) 

Prom  (8)  and  (9)  it  follows  that  A/<-<oe—  (0«+AO)  and  s  ■  k  4  p  - 
-  3.  In  this  case,  heat  transfer  from  the  disk 

0*=2.'t//oc^  (^oo — Otx)  •  (10) 

Let  us  express  heating  of  the  cooling  air  as  a  result  of  heat  trans¬ 
fer  with  the  disk  as  fellows: 

A  *  =-£-*=?*&*-««),  (11) 

Oef 

where  q—  — — —  is  a  dimensionless  parameter  and  Q  is  the  flow  of 

Gc„ 

cooling  air  on  both  sides  of  the  disk,  and  substitute  it  into  ex¬ 
pression  (9).  After  certain  transformations  we  finally  get 

t.  (>- (12) 

Assumption  (7)  is  valid  when  the  temperature  in  the  center  of 
the  disk  equals  that  of  the  cooling  air.  Such  a  condition  can  be 
strictly  satisfied  only  for  disks  which  are  constricted  toward  the 
center.  The  thickness  of  the  disk  can  also  be  expressed  as  H  » 

*  F  ^x*'"1,  which  means  that 

P  >  2.  (13) 

It  interesting  to  note  that  for  disks  of  constant  thickness, 
for  which  the  parameter  n  is  rather  high,  the  temperature  in  the  cen¬ 
ter  approaches  that  of  the  cooling  air.  When  n  *  5  the  temperature 
difference  is  several  degrees,  while  when  n  >  10  there  is  practically 
no  difference.  Consequently,  under  such  conditions  formulas  (6)  and 
(8)  can  be  used  to  calculate  temperatures  and  heat  flow  in  disks  of 
constant  thickness  as  well. 

To  illustrate  this.  Fig.  3  shows  a  curve  of  the  change  in  Q 
for  p  *  1  and  q  «  0,  the  ratio  of  the  values  of  heat  transfer  in  a 
flat  disk,  calculated  from  formulas  (4)  and  (10).  When  n  >_  10  the 
error  in  determining  the  heat  flow  frem  (12)  is  less  than  5f. 
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This  same  result  can  be  obtained 
by  expanding  the  functions  IQ(n)  and 
I(n).  Then 


Q- 


_3_  15 

8>t  +  1^»2  *  •  * 

_1_  ,  1 
Bn  +  12P/t3 


(14) 


Fig.  3.  Change  In  relative  error  whan 
detamlnlns  heat  transfer  from  a  disk 
it  constant  thickness,  using  simplified 
formulas,  is  a  function  of  values  of 
paraatters  r.  and  q. 


As  can  be  easily  seen,  this  relationship 
tends  toward  one  for  large  values  of  n. 


As  already  noted,  the  central  parts  of  gas-turbine  disks  have 
a  slight  taper.  For  fixed  and  portable  gas-turbine  plants  with  such 
disks,  n  >>  5  as  a  rule.  For  these  cases  the  additional  heat  trans¬ 
fer  along  the  turbine  shaft,  and  also  the  central  opening  in  the  disk 
(where  tnere  is  no  heat  transfer),  have  no  noticeable  effect  on  the 
distribution  of  temperature  along  the  radius  or  on  the  intensity  of 
heat  transfer  from  the  disk.  Consequently,  for  these  conditions,  ac¬ 
cording  to  [2],  we  can  use  the  familiar  dependences  for  a  disk  of 
constant  thickness  derived  relative  to  boundary  conditions  (2). 


Let  us  write  the  distribution  of  temperature  along  the  disk 
radius  relative  to  boundary  conditions  (2)  in  the  form 


Reproducing 


t  —  b 


7.0) 


i:(4 


(15) 


where  x(x)  is  some  function  of  a  dimensionless  coordinate. 


Having  at  our  disposal  solutions  which  describe  the  distribution 
of  temperature  along  the  radius  for  the  simplest  contours,  which 
individually  approximate  the  peripheral  and  central  parts  of  the  disks, 
with  the  given  boundary  conditions  we  can  find  an  expression  for  cal¬ 
culating  the  temperature  in  a  disk  of  more  complex  profile.  For  this 
purpose,  at  the  point  where  one  contour  becomes  another  it  is  neces¬ 
sary  to  combine  solutions  (15).  Designating  by  subscripts  1  and  2 
the  values  in  the  expression  for  the  peripheral  and  central  parts  of 
a  disk,  respectively,  for  the  transition  point  (x  =  5)  we  get 
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the  condition 


•li  |  „  *^3. 


Obviously,  for  t*'  fiven  values  of  tQg  and  00X  conditions  (16) 
are  mutually  exclusive.  In  order  to  observe  equality  of  the  heat 
flows  at  tha  transition  point,  temperature  00Xl  must  be  selected  so 
as  to  assure  the  required  heat  transfer  on  the  part  of  the  periphei’al 
section  when  replacing  the  poorly  approximating  contour  in  the  cen¬ 


tral  p?,  >*t.  We  designate  this  temperature  as  0 


analog 
trodu 
heat  t:- 


BX,3H 


We  proceed 


ly  Wt'-en  selecting  temperature  tQ02  and  correspondingly  in- 

\e  designation  t  *  .  r’hen  the  expression  for  calculating 

00 

?er  o'nd  .  jr  *.*  ?  rature  distribution  along  the  disk  radius 


assumes  ' ne  form 

Q  —  2rt//0&jl(/01-,|  »*)  •  ( - 

The  values  in  formula  (17)  are  the  equivalent  temperature  of  the 
cooling  medium  for  approximation  contour-  1 

•/.!  <!)  /q (?)  ^  -r  V.  (/.;  (;)  •/■". (0  -  •/■! (0  (01  ■  C 

-n  (•;)/.;(:) 


/,  0  J- . 

‘I  »fc  I 


-—-,7"  /!  {.V)  when  1  C 
•/.i  10 


=  **  T  Vi (-v)«hen  5  <  .v  <  0 

7.2  (0 


and  the  equivalent  temperature  on  the  crown  for  approximation  con¬ 
tour  2 

(20) 

Xj  ©Xi«!) 

In  a  number  of  cases,  when  peri’orming  actual  calculations  it  is 
possible  to  simplify  considerably  che  derived  formulas.  For  example, 
when  approximating  the  constricted  peripheral  sections  x^C1)  *  1* 
Then,  for  relatively  thin  disks,  with  high  values  of  the  heat-trans¬ 
fer  coefficients,  the  dimensionless  parameters  for  the  central  parts 
are  quite  large  (n2  >  10).  In  these  cases  00X  3K  and  to6  3K  tend 
toward  0n  and  then  relationships  ( 17)— ( 20 )  degenerate  into  elemental 

B  X 
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formulas  (8)  and  (10).  This  indidates  that  for  large  values  of 
parameter  n,  basic  heat  transfer  occurs  on  the  peripheral  sections 
of  the  disk,  and  the  shape  of  the  central  portions  has  no  substan¬ 
tial  influence  on  it . 


-  e  should  note  that  a  calculation  expression  of  the  type  of 
(6),  which  is  a  particular  case  of  expression  (8),  was  proposed  in 
Stodola's  work  [3],  and  at  that  time  was  widely  used  [4].  However, 
criticism  contained  in  [5]  resulted  in  the  fact  that  the  Stodola 
method  was  generally  relegated  to  obscurity.  This  was  not  quite 
justified.  • 


Figure 
2 


4  shows  curves  of  relative  heat  transfer  from  a  disk 


ir.r*£lt'4-6) 


In  conclusion,  let  us  estimate  the 
error  when  reducing  the  two-dimensional 
heat-conduct] on  to  a  one-dimensional  one. 


Let  us  examine  temperature  distribu¬ 
tion  in  a  flat  disk  under  the  assumption  that  the  temperature  of 
the  cooling  medium  does  not  change  along  the  radius.  The  initial 
heat'-conduction  equation  for  this  case 


0-1 


_ _  j_  dt  .  i  fit 

d.x?  '  x  dx  <V  5 


(21) 


where  x  =  rH,  y  -  h/H  are  dimension] ess  coordinates. 


Let  us  find  the  general  solution  for  the  following  boundary 
conditions : 
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Q/2ur  ot(to6  -  9)  as  a  function  of  parameter  n  for  various  values  of 


p .  With  an  Increase  in  n  the  influence  of 
constant  p  on  heat  transfer  decreases . 
Consequently,  the  Stodola  formula  and  cal¬ 
culation  methods  based  thereon  can  be  used 
in  cases  when  n  ^  i5.  Such  conditions,  as 
a  rule,  exist  for  the  disks  of  aviation 
gas  turbines . 


•^iTrvrir  — ' ri i  -jVwVi VT J ' ■  a? 


Designating 


1)  t  =»  /o4  wh«n  jc=.vf;5v 

2)  _±-=— (/-5'.whenw=I. 

Oy  \ 

cty  i _ ... _ 

Reproduced  from 
best  available  copy. 
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t  Bi*  (/r » •—  *>") 

■*  =  Br  <ycS— 6')  * 


after  the  solution,  presented  in  [6],  we  get  calculation  expressions 
for  determining  the  temperature  field  through  a  flat  disk  and  the 
value  of  heat  transfer  from  the  crown: 

/^/o6-2(/o6-V,,31'-V.  ;  (23) 

rt—l 

Q«4nr06).(45-V)BI'  £  - .  (2^> 


(Bi'  .Ua  —  P,3)  Hr, 

where  Mn—  — ^  p~  ~l  Prt  *fi_  •  p  *  u  /x  ,  is  an  expansion  param- 
/to ‘h/to  +  Bi-ch  ®  m  06  *  * 

eter:  u  is  the  value  of  the  roots  of  a  zero-order  Bessel  function 
ra 

of  a  real  argument. 

In  the  case  of  symmetrical  cooling  conditions  (Bi1  *  Bi”  *  Bi, 
k  *  1),  calculation  expression  (2l})  is  considerably  simplified: 


m 

Q=8xr;,',.(r,t-6) Bi  V  - - ■ 

—  (p-  — Bit  x 


Having  determined  heat  flows  Q2  from  formula  (25)  for  symmetrical 
boundary  conditions,  and  hatting  calculated  their  corresponding  heat 
flows  Q1  from  formula  (4),  we  can  construct  the  graph  *  f  x 

x  (Bi,  xq6)  (Pig.  5).  From  the  graph  we  see  that  reduction  of  the 
two-dimensional  problem  to  a  one-dimensional  one  leads  to  an  error 
which  becomes  substantial  at  large  Bi.  This  error  tends  to  increase 
the  value  of  heat  transfer  from  the  di3k  and  results  in  lower  tem¬ 
perature  values  on  the  crown j  when  Bi  <  0.2  the  error  does  not  ex¬ 
ceed  5%  and  need  not  be  taken  into  consideration. 
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Mg.  5.  Change  in  error  when  ae- 
termlnlng  heat  transfer  from  the 
crown  of  a  disk  of  constant  thick¬ 
ness  when  reducing  the  two-dimen¬ 
sional  heat-conduction  problem  to 
a  one-dimensional  one. 


from  the  crown. 


defined,  as  shown  in  [8] 


In  order  to  bring  into  accord  the  heat 
flows  obtained  from  formulas  (4)  and  (25) » 
it  is  necessary,  when  examining  the  one¬ 
dimensional  problem,  to  perform  the  calcu¬ 
lations  for  a  certain  effective  width  H  A  = 
=  H/n,  where  n  =  Q-^/Qg  a  correction 
factor. 

For  contours  which  constrict  toward 
the  center,  the  correction  for  two-dimen¬ 
sionality  will  be  substantially  smaller, 
since  the  thermal  resistance  nf  such  bodies 
will  tend  toward  its  minimum  value,  deter¬ 
mined  by  the  one-dimensional  scheme  of 
heat  flows  [?]. 

Using  expression  (24)  we  can  also  es¬ 
timate  the  influence  of  asymmetry  of  the 
conditions  for  cooling  of  the  lateral 
surfaces  on  heat  transfer  from  a  disk.  As 
we  know,  when  solving  a  one-dimensional 
heat-conduction  problem,  to  calculate  the 
asymmetry  of  cooling  of  the  lateral  sur¬ 
faces,  by  the  values  t,  0,  and  a  we  mean 
certain  average  temperatures  and  average 
heat-transfer  coefficients,  which  should  be 
,  by  the  formulas 


t 


u't’  —  oV 
‘  a'  — 


(26) 


where  t’,  t",  O’,  0",  a’,  and  a"  are  the  temperatures  on  the  surface 
of  the  disk,  the  temperatures  of  the  washing  medium,  and  the  heat- 
transfer  coefficients  on  one  side  of  the  disk  and  the  other,  re¬ 
spectively. 


Such  a  method  is  approximate,  and  therefore  the  correct  result 
will  be  obtained  only  when  the  temperature  on  the  lateral  surfaces, 
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averaged  by  the  indicated  method,  does  not  differ  substantially  from 
the  average  temperature  across  the  section. 

Figure  6  3hows  the  results  of  calculating  heat  flows  by  formula 
(24)  in  the  form  of  the  ratio  of  Q_,  obtained  for  the  case  Bi*  « 

“  Bi”  *  Bi,  to  the  value  of  Qm  found  for  Bi*  /  Bi”  with  Bi  *  idem 
and  unchanged  temperature  of  the  cc  ling  medium  (6  *  idem).  In  ac¬ 
cordance  with  (26),  Bi  *  (Bi*  +  Bi” )/2 . 

i 

From  the  graph  we  see  that  with  a  decrease  in  the  ratio  Bi”/Bi’, 
i.e.,  with  an  increase  in  the  degree  of  nonuniformity  of  cooling  of 
the  lateral  surfaces,  heat  transfer  from  the  crown  decreases,  where¬ 
upon  its  maximum  deviation  from  the  value  corresponding  to  symmetrical 
cooling  conditions  will  occur  when  Bi"  ■  0  and  Bi*  «  2Bi.  A  decrease 
in  heat  transfer  with  asymmetrical  cooling  depends  on  the  value  of 
Bi  and  also  on  the  relative  disk  radius  xo6,  and  can  attain,  for  real 
disks,  20-30*  of  the  heat  transfer  determined  from  the  averaged  Bi 
values.  Such  asymmetrical  conditions  can  be  found  in  schemes  where 
one  of  the  lateral  sides  of  the  disk  is  shielded  from  the  cooling 
air.  Therefore,  when  calculating  the  heat  transfer  from  a  disk  using 
dependences  of  the  type  of  (17)  it  is  necessary,  for  the  indicated 
cases,  to  introduce  a  correction  in  accordance  with  the  graph  in 
Fig.  6. 

Such  calculations-,  carried  out  for  9*  *  9" ,  Bi’  ■  Bi?  *  Bi,  and 
9  *  (9'  +  6")/2  *  idem,  showed  that  a  chapge  in, the  value  of  (tQ6  - 
-  9")/(to(5  -  9’)  within  limits  from  one  to  zero  leads  to  a  decrease 
in  heat  transfer  from  the  di3k  crown  by  no  more  than  1.5*  compared 
with  the  heat  transfer  corresponding  to  symmetrical  cooling  con¬ 
ditions  with  0.1  £  Bi  £  10  and  1  <  xq6  <  10.  Thus,  nonuniformity  of 
cooling,  caused  by  varying  temperatures  of  the  washing  medium,  does 
not  result  in  noticeable  errors. 

The  simple  dependences  obtained  can  be  used  when  calculating 
the  cooling  of  gas-turbine  rotors,  with  consideration  of  the 
variable  area  of  the  cylindrical  disk  section  and  heating  of  the 
coolant . 
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EXPERIMENTAL  STUDY  OF  THE  TEMPERATURE  STATE  OF  THE  ROTOR  OF  AN 
AXIAL  GAS  TURBINE 

A.  A*  Luzhin,  V.  S.  Petrovskiy,  M.  I.  Tsaplin,  and  V.  I.  Krichakin 


Experimental  study  of  the  temperature  state  of  a  gas-turbine 
rotor  is  of  primary  importance  for  estimating  its  reserve  strength, 
determining  The  service  life  of  the  entire  turbine,  explaining  de¬ 
ficiencies  in  tfce  rotor  cooling  system,  and  discussing  prospects 
for  raising  the  temperature  of  the  gaa  ahead  of  the  turbine.  Direct 
thermometry  allows  us  to  refine  a  vast  number  of  calculation  charac¬ 
teristics  used  as  the  basis  for  designing  the  turbine.  Besides, 
accumulation  of  information  on  the  actual  temperature  fields  in 
real  gas-turbine  designs  allows  us  to  improve  the  calculation  methods 
for  studying  the  temperature  states  of  the  rotors,  and  develop  new 
theoretical  methods  for  studying  boundary  conditions  in  operating 
turbines . 

In  this  article  we  give  the  results  of  thermometry  of  the  rotor 
of  a  small,  three-stage,  axial  gas  turbine. 

The  mean-mass  temperature  of  the  gas  ahead  of  the  turbine  during 
the  tests  reached  1100°K;  the  flow  of  air  for  cooling  the  rotor  was 
no  more  than  5%  of  the  gas  flow  through  the  turbine.  Air  for  cooling 
the  turbine  rotor  is  taken  at  the  compressor  outlet  and  its  tempera¬ 
ture  under  rated  conditions  *»  520-530°K. 
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To  reduce  the  temperature 
differential  along  the  disks  of 
the  first  stages,  in  the  lower 
parts  of  the  rotor  blades  we 
created  unique  heat  chokes 
in  the  form  of  cooling  channels. 
The  rotor  blades  of  all  three 
turbine  stages  have  upper  and 
lower  shrouds.  The  cooling 
channels,  as  can  be  seen  from 
Fig.  1,  are  formed  by  the  blade 
root,  the  lower  blade  shrouds, 
and  the  outer  surface  of  the 
disk  crown. 

Air  from  the  compressor  is 
fed  to  cool  the  turbine  through  an  annular  gap  between  the  flame 
tube  and  the  support  housing,  and  goes  to  cavity  I.  From  here  it 
is  directed  through  labyrinth  seals  to  the  rear  bearing  ind  to  the 
peripheral  part  of  the  stage  I  disk  crown,  and  also  through  6  open¬ 
ings  in  the  deflector  to  the  lateral  surface  of  the  disk  and  then 
through  central  openings  in  the  stage  I  disk  to  stages  II  and  III  of 
the  turbine. 

Air  enters  the  cooling  channels  c<f  stage  I  from  cavity  II, 
passing  through  the  labyrinth  seal;  only  a  certain  part  of  the  air 
enters  the  channels.  The  rest  passes  through  the  annular  gap  to  the 
turbine  blading  and  creates  unique  film  protection  at  the  bottoms  of 
the  blades. 

The  upper,  peripheral,  parts  of  the  rotor  blades,  just  as  the 
roots,  will  be  in  a  zone  of  reduced-temperature  gas  flow  since  the 
air  which  cools  the  stator  is  mixed  vrith  the  gas.  Thus,  the  middle 
sections  of  the  blades  undergo  the  greatest  heating. 

Since  the  distribution  of  gas  temperature  with  height  of  the 
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rotor,  of  the  thsracaatry  points.  The  Roam  nuatrals 
designate  tne  ssvitlss  which  tht  cooling  sir  enters . 

1  -  gap  between  flu*  tub*  and  transfer  ring  of  the 
upper  shroud  of  tha  first  neasl*  cascade ;  2  -  point 
of  location  of  the  tharoocoupla  connecting  atrlps 
for  Manuring  tj;  3  -  gep  between  flaae  tub*  and 
lower  shroud  of  first  nozel*  cascade ,  4  -  annular 
silt  through  which  cooling  air  passes  froa  tha  coa- 
pressor  to  the  rotors  5  -  openings  In  tha  deflector 
on  the  leading  aid*  of  the  disk  of  stag*  I  for 
feeding  cooling  air  to  stags*  II  and  III. 


r  ■  *  a 
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van*  cascade  determines,  to  a  considerable  extent,  the  temperature 
field  in  the  turbine  rotor  blades  and  disks,  the  results  of  rotor 
thermometry  can  be  correctly  explained  only  if  we  simultaneously 
measure  the  gas  temperature  at  various  sections  of  the  rotor  blades. 
This  was  taken  into  consideration  when  setting  up  the  experiment. 

The  following  basics  were  used  to  prepare  the  turbine  for 
thermometry . 

The  turbine  rotor  is  the  most  stressed  part,  and  therefore  we 
provided  for  temperature  measurements  at  all  of  its  most  important 
points,  from  the  standpoint  of  strength.  Basically,  these  were 
those  points  at  which  we  can  obtain  data  on  the  temperature  dif¬ 
ferential  along  the  disk  radius  from  the  periphery  to  the  center, 
through  the  disks,  and  in  the  blade  firtree  Joints.  It  was  also 
important  to  measure  the  temperature  of  the  blade  foil  at  various 
sections  along  its  height. 

Besides  the  turbine  rotor,  when  conducting  the  experiment  as 
a  whole  we  prepared  for  thermometry  the  following:  the  blades  of 
the  nozzle  cascades,  the  front  and  rear  turbine  shaft  supports,  the 
stator  shrouds,  and  a  number  of  other  items.  Below  we  will  examine 
questions  pertaining  only  to  the  thermal  state  of  the  rotor. 

To  obtain  the  required  results  from  study  of  cooling  effective¬ 
ness  it  was  necessary,  in  addition  to  measuring  the  rotor  tempera¬ 
ture,  to  determine  the  temperature  of  the  gas  at  the  turbine  inlet 
(t^),  the  temperature  field  of  the  gas  flow  along  the  height  of  the 
vane  cascades  of  all  three  stages ,  and  the  air  temperature  behind 
the  compressor  (tip  and  at  the  compressor  Intake  (tjj). 

In  addition,  it  is  necessary  to  know  the  air  pressure  at  the 
compressor  intake  and  outlet  (Pjj  and  p^)  and  the  gas  pressure  at  the 
turbine  inlet,  and  outlet  (p^  and  p^). 

The  positioning  of  the  thermometry  points  on  the  rotor  also  had 
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to  be  associated  with  data  from  a  theoretical  preliminary  determina¬ 
tion  of  the  temperature  field  in  the  turbine  disks  and  rotor  blades 
to  make  possible  later  comparison  of  the  results  of  theory  and  experi¬ 
ment.  The  final  scheme  for  positioning  the  thermometry  points  on 
the  rotor  is  shown  in  Pig.  1. 

Measurement  of  t^  was  carried  out  in  two  steps .  In  addition  to 
the  thermocouples  on  the  rotor,  12  thermocouples  were  placed  at  the 
turbine  inlet  around  the  periphery,  on  the  average  diameter  of  the 
blading.  To  measure  the  distribution  of  gas  temperature  along  the 
Cn  vane  cascades  [t*(x>]  we  first  installed,  in  the 

p  - lJ_ —  nozzle  cascades  of  al3  three  stages,  thermocouples 

-  r  .-I..-  on  one  nozzle  blade  of  each  stage,  as  shown  in 

— «  Fig.  2.  Processing  of  the  results  showed  that  it 

""  3  is  impossible  to  determine  t*(x)  from  one  blade 

E _ per  stage,  since  the  gas  temperature  is  nonuniform 

_ _  *  around  the  turbine  blading.  Therefore,  to  measure 

.  —  —  t*(x)  we  performed  additional  thermometry  for  which, 

besides  the  12  thermocouples  on  the  average  diam- 

?iB.  a.  Dugrm  of  ih»  eter,  we  placed  connecting  strips  of  5  thermocouples 

location  of  th«  thtrmo- 

coupiaa  on  the  nosiie  each  at  six  locations  around  the  periphery  ahead  of 

blade  (the  nuabers  are 

tionifofBnot*ie°eae-*c~  the  nozzle  cascade  of  stage  I;  then  in  the  nozzle 
cascade  of  each  stage  we  set  up  6  blades  with 
thermocouples  positioned  In  shown  in  Fig.  2.  The  resulting  distri¬ 
bution  of  T*(x)  at  the  outlet  from  the  flame  tube  and  in  the  nozzle 
cascade  of  each  stage,  shown  below  in  Figs.  4  and  5,  was  constructed 
from  points,  each  of  which  is  the  arithmetic  mean  of  the  readings  of 
6  thermocouples  placed  on  the  same  diameter  around  the  periphery  of 
the  blading. 


8  *"  — 


We  must  mention  the  following  relative  to  the  thermocouples  on 
the  nozzle  blades.  Since  the  bead  of  the  thermocouple  was  1.5  mm 
above  the  blade  surface,  it  was  in  the  boundary  layer  of  the  gas 
flow.  Measurement  of  the  temperature  of  the  boundary  layer,  combined 
with  the  damping  exhibited  by  the  thermocouple  bead  itself,  assured 
measurement  of  T*,  to  a  certain  extent.  As  additional  verification 
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sf tolled,  the  error  of  this  method  when  measuring  stagnation  tempera¬ 
ture,  did  not  exceed  the  errors  arising  from  other  causes  (limited 
accuracy  of  the  thermocouple  itself  and  the  recording  instrument, 
due  to  heat  transfer),  and  was  within  acceptable  limits. 

The  air  temperature  at  the  compressor  outlet  was  measured  with 
six  thermocouples.  This  «as  a  sufficient  number,  since  t£  varies 
only  insignificantly  around  the  periphery  of  the  compressor  blading. 

Four  thermocouples .were  installed  at  the  compressor  inlet  ahead 
of  the  protect! fO  cascade  to  measure  tg.  For  greater  experiment 
reliability,  ttlfe  turbine  rotor  contained,  as  shown  in  Fif..,  1,  dupli¬ 
cate  thermocouples  on  the  same  radii  as  the  main  thermocouples  but 
shifted  120°  and  180°  in  a  plane  perpendicular  to  the  turbine  axis. 

An  Indexing  System,  which  facilitated  operations  during  the 
experiment ,  was  developed  for  the  entire  thermocouple  system. 


To  measure  the  temperature  of  the  gas  and  the  rotor  parts  we 
used  Chrome 1-Alumel  thermocouples,  while  to  measure  t|  and  tH  we 
used  Chromel-Kopdl  thermocouples  of  class  1  accuracy.  The  diameter 
of  the  thermoelectrodes,  depending  on  the  measurement  site,  was  0.2- 
0.5  mm.  On  the  rotor  blades  and  disks  the  thermocouple  beads  were 
welded  to  the  surface,  while  the  therraoelectrode  leads  in  the  heat- 
resistant  insulation  were  covered  with  0.1-mra  thick  refractory  steel 
foil,  spot-welded  to  the  surface  of  the  blades  or  disks.  Figure  3 

shows  a  rotor  prepared 

/$$$'£*■  for  thermometry .  The 

s'  ••  ':-,u  v'; foil  assures  the  necess- 

f-  \V  V  ■  i ary  strength  and  protec- 
GZZZZ f.  . :! tlon  of  the  thermocouple 

"'O'*  ^jit^  against  the  action  of  the 

Fig.  3.  Turbine  rotor  disks  prepared  for  gas»  511(1  alao  guarantees 
thermometry.  measurement  accuracy. 


Fig.  3.  Turbine  rotor  disks  prepared  for  6  *  6 

thermometry.  measurement  accuracy, 

since  it  shields  the  thermocouple  beads  from  direct  action  by  the 

washing  medium. 
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The  thermoelectrodes  from  the  thermocouples  were  braided 
together  and  fed  to  a  mercury  slip  ring  attached  to  a  bracket  on 
the  exhaust  pipe.  Stage  III  of  the  turbine  was  used  to  drive  the 
slip  ring. 

The  mercury  slip  ring  is  intended  to  transmit  an  electrical 
signal  from  sensors  mounted  on  the  rotating  parts  to  fixed  recorders. 
It  has  22  pairs  of  mercury  contacts.  Each  pair  of  contacts  is  an 
individual  section;  the  moving  contact  is  located  on  the  rotor  shaft 
while  the  fixed  contact  is  imbedded  in  the  slip-ring  stator.  All 
sections  of  the  contacts  are  insulated  from  one  another  by  textolite 
washers . 

The  thermocouple  leads  fed  to  the  slip  ring  are  soldered  to  the 
contacts  on  a  sliding  insulating  ring  placed  on  an  intermediate 
shaft.  The  contacts  of  the  sliding  ring,  in  turn,  are  connected  to 
the  sliding  mercury  contacts  of  the  slip  ring.  The  intermediate 
ring  was  introduced  into  the  design  to  facilitate  resoldering  of 
the  thermocouples. 

The  cold  junctions  of  the  rotor  thermocouples  were  located  on 
an  insulating  (textolite)  washer  which  rotates  inside  a  recess  in 
the  bracket  which  holds  the  slip  ring.  A  jet  of  air  is  directed 
into  the  recess  to  balance  the  temperature  of  the  cold  junctions. 

The  temperature  of  the  junctions  varied  depending  on  the  injection 
rate,  the  temperature  of  the  ambient  air,  and  other  factors.  The 
result  is  recording  not  of  the  true  value  of  the  temperature  at  the 
investigated  point  but  the  difference  in  temperatures  of  the  hot 
and  cold  thermocouple  junctions.  To  obtain  the  true  temperature  it 
was  necessary,  when  processing  the  measurement  results,  to  add  the 
cold- junction  temperature  to  the  instrument  readings.  This  tempera¬ 
ture  was  considered  equal  to  that  of  the  injected  air,  and  was  mea¬ 
sured  by  a  separate  thermocouple. 

A  total  of  90  thermocouples  were  installed  on  the  rotor.  Since 
the  number  of  contacts  allowed  for  the  simultaneous  use  of  no  more 
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11  thermocouples,  some  of  them  were  periodically  resoldered  in 
ordpr  to  compare  the  measurement  results  and  compile  them  into  a 
single  temperature  field  such  as  would  be  obtained  if  all  thermo¬ 
couples  worked  simultaneously. 

As  the  measurement  Instruments  for  the  turbine  rotor  we  used 
single-point  high-speed  potentiometers,  accuracy  class  0.5,  with 
continuous  recording.  To  measure  the  cold-junction  temperatures 
tH,  t|,  and  t$  we  used  multiple-point  potentiometers  with  complete 
carriage  coverage  pf  the  scale  in  2.5  seconds.  The  high  speed  of 
the  carriage  of  the  discrete  potentiometers  and  the  great  number  of 
polpts  for  measprlpg  each  of  the  indicated  temperatures  made  it 
possible,  Khen  processing  the  measurement  results,  to  clarify  the 
trend  of  the  temperature  dependences  during  transitional  regimes 
as  well  (start-up,  stop). 

To  measure  the  distribution  of  gas  temperature  along  the  vane 
casqade  we  used  multiple-point,  discrete-recording  potentiometers. 

The  readings  for  identical  thermocouples  were  recorded  less  fre¬ 
quently  than  with  thermometry  of  the  rotor  but,  since  we  studied 
only  the  steady  state  of  the  turbine,  the  number  of  repetitions  was 
quite  satisfactory  for  obtaining  precise  results. 

The  correctness  of  the  instrument  readings  was  systematically 
monitored  using  an  additional  potentiometer  coupled  in  parallel  with 
all  the  others  to  measure  the  rotor  temperature.  If  necessary,  the 
signal  of  the  thermocouple  of  any  potentiometer  could  be  fed  to  this 
instrument  using  a  switch. 

The  switch-on  of  all  instruments  was  synchronized,  using  one 
start  button  on  the  operator's  control  panel. 

The  steady  turbine  operating  mode  under  various  loads  was  of 
particular  interest.  At  the  same  time  it  was  also  important  to  study 
the  nonsteady  temperature  states  of  the  rotor. 


.1. 


During  3tart-up,  the  turbine  operates  under  severe  thermal  con¬ 
ditions,  under  which  we  can  observe  high  temperature  differentials 
on  the  disks,  jumps  in  the  thermal  stresses  in  the  joint  section  of 
the  disks,  and  other  unfavorable  phenomena.  Study  of  nonsteady 
temperature  states  should  answer  the  question  of  the  number  of 
start-ups  which  limit  reliable  turbine  service  life. 

These  then  are  the  basic  concepts  for  compiling  a  graph  of 
turbine  operation  during  ther.iometry .  The  entire  period  from  start¬ 
up  to  shutdown  of  the  engine  was  divided  into  the  following  stages . 

1.  Start-up,  idle,  warmup,  and  measurement  of  all  parameters 
characteristic  of  turbine  operation. 


2.  A  gradual  increase  of  the  load  on  uhe  shaft  with  a  rise  in 
temperature  from  t*  to  t^max  in  the  intermediate  mode.  In  both  the 
intermediate  and  maximum  modes  the  measurements  were  made  after  a 
slight  delay,  for  the  turbine  to  achieve  the  steady  temperature 
state. 


3.  A  gradual  decrease  in  load  on  the  shaft  and  a  corresponding 
drop  in  gas  temperature  ahead  of  the  turbine  to  idle  through  an 

•J 

intermediate  load,  with  temperatures  measured  in  each  regime. 

The  results  of  all  measurements  were  obtained  in  the  form  of 
graphs  or  a  group  of  points  on  the  recorder  tapes.  Each  group  of 
tapes,  obtained  for  one  day's  testing,  wa3  processed  individually. 
All  recorder  tapes  for  one  day  were  coordinated  relative  to  the 
reference  point  and  the  check  marks  which  were  made  periodically, 
during  turbine  testing,  on  the  tapes  of  all  instruments  by  simul¬ 
taneously  shutting  them  off.  On  the  coordinated  tapes  we  noted  the 
steady  portions  and  the  data  obtained  were  compiled  ; nto  tables. 
Nonsteady  sections  were  noted  only  after  the  first  start-up,  when 
the  entire  engine  was  still  cold  and  its  temperature  was  equal  to 
that  of  the  ambient  medium  ^  The  discrete  values  of  t|,  t*, 

and  tjj  were  averaged,  and  from  them  we  constructed,  on  the  tapes. 
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the  continuous  dependences  t|(x),  tjj(x),  and  tH(x).  The  time  of  the 
nonsteady  regime  was  reckoned  from  the  moment  of  start-up  using  the 


tape  showing  the  temperature  of  the  cold  junction,  which  thus  served 


as  a  unique  time  scale. 


flS.  4.  D#p«nd«nc»  V*(x),  aoiurad 
tulnc  th«t»oco\»pl«*  «t  tht  exit  free 
tlN  riamt  tuba.  <*uh  point  lx  th« 
inrift  of  tha  ratting*  of  (lx 
thamocouplai  positioned  around  the 
blading.  1  -  Idle,  tfi  «  +30«2j  2  - 
Intaroadiata  reglae,  ty  «  +30’C; 

3  -  MXlavai  rtflu,  t,j  «•  *30*0;  k  - 
xmxIbub  ragiaa,  tj  *  +50°C. 

Figures  4  and  5  give  the 
resulting  graphs  of  the  change 
in  gas  temperature  with  height 
of  the  turbine  blading.  Mea¬ 
surements  were  made  in  three 
regimes:  idle,  with  an  inter¬ 
mediate  load,  and  with  maximum 
load  on  the  shaft  (with  maximum 
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mean-mass  temperature  of  the  gas)  at  t^  *  +30  and  +50°C.  Each  point 
cn  these  graphs  is  the  average  reading  from  six  thermocouples  posi¬ 
tioned  about  the  periphery  of  the  blading. 


Simultaneous  thermometry  of  the  gas  flow  at  four  points  of  the 
turbine  blading  made  it  possible  to  get  a  clear  picture  of  the  change 
in  T*(x)  from  the  outlet  from  the  flame  tube  to  the  nozzle  cascade 
of  stage  III.  Associating  the  obtained  dependences  with  the  turbine 
design  features,  it  is  easy  to  explain  the  reason  for  certain  anom¬ 
alies  in  the  gas-flow  temperature  field  end  to  explain  their  effect 
on  turbine  characteristics. 
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The  lowering  of  the  temperature  in  the  peripheral  and  root 
sections  of. vane  cascades  is  explained  only  by  the  admixture  of 
cooling  air.  This  is  indicated  by  the  comparatively  slight  dis¬ 
tortion  of  the  gas  temperature  field  at  the  flame-tube  outlet.  A 
certain  reduction  in  temperature  in  the  peripheral  region  can  be 
explained  here  by  the  flow  of  air  from  the  region  around  the  stator  - 
into  the  turbine  blading  through  the  gap  between  the  flame,  tube  and 
the  inner  surface  of  the  transfer  ring  of  the  upper  shroud  of  the 
stage  I  nozzle  cascade  (see  Fig.  1). 

The  intensity  of  flow  through  the  gap  increases  with  an  in¬ 
crease  in  the  pressure  differential  between  the  air  cavity  and  the 
inlet  to  the  stage  I  nozzle  cascade.  This  is  observed  when  there 
is  an  increased  load  on  the  turbine  (an  increase  in  shaft  load)  and 
is  confirmed  by  deformation  of  the  temperature  field,  which  is  evi¬ 
dent  in  Figs.  4  and  5  for  stage  I  upon  transition  from  idle  to  maxi¬ 
mum  take-off.  A  reduction  in  temperature  of  the  gas  of  the  peri¬ 
pheral  region  at  the  outlet  from  the  flame  tube  and  the  stage  I 
nozzle  cascade  does  not  appreciably  influence  turbine  performance, 
since  it  has  an  insignificant  effect  on  the  mean-mass  gas  tempera¬ 
ture  and  it  corresponds  approximately  to  the  calculated  temperature. 

A  slight  reduction  in  gas  temperature  is  also  noted  in  the  root 
section  of  the  first  nozzle  cascade  (see  Fig.  5).  This  can  be  ex¬ 
plained  by  the  leakage  of  air  into  the  turbine  blading  at  the  point 
where  the  flame-tube  outlet  and  the  inner  shroud  of  the  stage  I 
nozzle  cascade  meet. 

Distribution  of  gas  temperature  at  the  inlet  to  the  stage  II 
nozzle  cascade  is  deformed  in  both  the  peripheral  and  the  root 
sections  of  the  cascade.  Cooling  in  the  peripheral  region  is  ex¬ 
plained  mainly  by  the  influence  of  the  air  which  reaches  the  blading 
through  the  gap  along  the  outer  perimeter-  of  the  flame-tube  outlet, 
and  then  by  the  addition  of  air  which  cools  the  peripheral  shroud 
of  the  stage  I  nozzle  cascade. 


'  A  lowering  of  the  temperature  in  the  root  •,  section  is  explained 

*  *  f 

by  the  feed,  to  the  blading,  of  air  used  to  cool  the  rotor.  The 
trend  of  curves  T*(x)  for  the  nozzle  cascades  of  stages'll  and  III  ' 
attests  to  the  fact  that  the  basic  rotor-cooling  air  is  expended, 
in  stage  X. 

The  nature  of  the  dependences  T£(x)  for-  all  three  stages  allows 
us  to  conclude  that  the  real  change  in  mean-mass  T*  along  .the  tur¬ 
bine  stages  corresponds  to  the  calculated  change. 


An  explanation  is  required  for  the  /certain  discrepancy  in  gas 
temperatures  at  the  flame-tube  outlet  and  at  the  inlet  to  the  stage . 
X  nozzle  cascade,  although  the  maximum  ga3  temperature  in  both 
cases  should  have  been  practically  identical  or  somewhat  higher  at 
the  flame-tube  outlet.  But  the  measurement  results  show  Just  the 
opposite:  the  gas  temperature  at  the  inlet  to  the  first  nozzle  cas¬ 
cade  is  somewhat  higher  than  at  the  flame-tube  outlet.  This  can  be 


explained  by  the  nonuniformity  of  the  gas  temperature  around  the 
turbine  blading  and  by  a  slight  displacement  of  the  thermocouple 


connecting  strips  relative  to  the  stage  I  nozzle  cascades; 


The  results  of  measurements  of  rotor  temperature  in  steady 
regimes  are  shown  in  Fig.  6  as  functions  t^(T^)  and  t^(T^).  Here 
T^  is  the  mean-mass  temperature  at  the  inlet  to  the  first  nozzle 
cascade,  obtained  after  processing  the  results  of  gas-temperature 
measurements  by  12  thermocouples  simultaneously  with  thermometry 
of  the  turbine  rotor.  Each  line  in  Fig.  6  was  constructed  from  the 
readings  of  a  specific  thermocouple  during  operation  of  the  turbine 
in  various  regimes.  Naturally,  during  the  tests  did  not  exceed 
the  maximum  permitted  value.  All  the  lines  are  extrapolated  to 
the  region  of  high  values  of  T^.  We  can  assume  that  such  extrapo¬ 
lation  is  permissible  for  a  specific  turbine  design.  In  this  par¬ 
ticular  case,  with  an  increase  in  T^  the  nature  of  the  function 
T*(x)  for  the  entire  turbine  blading  will  hardly  change,  since  the 
distribution  of  the  cooling  air  over  the  stator  and  rotor  does  not, 
in  all  prob?'  lity,  depend  on  T|.  Therefore,  with  an  increase  in 
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Fig.  6.  Dependence  ct  rotor  temperature  on  gas  temperature  at  the 
turbine  Inlet  T?  (the  lines  correspond  to  the  points  of  temperature 
measurement  as  shown  in  Fig.  1.  The  rcints  on  lines  3,  9,  lo,  19, 

28,  and  29  designate  temperatures  for  which  we  calculated  the  re¬ 
serve  strength  of  o  blade  in  a  given  section  for  rated  turbine  rpm 
(see  Fig.  10). 

in  the  similarity  of  the  temperature  fields  in  the  rotor  parts 
will  be  preserved.  The  possibility  of  extrapolation  and  its  suf¬ 
ficient  accuracy  are  also  confirmed  by  the  linear  nature  of  the  de¬ 
pendences  t^(T|)  and  tfl(T^). 

Extrapolation  of  t^  and  beyond  =  1150°K  allows  us  to 
more  accurately  establish  the  connection  between  rotor  strength 
characteristics  and  the  temperature  field  and  find  the  most  likely 
places  which  limit  the  maximum  temperature  T^. 

The  linear  nature  of  the  graphs  of  the  dependence  of  blade  and 
disk  temperature  on  gas  temperature  at  3  turbine  inlet,  and  their 
varying  slopes,  confirm  the  results  of  a  theoretical  study  of  tem¬ 
perature  fields  in  gas-turbine  rotors,  presented  in  [1]  and  others. 

Using  the  lines  in  Fig.  6  we  can  find  the  radial  temperature 
field  in  the  disk  and  blades  for  each  stage.  This  was  done  for 
temperatures  T|  =  1100  and  1200°K  (Fig.  7).  Comparison  of  the  lines 
makes  it  possible  to  estimate  the  effect  of  cooling  of  each  rotor 
stage.  The  very  high  temperature  gradient  from  the  middle  of  the 
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Pig.  7.  Radial  temperature  distri¬ 
bution  in'  a  turbine  rotor  in  the 
etationary  mode,  from  thermometry 
results.  Tc  the  left  are  the 
measurement  points  as  In  Pig.  1. 


blade  to  the  base  and  tip  for  stage  I 
attests  to  the  intensity  of  blade 
cooling  in  the  peripheral  and  base 
sections.  This  result  agrees  with 
the  distribution  of  gas  temperature 
along  the  nozzle  vane  cascades  of 
stages  I  and  II  (see  Fig.  5)*  Such 
an  effect  of  cooling  is  not  observed 
for  the  rotor  blades  of  stages  II  and 
III.  The  more  uniform  heating  of  the 
blades  of  stage  II  and  the  high  tem¬ 
perature  differentials  along  the  disk 
indicate  that  the  main  portion  of  air 
fed  to  the  rotor  from  the  compressor 
is  used  to  cool  stage  I  and  the  rear 
support  bearing.  This  can  be  con¬ 
sidered  as  expedient,  since  stage  I 
operates  under  the  severest  tempera¬ 
ture  conditions.  Calculations  also 
show  that  the  temperature  different¬ 
ials  along  the  stage  II  disk  under 
operating  conditions  are  quite  ac¬ 
ceptable,  since  the  reserve  strength 
in  this  case  does  not  drop  below 
standard. 


Judging  from  the  lines  in  Fig. 

7,  from  the  standpoint  of  cooling,  stages  II  and  III  are  under  prac¬ 
tically  identical  conditions,  if  we  consider  that  somewhat  more  heat 
enters  the  st$ge  II  than  the  stage  III  disk  from  the  gas. 


The  change  in  temperature  along  the  leading  and  trailing  sides 
of  the  disks  of  all  three  stages  also  confirms  the  weak  cooling  of 
stage  II.  Nevertheless,  this  cooling  does  exist,  as  attested  to  by 
the  higher  temperature  of  the  rear  side  of  the  stage  III  disk. 
Cooling  of  the  fronts  and  backs  of  the  stage  II  disk  is  symmetrical. 
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since  there  is  practically  no  temperature  differential  here. 

Of  great  interest  is  the  break  in  the  temperature  dependences 
at  the  junction  of  the  blade  and  disk;  this  can  be  explained  by  the 
entrainment  of  part  of  the  heat  by  the  cooling  air  in  the  steady 
state  and,  in  addition,  by  the  thermal  resistance  of  the  blade  fir- 
tree  joints.  Naturally,  thermal  resistance  of  the  joint  is  charac¬ 
terized  by  a  temperature  differential  from  the  stem  of  the  blade  to 


Fig.  8.  Dependence  of  temperature  dif¬ 
ferential  in  a  firtree  Joint  on  time 
when  atartlng-up  a  turbine. 


the  projections  on  the  disk  (Fig.  8), 
but  the  combined  effect  of  both  fac¬ 
tors  is  more  complex.  Figure  9  shows 
the  temperature  dependences  t^(r)  and 
t^(x)  under  nonsteady  conditions, 
determined  theoretically  for  stage  I 
of  a  turbine.  Figure  9a  also  gives 
an  experimental  dependence  (curve  A). 
Comparison  of  Fig.  9a  (with  cooling  of 
the  lower  part  of  the  blade)  and  Fig. 
9b  (without  cooling)  shows  the  great 
part  played  by  the  cooling  channels. 

We  also  see  the  Influence  of  the 
thermal  resistance  of  the  joint  on 
the  temperature  function  with  trans¬ 
ition  from  the  blade  to  the  disk 


3ince,  despite  the  lack  of  cooling  (see  Fig.  9b),  there  is  still  a 
break  in  the  lines.  We  can  conclude  that  thermal,  resistance  of  a 


firtree  joint  with  no  air  passing  through  its  channels  leads  to  a 
10-20°C  temperature  differential  on  the  blade/disk  boundary. 


The  strength  of  the  rotor  disks  is  determined,  to  a  considerable 
extent,  by  the  temperature  differential  from  the  center  to  the  peri¬ 
phery  of  the  disk.  Table  1  gives  temperatures  of  the  centers  and 
peripheries  of  the  disks  of  all  three  stages  of  a  turbine  rotor  under 
steady  conditions  .and  with  T|  *  1100°K,  and  the  corresponding  tem¬ 
perature  differentials.  As  is  to  be  expected,  the  maximum  differ¬ 
entials  are  observed  in  the  stage  II  disk  .  When  designing  a  disk 


'Km — TTrTTr'r 

i  I  i  • 

r  •»  _ _ I  » .  _i_ 


1 1  If  ■"/*  l  l 

a*  !  V  Y\  I  \* 


•  'I  j  I  j 

SOD  ^  ■*"  — ~ 

fr//5l 

5:2  ■ 

•i  my  vi: 

4*c  JLi/  ?y* «'  /]  j  _ 

fT/Vi"  j  j  1 

'  _ 1 


■  optn-  root  !  X)  t  -Xj> 
ln*  disk  I  blad* 


II  J 


»k:M- 


TMWiKu 

-:..!  yt.i /-i/li . J 

1  lyf  l/i  // 


5WK* 


WT77T 

l/i  ■/  i/l  ' 
W  /  /  : 


1  !  g  y 


'  .  ytn«  roct  :  1  x«  1 
in*  M  '> ]bUd« 


Pig.  9*  Radial  temperature  distribution  In  a  turbine  rotor 
at  various  moments  of  time  after  start-up  (shown  on  the 
.  curves  In  minutes  and  seconds):  a  -  blowing  of  cooling  air 
through  channels  in  the  lower  parts  of  root-shrouded 
blades;  b  -  without  air  blowing. 
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for  “strength. ;we> used  temperature  differentials  higher  than  thoafe' 
measured;,  ,  Ae  -at  result  *  for  a  stage  II  disk  operating  under  more* 
severe  temperature  conditions,  the  reserve  strength  with  given  a 
was  3 •'8.  for  the  center  of  the  disk  aud  3*^  for  the  periphery. 


Table  2. 


Position  of 
thermometry 
points  on 
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Table  2  gives  the  blade  temperatures  for  the  stationary  regime 
and  T*  *  1100°K.  Because  of  high  heating,  the  middle  parts  of 
stage  I  disks  have  the  least  strength.  Figure  10  shows  dependences 

of  the  coefficient  of 
reserve  strength  on  T^, 
which  show  that  a  rise 
in  gas  temperature  ahead 
of  the  turbine  is  limited 
to  the  root  section  of 
the  stage  III  rotor  blade. 
This  can  be  explained  by 
the  fact  that  the  material 
of  the  stage  III  blades  is  less  heat-resistant  than  the  stage  I  and 
II  blades.  However,  the  stage  III  vane  cascade  is  under  rather 
stable  temperature  conditions,  and  therefore  the  reserve  strength 
of  the  blades,  even  at  maximum  gas  temperatures,  completely  assures 
reliable  efficiency  of  the  turbine. 

Practice  has  shown  that  the  most  important  parts  of  vane  cas¬ 
cades  Is  the  midsection  of  the  stage  I  rotor  blades,  which  undergo 
the  action  of  the  nonuniform  temperature  field  of  the  gas  flow,  both 
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pig.  10.  Dependence  of  strength  reserve  on  tempera¬ 
ture  T*  fox'  some  of  the  most  stressed  cross  sections 
of  turbine  rotor  blades.  The  points  indicate  the 
calculation  results. 
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along  as  well  as  around  the  turbine  blading.  Therefore,  here  the 
possibility  of  short-term  heating  of  the  midsection  of  the  blade  is 
very  likely  and,  during  turbine  operation  with  maximum  values  of  Tjj, 
the  reserve  strength  of  stage  I  blades  in  the  midsection  will  vary 
and  can  assume  lower  values  than  in  stage  III.  In  addition,  we  cannot 
help  but  consider  the  possibility  of  overheating  of  the  thermally- 
stressed  leading  and  trailing  edges  of  the  blade  foil. 

Thus,  we  can  consider  that  the  rise  in  gas  temperature  ahead 
of  the  turbine  is  limited  by  the  reserve  strength  in  the  midsection 
of  the  stage  I  rotor  blades.  The  influence  of  variations  in  the  tem¬ 
perature  field  around  the  blading  on  the  temperature  and,  conse¬ 
quently,  on  the  strength  characteristics  of  the  blades  requires 
additional  study. 

When  processing  the  thermometry  results  we  also  examined  the 
nonstationary  temperature  state  in  the  process  of  heating  of  the 
rotor  during  start-up.  Some  of  the  more  typical  of  the  obtained 
graphs  for  the  temperature  dependence  are  given  in  Pigs.  11-13  for 
stages  I,  II,  and  III  (the  marking  of  the  measurement  points  is  as 
In  Pig.  3). 

The  change  In  temperature  of  a  stage  I  disk  and  blades  is 
shown  in  Fig.  11a.  The  first  start-up  (see  Pig.  11a)  was  not  stable, 
and  was  characterized  by  sharp  variations  in  t|  right  up  to  the 
moment  of  transfer  from  start-up  fuel  to  main  fuel  at  the  13th 
second.  Such  a  turbine  operating  regime  allowed  us  to  explain  that 
the  thermal  inertia  of  the  blades  is  very  low:  variations  in  blade 
temperature  follow  the  variations  in  the  gas  temperature.  This  in¬ 
dicates  the  significant  predominance  of  the  amount  of  heat  from  the 
gas  over  the  amount  of  heat  fed  from  the  working  part  of  the  blades 
to  the  disk  and  the  correspondingly  slight  difference  between  the 
temperature  of  the  blade  midsection  and  that  of  the  gas. 

Above  it  was  noted  that  the  rear  of  the  stage  I  disk  is  hotter 
than  the  front.  The  trend  of  the  temperature  at  points  1,  3,  5,  2, 


4,  and  6  of  the  disk  in 
the  nonstationary  regime 
makes  it  possible  to  . 
trace  how,  from  the  very 
first  seconds,  there  is 
a  discrepancy  in  the 
thermocouple  readings; 
the  rear  of  the  disk  is 
more  strongly  heated. 
This  indicates  the  main¬ 
tenance  of  the  distri¬ 
bution,  in  the  turbine, 
of  the  cooling  air  by 
stages,  in  both  steady 

The  maximum  rate  of  heating  of  the  peripheral  part  of  the  disk 
in  the  first  4  seconds  in  60-90°C/s. 

Let  us  mention  the  feature  of  function  t|(x)  noted  on  the 
graphs  for  ail  stages.  If  the  nature  of  the  change  of  t^  with  time 
can  be  expressed  approximately  by  an  exponential  function,  tempera¬ 
ture  t|  changes,  in  the  first  4-5  second,  by  a  parabolic  law,  and 
then  the  rate  of  its  increase  decreases  and  it  asymptotically  ap¬ 
proaches  the  rated  value.  These  regularities  must  be  taken  into 
consideration  when  stating  problems  of  heat  conduction  for  a  gas- 
turbine  rotor. 

The  graphs  given  in  Pig.  11  were  used  to  determine  the  depen¬ 
dences  of  the  temperature  differentials  along  the  disk  on  the  time 
for  all  stages  (Pig.  14).  The  graphs  3how  that  during  start-up  the 
temperature  differentials  along  the  radius  of  the  stage  I  disk  are 
somewhat  higher  than  in  the  stationary  temperature  state.  If  in 
steady-state  regimes  t^  -  t^  is  60°,  during  start-up  this  difference 
reaches  95°,  while,  the  differential  tg  -  tg  reaches  140°  instead  of 
87°  in  the  stationary  state.  We  can  observe  anomalies  in  which  the 
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Fig.  13.  Change  in  disk  and  blade  teaparature  of  stage 
IZZ  during  start-up. 


and  transitional  regimes. 


Fig.  14.  Change  In  temperature  differential  from  the 
periphery  of  the  disk  to  Its  center  during  turbine 
start-up. 
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differential  increases  continually  with  certain  variations  (line 
fc5  “  ^i*)»  but  they  cannot  be  considered  typical.  The  nature  of  the 
curve  tjj  -  t2  indicates  that  because  6f  the  more  rapid  heating  of 
the  disk  periphery  on  the  front  side  (lines  5  in  Fig.  11)  the  differ¬ 
ential  t^  -  t2  is  first  greater  than  zero  and  only  after  the  15th 
second  does  it  become  negative,  i.e.,  t^  <  tg.. 

The  trend-  of  dependences  tg  -  t^  and  t^  -  t2  shows  the  absence, 
during  start-up,  pf  a  maximum  temperature  differential  for  alternate 
points . 

The  maximum  temperature  differential  for  a  stage  I  disk  durin* 
start-up  is  lower  than  that  which  served  as  the  initial  one  for  the 
strength  calculations  with  steady  heating.  Consequently,  the  trans¬ 
itional  temperature  state  has  no  influence  whatsoever  on  the  effi¬ 
ciency  and  reliability  of  a  disk.  The  time-dependence  of  the  tem¬ 
perature  of  stage  II  blades  and  sisks  3hows  the  low  thermal  Inertia 
of  the  blades  (curve  18  in  Fig.  12).  Thanks  to  stable  start-up  we 
easily  see  that  the  rate  of  heating  of  the  central  part  of  the  disk 
is  somewhat  less  than  that  of  the  peripheral  part. 

Curves  11,  13,  12,  and  14  show  that  in  the  first  seconds  after 
start-up  the  center  of  the  stage  II  disk  is  heated  more  rapidly 
than  its  middle  part.  To  explain  this  we  must  simultaneously  analyze 
the  temperature  dependences  for  stage  I  and  II  disks.  Above  we  noted 
that  air  from  the  compressor  is  almost  completely  expended  on  cooling 
only  the  leading  edge  of  the  stage  I  disk.  Therefore,  we  cannot  ex¬ 
plain  the  higher  heating  at  points  11  and  12  than  at  points  13  and 
14  (see  Fig.  12)  in  the  first  seconds  after  start-up  by  the  fact 
that  the  center  of  the  stage  II  disk  is  washed  by  hot  air  from  the 
compressor  during  this  time,  all  the  more  so  since  in  this  case, 
point  1  of  the  stage  I  disk  (see  Fig.  lib)  would  be  heated  most 
rapidly  of  all  in  this  case.  However,  its  temperature  is  lower 
than  at  point  2,  both  during  start-up  as  well  as  during  subsequent 
operation  of  the  turbine. 
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Such  a  heating  peculiarity  can  be  explained  as  follows.  Cooling 
air  on  its  way  to  exhausting  into  the  turbine  blading  between  the 
first  rotor  and  the  second  nozzle  cascade  loses  so  much  head  that  it 
cannot  counteract  the  pressure  on  the  part  of  the  gas  flow.  At  the 
same  time,  the  difference  in  pressures  at  the  inlet  to  and  outlet 
from  the  root  section  of  the  second  and  third  nozzle  cascades  (Pig. 

1)  leads  to  suction  of  the  gas  into  the  gap  between  the  first  rotor 
and  the  second  nozzle  cascade.  ■  The  temperature  of  the  gas  entering 
the  gap  is  low,  since  in  the  lower  part  of  the  first  nozzle  cascade 
th-re  is  mixed,  with  the  gas  flow,  cooling  air  from  the  gap  between 
the  first  nozzle  cascade  and  the  first  rotor.  When  the  gas  enters 
the  gap, its  temperature  is  reduced  still  more  due  to  mixing  of  the 
air  coming  from  the  cooling  channels  in  the  lower  part  of  the  first 
rotor  cascade.  The  gas,  mixing  with  the  cooling  air,  penetrates 
through  the  radial  openings  in  the  skirt  of  the  stage  I  disk  and 
washes  its  rear  side.  This  is  attested  to  by  the  synchronous  vari¬ 
ation  of  temperature  dependences  2  and  4.  A  small  amount  of  the  gas- 
air  mixture  leaks  through  the  labyrinth  seal  between  stages  I  and 
II,  and  also  through  the  radial  channels  in  the  skirt  of  the  stage 
II  disk,  and  enters  the  blading  of  the  turbine  between  the  second 
nozsle  cascade  and  the  second  rotor.  The  fact  that  this  amount  of 
gas-air  mixture  is  really  slight  is  confirmed  by  the  lower  tempera¬ 
ture  at  point  13  than  at  point  11  during  start-up  (see  Pig.  12). 

The  gas-air  mixture,  reaching  the  central  recess  between  the 
stage  I  and  II  disks,  passes  to  stage  III,  washing  the  walls  of  the 

opening  of  the  second  disk  and  heating  its  central  part.  Along  the 

way  it  is  again  mixed  with  a  certain  quantity  of  cooling  air  from 
the  central  opening  of  the  first  disk,  and  its  temperature  continues 
to  drop.  Then  from  the  cavity  between  the  second  and  third  disks 
the  gas-air  mixture  exits  to  the  blading  through  the  radial  channels 
in  the  skirt  of  the  stage  II  disk.  Here  the  rear  surface  of  this 
disk  and  the  central  part  of  the  third  disk  are  washed,  which  is  well 
illustrated  by  Pig.  13*  The  temperature  at  point  21  first  rises 
rather  rapidly  but,  with  influx  of  heat  from  the  peripheral  part, 

heating  on  the  part  of  the  gas-air  mixture  plays  a  lesser  role,  and 
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t^;/temperatjare- at  points?  22a  23^.  24  becomes  higher  tbaii  at' 

pqint,.2JL«  .  *.  .  ' 

,  It  is  important  to  note  that  circulation  of  the  ga's-air  mixture, 
in  final  analysis,,  «cools  stages  II  and  III’  of  the  turbine  rotor. 

The  temperature  of  the  mixture  trashing  the  rear  surface  of  the  first 
disk  does  not  exceed  390-400°C,  while  that  washing  the  rear  surface 
of  the  second  disk  and  the  front  surface  of  the  third  disk  does  not 
exceed  310-320°C.  „ therefore,  the  gas-air  mixture  heats  the  disks 
only  during  the  first  few  seconds  after  start-up.  As  the  disks  be-  - 
come  heated  from  the  periphery,  the  mixture  acts  as  a  coolant, 
assuring,  under  stationary  turbine  operating  conditions,  a  stable 
temperature  field  for  the  disks,  which  is  t^'en  into  consideration 
in  strength  calculations  of  disks.  In  thv  aenoe  of  such  circula¬ 
tion,  heating  of  the  central  parts  of  stage  II  arid  III  disks  would 
be  determined  by  the  temperature  of  the  air  coming  from  the  compres¬ 
sor,  which  would  lead  to  an  increased  temperature  differential  along 
the  disk  radius  and,  consequently,  to  increased  thermal  stresses. 

Heating  of  the  central  part  of  rotor  disks  also  aids  in  rapid 
heating  of  the  rotor  and  passage  to  the  stationary  temperature 
regime . 

Judging  from  the  graphs  of  nonstationary  heating  (see  Pig.  11), 
the  washing  medium  begins  to  cool  the  stage  I  disk  approximately 
20  seconds  after  start-up,  and  the  disks  of  stages  II  and  III  -  after 
30  seconds . 

The  rule  for  the  change  in  temperature  differential  in  the  stage 
II  disk  from  the  periphery  to  the  center  during  start-up  is  the  same 
as  for  the  stage  I  disk,  but  the  magnitude  of  the  differentials  in¬ 
creases  (see  Pig.  14).  In  connection  with  the  fact  that  the  disk  is 
cooled  symmetrically,  the  maximum  temperature  differential  in  the 
first  seconds  of  turbine  operation  is  also  observed  for  alternate 
points.  For  strength  calculation  of  the  stage  II  disk  we  assigned 
a  temperature  differential  of  300°C,  and  therefore  the  maximum 
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occurring  during  start-up  for  t^g  -  t 12  211(5  ^5  ~  is  not  defini¬ 
tive.  A  negative  differential  can  occur  at  the  very  beginning  of 
scart-up,  as  already  stated,  due  to  heating  of  the  central  part. 


The  rules  for  temperature  distribution  in  stage  III  disks  and 
blades  during  turbine  start-up  are  the  same  as  for  stages  I  and  II 
(see  Pig.  13). 


The  temperature  differentials  in  a  stage  III  disk  are  shown  in 
Fig.  14.  From  the  graphs  we  3ee  that  heating  of  point  21  in  the 
first  few  seconds  after  start-up  all'ws  us  to  make  the  temperature 
increase  monotonic.  As  a  result,  function  t -  t21  became  almost 
monotonic,  just  as  c2g  -  t2-,  • 
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THE  NONaTATiOHAftY  HEAT  CONDUCTION, OF  BODIES  OF  COMPLEX  CONFIGURATION 


M,  N.  Galkin  and  S.  P.  Kolesnikov 


Many  parte  of  airborne  vehicles  and  their  engines  have  sharp 
edges  and  projections  in  the  form  of  dihedral  and  polyhedral  angles. 
During  operation j  these,  parts  are  subjected  to  short-term  heating  or 
cooling,  and  from  a  thermal  standpoint  can  be  examined,  as  semi- 
bounded  bodies. 

In  Ghis  work  we  give  our  general  theoretical  solutions  to 
problems  on  determining  the  temperature  field  and  the  heat-retaining 
capabilities  of  semi-bounded  bodies  having  the  form  of  dihedral  and 
polyhedral,  acute  and  obtuse  angles.  The  solutions  are  approximate 
and  have  been  obtained  for  the  case  of  constant  temperature  on  the 
surface  of  the  bodies. 

The  essence  of  the  solution  of  the  multidimensional  nonstation¬ 
ary  problem  is  first  examined  in  detail  u§ing,  as  our  example,  a 
dihedral  angle  (the  two-dimensional  problem),  and  is  then  extended 
to  polyhedral  angles. 


The  particular  solution  to  the  problem  for  the  temperature  field 
of  a  dihedral  right  angle  (Fig.  1)  is  the  product  of  the  solution 
for  two  seml-btpunded  bodies  (A.  V.  Lykov:  Teoriya  teplcprovodno8ti 
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s  .|  (The  theory  of  heat  conduction) ,  GITTL,  1952): 

y  ^  T"  » =  1  —  erf  / — erf  ,  . 

*  \  •  \2ph»t  J  |/«t  / 

\  j  or  ft=l  —  erf/ — 7=r-)er(/ — ^z=~\  <  (1) 

V/  .  \  2J  Fo.,  /  V  2  I  Fo.r  / 

y|  where  8=  ^x'  p'  is  the  relative  temperature  of 

PI*.  1.  Otter®-  nln*  the  ^  ^ 

d*"Sdr«iu^su%1-°9o«1?*  the  body  at  a  point  wi.th  coordinates  x  and  y  at 

time  t;  tQ  is  the  constant  temperature  of  the 
body  surface;  tQ  is  the  initial  temperature  of  the  body;  a  =  X/cp 
is  the  coefficient  of  thermal  conductivity  of  the  body  material; 

X,  c,  and  p  are  the  heat  conduction,  specific  heat,  and  density  of 
the  body  material;  Fox  =  aT/x  is  the  Fourier  criterion;  y  =  y/x; 
erf(U)  is  the  Gauss  error  function,  which  is  equal  to  zero  when  U  = 

=  0  and  is  practically  equal  to  one  when  U  >_  2. 


PI*.  1.  Deter®'  nln*  th* 
temperature  fit  -1  of  th* 
dihedral  angle  f  •  90°. 


To  derive  the  approximate  general  solution  we  used  the  scheme 
shown  in  Fig.  2.  The  temperature  at- point 
x,y  of  investigated  angle  <p  was  determined  yx 

on  the  basis  of  the  equivalent  right  angle,  /  \ 

at  which  the  examined  point  has  coordinates  \ 

Xy.  Considering  this,  equation  (1)  can  be  yyfesyj - jf  ■ 


represented  in  the  form  rr - t — - - — 

J  Reproduced  from 
I  best  available  copy, 

M, 

V  2Vrflt  /  \2/ut/ 

or 

fi  —  l  —  erf  (  erf  f — • 

\2VFox)  \2)'fo,/ 


Pig.  2.  Determining  the 
temperature  field  of  di¬ 
hedral  angle  e  of  any 
magnitude. 


Equation  (2)  allows  us  to  approximately  calculate  the  tempera¬ 
ture  field  or  a  body  with  any  angle  <p,  while  it  gives  a  precise  so¬ 
lution  in  the  two  particular  cases  when  angles  <p  are  180°  and  90°. 
If  <p  =  180" ,  then  tan  90°  =  00 ,  erf(°°)  =  1  and 


Equation  (3)  is  the  known  precise  solution  of  the  problem 
a  semi-bounded  body. 


If  9.  *  90°,  then  tan  45°  *  1  and  equation  (2)  is  simplified,- 
assuming  the  fora  of  equation  (1),  which  is  fthe  precise  solution  . to 
the  problem  for  a  right  angle.' 

When  calculating  bodies  with  small  angles  9  and,  in  particular,' 
when  angle  9  ' 0,  i.e.,  the  sides  of  the  body  are  parallel,  an 
indeterminacy  arises  in  equation  (2).,  since  tan  0*0,  erf(0)  *  0, 
•mid  *  *  1.  To  expose  this  indeterminacy  in  equation  (2)  we  must  ex¬ 
clude  the  coordinate  x: 

rw7.  ' 


-  l—erf  /  — —Wf/-2 — 

\*v»y  \2  y'ax) 


where  Y  is  half  the  distance  between  sides  (or  half  the  plate  thick¬ 
ness). 

To  estimate  the  accuracy  of  approximate  solution  (<;)  we  experi¬ 
mentally  studied  the  temperature  field  of  bodies  haying  various 
thermophysical  properties  and  angles  9  from  0°  to  180°.  The  maximum 
deviation  of  the  calculation  data  from  experimental  data  occurred 
on  the  bisector  of  angle  9  and  in  the  range  9  *  50°-l80°,  and  did 
not  exceed  ±5%.  The  deviations  were  estimated  using  the  criterion 
(fcnp  “  ”  to^;*  where  fcnp  “  is  the  excess  temperature, 

calculated  from  formula  (2),  and  t  -  tQ  is  that  found  from  experi¬ 
ment  (tnp  is  the  approximate  temperature,  t  is  the  real  temperature) 


;-:T 


Figure  3  shows  the  values  of  these 
criteria  for  angles  9  from  0°  to  180° 
(curve  1). 


[ _  N/  |  From  Fig.  3  it  follows  that  theo- 

o  ij  «3  j3  *3  30  es  73  so  retical  solution  (2)  gives  an  accurate 

pi*.  3.  Dependence  or  the  criterion  result  for  angles  9  *  180°  and-90°.  an 

(tnp  -  *  V  on  ***•  *»  ehlu>  * 

•ottriiin*  th»  maxiDua  devincioru  of.  overvalued  result  for  angles  180°  >  9  * 

the  calculated  teaperature  field  froa 

^uSSd^S^uuiin  I°-  >  90°  (the  maximum  error  is  obtained  at 

fro*  aquation  <5>.  9  *  140°  and  reaches  5%),  and  an  under¬ 

valued  result  for  90°  >  9j>  0°  (maximum  error  is  obtained  when  9*0° 
and  reaches  .7%).  - 


Pi*.  3.  Dependence  or  tha  criterion 

<»np  *  V/(t*  *  V  on  •‘W1*  eheu> 

aotarltln*  tha  aaxloua  deviation*  of. 
the  calculated  teaperature  field  froa 
the  experlaantal  one.  Curve  1  ie  con- 
atructed  froa  equation  (2),  curve  2  - 
froa  equation  (3). 
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These  deviations  occur  since,  according  to  the  scheme'  (Fig.  2) 
used  for  deriving  equation  (2),  we  get  an  undervalued  heat  flow  along 
the  bisector  of  the  angle  when  180°  >.  <p  >  90°  and  an  overvalued  flow 
when  0°  <  9  <  90° .  • 


To  increase  the  accuracy  in  calculating  the  temperature  field 
of  bodies  with  small  angles  <p ,,  it  is  expedient  to  seek  the  theoreti¬ 
cal  solution  on  the  basis  of  the  scheme  shown 
in  Fig.  4,  which  generally  excludes  heat  trans¬ 
port  along  the  bisector  of  the  angle.  In  such 
a  case,  each  cross  section  x  of  angle  <p  at 


Mg.  *.  thriving  th«  »pproxi-  depth  Y  can,  from  a  thermal  standpoint,  be  ex- 

m&tt  colutlon  to  th«  problem 
for  the  temperature  field  of 
bod'.et  wltn  smell  angles  (f  < 

*  W 5"  which  the  solution  (see  Lykov’s  book  mentioned 


amined  as  half  the  thickness  of  a  plate,  for 


mi 


■Mi 


above)  has  the -following  form: 


*“1  “i  *{-£§72}'  (5) 

n  *»1 

where  .4,  —  —  ( —  l.'-*  is  the  initial  thermal  amplitude  I  —  1) . 


In  the  particular  case  when  <p 
elude  coordinate  x: 


0,  in  equation  (5)  we  must  ex- 


Reproduced  from 
best  available  copy. 


(6) 


Equation  (6)  is  the  precise  solution  of  the  problem  for  an  un¬ 
bounded  plate.  Curve  2  in  Fig.  3  corresponds  to  the  relative  ex¬ 
cess  temperature  — — --  calculated  from  formula  (5)-  Figure  3  shows 

ti—  '9 

that  for  cp  =  0°-40°  we  must  use  approximate  solution  (5),  while  for 
<p  —  40°-l80°  we  use  solution  (2).  Here  the  maximum  temperature-cal¬ 
culation  error  does  not  exceed  9%  (on  the  bisector,  and  only  for  a 
body  with  angle  cp  =  40°). 


In  Fig.  3  we  also  see  that 
mate  solutions  (2)  and  (5)  dees 


the  calculation  error  using  approxi- 
not  exceed  2%  for  angles  9  »  0°-20° , 
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9  •  70°-110° ,  and  9  -  l60°-l80° 


0  0,1  c.n  c,$  o,3 

nt.  S;  Cvmi  of  ths  distri¬ 
bution  of  rslstlrs  sxosm  tta- 
psraturs  in  dihsdral  anglss 
SO*  <  »  <  180*  u  a  function 
of  tEs  orltsrlon  Po/tan**/2. 


0,2  Ofl  CJt  0,6  Y/ Igj 

tit .  6.  Curva*  of  ths  distri¬ 
bution  of  relative  axccas  taa- 
paratura  In  dihsdral  anjlss 
0»  «  a  <  <0*  as  a  function  of 
ths  erl'F’srlon  Foj/tan**/2- 


Figures  5  and  6  show  curves  of  the  distribution  of  relative 
excess  temperature  in  dihedral  angles  for  various  values  of  the 

p 

criterion  Fox/tan  9/2,  constructed  from  solutions  (2)  and  (5),  re¬ 
spectively. 


The  specific  heat  flow  arriving  at  the  body  at  time  t  per  unit 
area  of  heating  surface,  located  at  distance  x  from  the  vertex  of  the 
— 2 _ H  angle  (Fig.  7),  is  determined  from  the 


imiiiiimiiLiiMi  i  nv 


Kiv't 


Fourier  law  of  heat  conduction: 

I  OH  J.v-r' 

where  is  the  temperature  gradient 


\  on  surface  ox  of  the  dihedral  angle  at  time 

t,  determined  with  consideration  of  formula 
(2)  or  (5)  as  a  function  of  the  value  of 

Pl«.  7.  Dlxirsa  of  ths  host  flow 

on  ths  surf  set  of  a  dlhsdral  dngJ.e  9. 


It  should  be  noted  that  formulas  (2)  and  (5)  give,  for  all  prac¬ 
tical  purposes,  an  accurate  temperature  at  the  body  surface.  There¬ 
fore,  the  value  of  the  heat  flow  found  with  consideration  of  formulas 


■  i;- 


(2)  or  (5)  will  fee,,  f dr  all  .practical  purposes,  accurate. 
Let  us  find  the  heat  flow  when  9  «  lSO0-^0. 


Having  determined  the  gradient  [3t(x ,y ,x)/3y]y=0  with  con¬ 
sideration  of  formula  (2)  for  the  body  surface  (y  *  0),  we  find 


<*~7 irl''-'Jcd(iSr)- 


where  b  *  /XTcp  is  the  coefficient  of  accumulation  of  heat  by  the 
body  material. 


The  quantity  of  heat  passing  in  time  interval  t  across  a  unit 
of  surface  located  at  distance  x  from  the  vertex  of  the  angle 


Substituting  qx  from  formula  (7),  after  integration  we  get 


ferff  x'*?'2) 

* ’  *8  -  2  pj 

r-  *J|gV2M 

IUI'«  ' 

j'.Tsr^  1 

J  l  Ul'FOx  '  2 1  .T  |  l-o.r  V  4Fcy 


The  quantity  of  heat  entering  the  body  on  segment  X  of  unit 

width  in  finite  time  x  _ 

.v  I  Reproduced  from 


.v  !  Reproduced  from 

Qx  \  Qx  if.x,  {  best  available  copy. 


while  after  substitution  of  Qx  from  formula  (8)  and  integrating 


0-v~-iTft— ,"u''/f[e"Tnr+ 


’ifivsr 


I—/'  ,fl' 
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,  2M 

\  —  C 
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p 

where  Fo^  *  a t/X  Is  the  Fourier  criterion. 

Let  us  find  the  heat  flow  for  angle  q>  *.  40°-0°. 

Having  determined  the  gradient  [3t(x,y,x)/3y]yjlt0  with  con¬ 
sideration  of  formula  (5)  for  the  body  surface  (y  *  0)  we  get 

t^-v—rS  oxp(-|‘J^S3-)-  <12> 

«»t 

After  substituting  the  value  of  qx  from  formula  (12)  into  form¬ 
ula  (7*)  and  Integrating,  we  get 

UJ 

here  the  average  relative  excess  temperature  on  the  surface  of  a 
dihedral  angle  in  section  x 

S\»A=%!U,V  4  „p  CU) 

t<-tj  ii  t*i  V  "i fPsfll  ' 

#1—1 

The  quantity  of  heat  entering  the  body  in  segment  X  of  unit  width 
in  finite  time  x  _ —rr~ - 

_  .  I  ReDroduced  Horn 

Q,~  (  Q,dx.  ^ 

6 

Considering  formulas  (13)  and  (14),  after  integration  we  get 


-  2Fo,.  'pi  /  ,  F«v  \ 

A-  V  Ei  (  — -A-  )  ta  y  .» 

lg3?/2  U  \  "  Ig-  =  2  • 

e-t,)X\‘x~ - 5^= - - 

2  V  r"x 


Prom  equation  (8)  it  follows  that  when  x  -*■  <*>,  heat  flow  Q? 
reaches  its  maximum  value 


Q*« (16) 

Equation  (16)  is  the  known  particular  solution  to  the  problem  of 
the  heat-retaining  capability  of  a  semi-bounded  body  (see  Lykov). 
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The  ratio  Qx/Q*  -  the  actual  amount  of  heat  to  the  maximum 
amount  -  is  the  parametric  criterion  Kx* Which  characterizes  the  in¬ 
fluence  of  the  angle  (bilateral  heat  feed)  on  the  heat-retaining’ 
capabilities  of  the  body  in  section  x. 

From  equations  (9)  and  (16)  it  follows  that  for  angles  40°  < 

<  <P  <  180°  the  criterion 
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*  V  2 1  r..,  /  2|  .t|  Fn,  \  4^.1 


while  for  angles  0  <  9  <  lJ0°,  from  equations  (13)  and  (16)  we  get 
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Figure  8  shows  the  dependence  of  criterion  K  on  criterion 

2  * 
l/Fo  *  x  /at,  calculated  from  formula  (18),  for  various  values  of 

angle  9/2  within  limits  from  5°  to  20°,  while  Figs.  9  and  10  show 

analogous  dependences,  calculated  from  formula  (17)  for  angles  <p/2 

=  20°-1»5o  and  (p/2  *  45°-75°,  respectively. 


From  the  graphs  we  see  that  the  influence  of  the  angle  on  the 
heat-retaining  properties  of  a  body  increases  with  decreasing  cri¬ 
terion  l/Fox,  i.e.,  as  the  examined  section  x  approaches  the  vertex 
of  the  angle,  while  with  decreasing  angle  <p  this  influence  increases, 


From  equations  (17)  and  (18)  we  see  that  if  criterion  l/Fo 

A 

vanishes,  for  any  value  of  angle  9/2  (except  9/2  *  90°,  when 
tan  9/2  *  °°)  the  criterion  K  *  0.  Consequently,  heat  flow  Q  = 

a  A 

=  Q*KX,  passing  through  the  vertex  of  any  angle,  is  equal  to  zero. 

If  criterion  l/Fox  ”,  in  accordance  with  equations  (17)  and 
(18)  Kx  ■*  1.  In  this  case  the  heat  flow  reaches  the  value  Q*  which 
characterizes  the  heat-retaining  capability  of  a  semi-bounded  body 
(16). 


Theoretical  analysis  shows  that  section  xQ,  behind  which  the 
angle  is  heated  practically  as  a  semi-bounded  body,  can  be  found 
from  the  condition 


1  Y'J  "•■  ;W 

...•/if  ■,  1 


■-  --■' *t  -Vffy  ^  „  V tT^:«5r  t$j 

I 


«•  ;-j'f?-— --~^2.  *  ~  -.  :  ‘  (19) 

Actuslly,  when  the  argument  r==— 2.,  function  \erf  z  is  practi¬ 
cally  equal  to  one,  while-  function  -Ei(.-z2) . *  0.  V-' 
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In  this  case,  from  equation  (if’)'  it  follows  that  K  *  1. 

An  analogous  conclusion  is  also  obtained  from  equation  (18). 

For  example,  as  Lykov  showed,  when  z  >  2  formula  (14)  for  the  average 
temperature  can  be  given  in  the  form 

c-i-  ‘  sm.. 

■  V«  lg?/2 

Considering  this,  from  formula-  (18)  it  follows  that  K  =  1. 


Using  criterion  K  ,  formulas  (9)  and  (13)  can  be  written  in  a 

A 

general  form,  convenient  for  practical  calculations: 

C?,=^i(/r-/0)VT/C.r.  (20) 

In  many  practical  cases  it  is  important  to  know  the  total  quan¬ 
tity  of  heat  Qy  absorbed  by  the  surface  of  a  body  in  a  certc  n  sec- 
tion  X  (see  Fig.  7).  This  heat  can  be  expressed  in  percents  Kx  of 
the  heat  flow  passing  in  this  same  section  X  to  a  semi-bounded 
body,  in  accordance  with  equation  (16): 

Q'x=-2L(tc-(,),\ Yr.  (21) 

The  value  Kx  =  QX/QJ  is  a  parametric  criterion  which  character¬ 
izes  the  influence  of  the  angle  (bilateral  heat  feed)  on  the  heat- 
retaining  properties  of  a  body  In  section  X. 

* 

From  equations  (11)  and  (21)  it  follows  that  for  angles  iJ0°  .5 
<  <9  <  180° 

/Cv=crJ  (  — VFo-v  ( - 1) - -!£iJL=-  Hi  {—  (22) 

A  \2]'Fov  /  t"=-2)’a  \  /  ■!)  -t|'F«y  \  4  Fo  v  j 
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while  from  equations  (15)  and  (21),  ,for  angles  0a  <  <p  <  ,40° '  we  geft 
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Figure  11  shows  the  dependence  of  criterion  Kx  on  criterion 
1/Fo^,  calculated  from -formula  (23),  for  various  values  of  the  angle 
<p/2  from  5°  to  20°.  “  I 

Figures  12  and  13  give  analogous  dependences,  calculated  from 
formula  (22)  for  angles  q>/2  *  20°-45°  and  <p/2  *  45°-75°,  respectively. 


Criterion  Kx  allows  us  to  write  formulas.  (11)  and  (15)  in  a  form 
which  is  general  and  convenient  for  practical  use; 

2^  ,  A  A  \  \P  «  .'“‘r/  A  _  S.  v 


<?*=-=  (/,-/<,)*  I  T Kx. 


The  theoretical  dependences  obtained  were  verified  experimentally 
on  a  devict  which  is  shown  schematically  in  Fig.  14.  Into  a  bath  of 
"iPf-09  molten  pure  aluminum,  having  a  temperature 

_ J  |  of  662  ±  2°,  i.e.,  superheated  by  3  ±  2°,  we 

E  immersed  for  60  seconds  prisms  made  of  re- 

- rn  fractory  magnesite  brick  whose  coefficient 


of  thermal  conductivity  a  *  8*10“^  m2/hr. 

The  working  element  of  the  prisms  was  a  di¬ 
hedral  angle  which  varied  within  limits 
<p/2  =  15°-90°  at  15°  intervals.  The  other 
faces  of  the  prisms  were  insulated  by  sheet 
asbestos  (a  =  0.712*10“^  m2/hr)  20  mm  thick. 


XT/!  /  At  the  control  point  of  each  prism  we 

installed  a  Chromel-Alumel  thermocouple 
whose  readings  were  recorded  by  an  EPP-09 
JtuiyiJi  tMMMtm  &iror  electronic  potentiometer.  The  thermocouple 
«uh#dr»*  »n«i«».  leads  (0.5  mm  in  diameter)  were  insulated 

near  the  hot  Junction  by  a  two-channel  porcelain  tube  3  mm  in  diameter 
and  a  cambric  covering. 


A* 

•»&: 

X* 

1 

iv 

J 

i 

Fig.  14.  oiagraa  of  the  device  for 
■tuning  tne  teaperature  field  of 
dihedral  angle*. 
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Based  bn  the  location  of  the  control  point  ,r£i.ative  to  the 
vertex  of  the  angle,  all  the  prisms  were  divided  Into  t^o  groups : 
the  first  group  had  angles  <?s/2  ■  15° ,  30° ,  and  45°  j  t,V,e  secortd  group 
had  angles  qp/2  »  45° ,  60°,  75°  >  arid  -90°  •  Tht>  control^  points^  had  the 
following  coordinates:  group  1  -  x  '»  22  mm,  y  *v-,5..5  nun  (y  *  0.25)  j 
group  2  -  x  ■  6.8  mm,  y  ■  6.8  ram  (y  «  1). 


0,CS  0.10  0,15  0,20  0.2s'  F0X 

50%,  • 

Fig.  15.  SxptrlMntal  and  ctloulatlon 
graph*  of  ttM  chang*  in  relative  oxoos* 
taaparatur*  at  tha  control  point  of 
grout.  1  prl*a*. 


0,5  J  1,5  Z  2,5  rax 
10  "  20  to  ~UQ  50  %>  ■ 

Fig.  1$.  Sun  a*  Fig.  15.  group  2 
prints. 


Figures  15  and  16  show  the  ex¬ 
perimental  points  and  calculation 
curves  of  the  change  in  relative  ex¬ 
cess  temperature  at  the  prism  con¬ 
trol  point  as  a  function  of  criterion 
Fox  and  time  t.  From  the  graphs  we 
see  that  the  calculation  results  are 
in  good  agreement  with  experimental 
data. 

Now  let  us  examine  the  tempera¬ 
ture  ''ield  and  heat  flow  on  the  sur¬ 
face  of  a  trihedral  angle. 


B  Figure  17  shows  the  trihedral 

Fig.  17.  Dlggraa  of  eh*  distribution  _ .  _ _  .  ,  ,  _ ,  _  , 

or  h**t  flow  on  th*  iurf*o*  or  *  tn-  angle  OABC  whose  temperature  field 

h*dr*l  «ngl*. 

can  be  approximately  described  by 

the  following  equation,  derived  from  the  scheme  given  shove  (Fig.  2): 


iSCVi,  ¥g. 
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where  i t/1*  is  the  relative  temperature  of  the  body  at  point 

/f  —  /j 

a  with  coordinates  x,y,z  at  time  t;  9*  and  9'  —e  the  plane  and  di¬ 
hedral  angles  of  the  trihedral  angle;  z  »  z/x. 

The  heat  flow  entering  the  trihedral  angle  in  time  t  across  a 
unit  area  with  coordinates  xy  can  be  found,  approximately,  from  thb 
following  expression: 

Qlxt)  ~y=(tc-ta)  Vx-K^K,,  (26) 

where  K  and  K  are  criteria  which  characterize  the  influence  of 

^  V 

angles  9'  and  9"  on  the  heau-re Gaining  properties  of  a  trihedral 
angle  at  point  xy.  They  are  functions  of,  besides  angles  9*  and  9”, 
criteria  1/Fo  and  1/Fo  ,  and  are  determined  from  the  graphs  given 

x  y 

above  in  Figs.  8-10. 


Thus,  the  formulas  obtained  from  theoretical  studies  make  it 
possible  to  calculate  the  nonstationary  temperature  field  and  the 
heat-retaining  properties  of  bodies  having  the  form  of  di-  and  tri¬ 
hedral  angles. 


The  results  of  an  experimental  verification  of  these  formulas 
confirmed  the  fact  that  they  are  correct,  which  allows  us  to  recom¬ 
mend  them  for  actual  calculations . 


HEAT  TRANSFER  IN  A  FLUID  FLOW  WITH  CONSIDERATION  OF  THE  INFLUENCE 
OF  DISSIPATION 


M„  N.  Galkin  and  S.  V.  Lomazov 


When  turbine  blades  are  cooled  by  convection  heat  transfer, 
the  cooling  fluid  flows  at  high  velocities  through  very  narrow 
channels.  In  this  case  the  distribution  of  temperature  along  the 
flow  is  influenced  by  the  dissipative  heat  released  as  a  result  of 
viscous  friction. 


With  laminar  flow  of  a  Newtonian  fluid,  the  dissipation  of 
mechanical  energy  per  unit  time  [1]  ie  defined  by  the  formula 


t  [’  /  dvi  ^  dv*  y  ii  f 
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(1) 


where  n  is  the  coefficient  of  dynamic  viscosity  and  av^/Sx^  are 
derivatives  of  the  velocity  components  with  respect  to  the  coordi¬ 
nates  . 


In  the  integrand,  according  to  the  symbolic  notation  of  tensor 
algebra,  there  is  summation  over  the  repeating  indices  in  the  pro¬ 
duct  from  1  to  3. 

The  dissipated  mechanical  energy  is  converted  into  heat,  and 
therefore 


where  dQ  /dx  is  the  heat  released  per  unit  time  in  the  examined 
A 

volume . 


Let  us  examine  flow  of  a  fluid  along  rectangular  and  round 
channels.  We  direct  coordinate  z  along  the  channel,  while  coordi¬ 
nates  x  and  y  are  in  its  cross-sectional  plane .  With  steady  motion 
in  a  channel  of  constant  cross  section,  the  velocity  will  be  a 
function  only  of  coordinates  x  and.y.  In  particular,  in  a  rectangular 
channel,  with  small  height/width  ratio 

2h{b«  1.  (3) 


where  h  is  one-half  the  height  and  b  is  the  width  of  the  channel 
section,  we  can  consider  that  the  velocity  changes  only  with  the 
height  of  the  section,  disregarding  the  derivative  of  velocity  with 
respect  to  coordinate  x.  In  this  case  [1]  the  velocity  will  be  a 
parabolic  function  of  coordinate  y: 


(4) 


where  dp/dz  is  the  change  in  pressure  per  unit  flow  length. 


As  fluid  moves  along  a  cylindrical  tube  of  radius  R,  in  the 
case  of  Poiseuille  flow  we  have  the  familiar  analogous  formula 


(5) 


Let-  us  average  the  velocity  across  the  channel.  For  a  rectan¬ 
gular  se  ;ion 
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for  a  round  section 
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(6) 
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Considering  (6)  and  (7),  formulas  (*J)  and  (5)  assume  the  form 


(7) 


.i  "Let  us  ^determine  the  quantity '"of  heat  released  during  dissipa-  . 
tion  of  mechanical  energy  per  unit  length  of  the  channel.  For  this' 
let  us  integrate  equation  (1)  across  the  section,  with  consideration 
of  (8)  for  a.  rectangular  section  dnd  (9)  for  a  round  one.  Let  us 

designate  by  dq  /dr  the  quantity  of  heat  released  per  unit  time  per 

»■ '  *  A 

unit  channel  ength.  Thqn,  for  the  rectangular  section 

'  ■  ' 

»  h 


(10) 


dx 


(11) 


Formulas  (10)  and  (11)  show  that  the  heat-release  power  per 
unit  channel  length  is  proportional  to  the  square  of  the  average 
velocity;  for  a  round  section  it  does  riot  depend  on  the  channel 
radius.  , 

As  a  fluid  mcives  .along  narrow  channels,  under  conditions  of 
insignificant  relative  intensity  of  heat  transfer  between  the  fluid 
and  the  channel  walls,  the  temperature  differential  along  the  flow 
section  can  be  disregarded  as  a  small  value  compared  with  the  tem¬ 
perature  head  on  the  surface  (Bi  <<  1).  In  this  case  the  thermal 
state  of  the  flow  can  be  characterized  by  uhe  temperature  which  de¬ 
pends  only  on  one  coordinate,  z. 


Let  us  assume  that  heat  transfer  between  the  fluid  and  the 
channel  walls  is  characterized  by  neat-transfer  coefficient  a.  Then 
the  change  in  heat  content  of  a  flow  element  dz  is  expressed  by  the 
algebraic  sum  of  the  heat  of  dissipation  and  the  exchange  heat  be¬ 
tween  the  channel  and  the  fluid.  The  corresponding  heat-balance 
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equation 
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(12) 


where  p,  c  are  the  density  and  specific  heat  cf  the  fluid;  S  and  P 
are  the  area  of  the  section  and  the  perimeter  of  the  channel;  tc  is 
the  characteristic  temperature  of  the  medium. 


After  transformation  we  get 


±  _L  =,0. 

dx  qX  dx  QfS 


(13) 


The  differential  equation  compiled  characterizes  the  change, 
with  time,  of  the  temperature  of  an  isolated  flow  element.  With 
stationary  flow  we  can  turn  to  an  examination  of  the  temperature 
field  of  the  entire  flow  by  exchange  of  variables  from  the  relation¬ 
ship 


ds 

uT  ‘  — ■  • 


(14) 


Substituting  (14)  into  (13)  and  separating  the  variables,  let 
us  integrate  the  obtained  expression: 


dt 


c  F  l  <tg-  ' 

v  ecS  c  ucS  dx 


( 15 ) 


or,  in  dimensionless  coordinates. 
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dt 
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(16) 
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where  Z  i3  the  channel  length;  t*  «  Z/v  is  the  characteristic  time 
for  the  total  volume  of  fluid  in  the  channel. 


Precise  integration  of  equation  (15)  is  possible  if  we  know  the 
dependences  of  fluid  density,  specific  heat,  and  viscosity  on  the 
temperature . 


To  explain  the  influence  of  the  heat  of  dissipation  on  the 
nature  of  the  temperature  field  of  the  flow,  in  first  approximation 
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U^mv^.v1  .?  *  const,  t}5**  const.  (IT.) 

These  physical  values  can  be  replaced  by  their  average ‘ values , 
if  the  total  change  in  temperature  as  the  fluid  moves  along  the 
channel  is  insignificant  and  is  not  accompanied  by  phase  transitions 
In  this  case.  Integration  of  equation  (15)  gives 
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*  Let  us  expand  the  right  side  in  a  series,  restricted  to  the 
linear  approximation: 


dx  _  j _ t_  aF 
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After  transformations  we  get 


Equation  (20)  solves  the  problem  in  the  approximation  used. 
For  a  rectangular  channel,  considering  (10), 

,  x  r  2a  (2ft  -*•  t\  “ 


For  a  round  channel,  using  (11),  we  get 


From  equations  (21)  and  (22)  we  can  determine  the  conditions 
under  which  dissipation  becomes  substantial.  Obviously,  in  thi3 
case  the  terms  in  brackets  should  be  of  the  same  order,  for  a  roct 
angular  section  this  condition  (with  b  >>  h)  will  be 


v*  A 


from  which 


(23) 


(24) 


Besides  conditions  (23)  and  (24),  which  determine  the  order  cf 
the  flow  velocities  v  for  which  the  dissipative  processes  become 
substantial  under  the  given  conditions,  we  must  take  into  account 
the  limit  of  applicability  of  formulas  (21)  and  (22),  associated 
with  retention  of  flow  laminarity,  i.e., 

Res-SRc  .  (25) 

Condition  (25)  determines  the  upper  boundary  of  velocity  v,  up 
to  which  the  obtained  dependences  remain  valid. 


If  the  flow  velocities  exceed  the  critical  values,  corresponding 
to  condition  (25),  dependences  (21)  and  (22)  become  invalid,  since 
the  heat  of  dissipation  dqfl/dx  will  no  longer  be  defined  by  formulas 
(10)  and  (11),  which  are  quadratic  functions  of  the  average  velocity. 
In  a  turbulent  regime  we  determine  the  heat  of  dissipation  from  the 
conditions  of  conservation  of  energy,  considering  that  the  work  of 
external  forces  is  completely  converted  into  heat: 

dq^dpS,  (25*) 

where  dp  and  dqfl  are  the  drop  in  pressure  and,  correspondingly,  tf . 
average  release  of  dissipative  heat  in  channel  section  dz. 


* 


* 


From  known  hydraulics  formulas,  the  drop  in  pressure  is  deter¬ 
mined  from  the  relationship 


dp  — 


(26) 


where  R_  =  S/F  is  the  hydraulic  radius;  X  =  f(Re)  is  the  friction 
resistance  coefficient. 


Substituting  (25)  into  (25')  and  dividing  both  sides  of  the 
equation  by  dx,  we  get 


The  expression  obtained  for  the  dissipative  function  is  ?ub- 
stituted  into  equation  (20): 


Equation  (28)  defines  the  temperature  field  of  the  flow  under 
conditions  of  turbulent  flow.  Here,  consideration  of  the  dissipative 
heat  is  necessary  with  flow  velocities  defined  by  the  relationship 
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or 


(30) 


The  conclusions  arrived  at  are  valid  if  the  physical  parameters 
(p,  c,  n)  are  constant,  which  occurs  under  flow  conditions  close  to 
adiabatic  (a  *  0)  and  with  low  temperature  heads  on  the  contact  sur¬ 
face  between  the  fluid  and  the  channel  walls  (tQ  -  tQ  *  0).  Let  us 
determine,  under  these  conditions,  the  influence  of  fluid  viscosity 
on  the  nature  of  heat  transfer  with  laminar  flow.  Examining  a  round 
channel  section,  let  us  transform  somewhat  equation  (22),  expressing 
velocity  v  in  terms  cf  viscosity  n  [see  formula  (7)]  and  the  pressure 
gradient  dp/dz: 


Hiflut)  , 


From  equation  (31)  it  follows  that  a  change  in  flow  temperature 
under  adiabatic  conditions  (with  a  »  0)  does  not  depend  on  the  vis¬ 
cosity  of  the  flowing  fluid: 


Equation  (32)  could  have  been  obtained  directly  from  the  law 
of  conservation  of  energy  for  a  moving  fluid  with  constant  viscosity 
in  the  absence  of  heat  transfer  with  the  ambient  medium. 


If  there  is  heat  transfer  between  the  fluid  and  the  channel 
walls,  as  follows  from  equation  (31)  the  viscosity  increases  its 
role  in  the  total  energy  balance,  and  an  increase  in  n  with  constant 
moving  flow  forces  (dp/dz  *  Ap /l  *  const)  leads  to  a  relative  re¬ 
duction  in  the  role  of  the  dissipative  processes. 


We  should  once  more  stress  that  these  conclusions  pertain  to 
small  a  and  small  initial  differences  in  temperatures  t  -  tQ.  Under 
conditions  of  intense  heat  transfer,  equation  (20)  is  invalid,  while 
with  large  initial  differences  in  temperatures  between  the  fluid  and 
the  channel  walls  equation  (19)  is  also  invalid,  since  the  temperature 
changes  of  density,  specific  heat,  and  viscosity  of  the  fluid  begin 
to  play  an  essential  role.  For  drop  fluids,  of  the  indicated  physical 
parameters  it  is  the  visco^^fy  which  changes  particularly  greatly. 
Consideration  of  the  temperature  dependence  of  viscosity  leads  to 
considerable  complications  in  integrating  equation  (15),  while  the 
average  velocity  can  no  longer  be  calculated  by  formulas  (6)  and  (7), 
which  result  from  solution  of  the  hydrodynamics  equations  for  fluids 
with  constant  viscosity.  Since  the  flow  temperature  is  a  function  of 
coordinate  z,  the  viscosity  must  be  treated  as  a  function  of  this 
same  argument.  The  pressure  gradient  along  the  flow  will  also  be 
variable.  Therefore,  equations  (6)  and  (7)  must  be  examined  as  dif¬ 
ferential  equations,  while  the  average  velocity  Is  determined  by 
integrating  them  along  the  channel.  For  a  rectangular  channel,  from 
equation  (6)  we  get 


from  which 
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For  a  round  section 
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For  a  rectangular  section  the  rate  of  release  of  the  heat  of  dis¬ 
sipation,  defined  by  equations  (10)  and  (11)  with  consideration  of 


For  a  round’  section 
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(36) 


Equation  (13),  which  describes  the  temperature  field  of  a  fluid 
flow,  with  consideration  of  (14),  (35),  and  (36)  assumes  the  form 
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For  a  rectangular  section  the  coefficient 

A- 3/A1,  (38) 

for  a  round  section 

A«8//?*.  (39) 

If  the  temperature  dependence  of  viscosity  n(t)  is  known,  the 

J 

integral  \  \)dz  can  be  calculated.  For  this  we  must  know  the  tempera- 
o 

ture  field  of  the  flow  z(t).  Then,  replacing  the  variables  under 
the  definite  integral,  we  get 

‘  'i 

(40) 

o  1  !i 

where  t^  is  the  temperature  of  the  flow  at  the  channel  outlet. 

Substitution  of  (40)  into  (37)  results  in  a  nonlinear  integro- 
differential  equation  for  determining  the  dependence  z(t).  To  ob¬ 
tain  an  approximate  solution  to  equation  (37)  let  us  assume  that  the 
density  arid  specific  heat  of  the  fluid  are  constant  and  equal  to 
their  average  values  in  the  examined  temperature  range.  For  drop 
fluids  p  and  c  change  considerably  less  than  viscosity  n  with  a 
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change  in  temperature,  and  therefore  euch  an  assumption  is  Justified. 
Further,  using  the  theorem  of  the  mean,  let  us  represent  integral 
(*J0)  in  the  form 

ftl'/z-rff,  (Hi) 

where  n*  Is  an  intermediate  value  for  viscosity  in  the  interval 
Cn0>  n j ] ;  n0  and  are  the  vis oosi ties  of  the  fluid  at  the  channel 
inlet  and  outlet,  respectively. 


Let  us  rewrite  equation  (37)  considering  (Hi): 

dz - S*£i - . 


(••2) 


Dimensionless  viscosity  n/n*  is  a  temperature  function,  equal 
to  one  within  the  temperature  range  [tg,  t^].  In  linear  approxima¬ 


tion 


»i  \t  . 


(^3) 


coefficient  h  is  selected  such  that  at  point  t  ■  t#  linear  function 
(H3)  has  a  general  derivative  with  function  n/n*  (see  the  figure): 


<**) 


Substituting  expression  (H3)  into  (H2)  we  get 


t/z  - 


c  cut 


(15) 


Integrating  the  obtained  equation,  we 
arrive  at  the  following  expression  for  the  flow  temperature  field: 


Tu/  ir  •"/• 


J  #p 


<1 


-  *  ~] 


AJ 


exp  (  —c 


£  I  * 

-  - -  ~~T  ■ 

S  Sp  I  I  ‘ 


(K6) 


When  z  *  l  expression  (H6)  is  the  equation  for  deterndning  the 
temperature  of  the  fluid  at  the  channel  outlet  (t^).  However,  so¬ 
lution  of  this  equation  is  possible  if  we  know  the  values  of  n*  and 


n.  These  parameters  can  be  determined  by  the  method  of  suooessive 
approximations.  In  first  approximation  we  oan  set  nj  ■  n0  and  - 
-  -X.  After  this,  from  equation  (46)  we  determine  t ^  and  ds/dt} 
then,  using  equations  (40)  and  (41),  knowing  the  temperature  depen¬ 
dence  of  viscosity  n(t) ,  we  find  the  second  approximation  for  n 
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Solving  the  equation  n  ■  n(t)  for  l  (or  graphically),  let  us 
determine  the  temperature  tJJ  corresponding  to*viscoaity  Sub¬ 

stituting  this  temperature  into  equation  (M3),  we  find  the  second 
approximation  for  coefficient  h: 

.  i'n  • 

( M8 ) 

1  ,/r 

The  values  of  and  h2  are  substituted  into  equation  (46)  and 
when  z  *  l  we  find  the  second  approximation  for  the  flow  temperature 
at  the  channel  outlet.  t«2.  To  determine  the  third  approximation, 
the  entire  calculation  process  is  repeated. 


The  method  described  allows  us  to  construct  the  temperature  field 
of  a  flow  of  fluid  as  it  moves  in  narrow  channels,  with  consideration 
of  the  temperature  dependence  of  viscosity  and  the  dissipation  energy* 
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